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Abstract 



Modified Gravity 
and the phantom of dark matter 
by 

Joel Richard Brownstein 



Astrophysical data analysis of the weak-field predictions support the claim that modi- 
fied gravity (MOG) theories provide a self-consistent, scale-invariant, universal description 
of galaxy rotation curves, without the need of non-baryonic dark matter. Comparison to 
the predictions of Milgrom's modified dynamics (MOND) provide a best-fit and experimen- 
tally determined universal value of the MOND acceleration parameter. The predictions of 
the modified gravity theories are compared to the predictions of cold non-baryonic dark 
matter (CDM), including a constant density core-modified fitting formula, which produces 
excellent fits to galaxy rotation curves including the low surface brightness and dwarf 
galaxies. 

Upon analysing the mass profiles of clusters of galaxies inferred from X-ray luminosity 
measurements, from the smallest nearby clusters to the largest of the clusters of galaxies, it 
is shown that while MOG provides consistent fits, MOND does not fit the observed shape 
of cluster mass profiles for any value of the MOND acceleration parameter. Comparison 
to the predictions of CDM confirm that whereas the Navarro-Frenk- White (NFW) fitting 
formula does not fit the observed shape of galaxy cluster mass profiles, the core-modified 
dark matter fitting formula provides excellent best-fits, supporting the hypothesis that 
baryons are dynamically important in the distribution of dark matter halos. 
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Chapter 1 



Introduction 

"Whoever undertakes to set himself up as a judge of Truth and Knowledge is 
shipwrecked by the laughter of the gods. " 

Albert Einstein 

As Einstein observed, one's perspective into the nature of physical reality determines the 
degree to which one can understand physics within existing notions, or whether a shift of 
paradigm is needed. 

1.1 Motivation and Objectives 

Motivated by ongoing advances in gravity theory, the basic objective of this thesis is to 
present the current state of the art in astrophysical models of the astronomically observed 
data and to explore the frontier of cosmology by presenting a model universe without the 
necessity of dominant non-baryonic dark matter. 

The central computation - in the question of dark matter - is the spatial distribution 
of the unseen component. Application of the Newtonian 1/r 2 gravitational force law in- 
evitably points to dark matter halos dominating disks of visible baryons within galaxies 
and clusters of galaxies throughout the cosmos (Oort, 1932). 

Alternatively, a modified gravity field theory, which is sourced only by baryons of the 
standard model of particle physics, may provide the solution to the unseen component. 
Regardless of how minuscule the physical effect must be at planetary and solar scales, 
the inclusion requires a fundamental modification of the known interactions, resulting in 
additional configuration variables and new couplings. The addition of scalar- vector-tensor 
fields and their couplings to the action for pure gravity will be identified as a measurable 
fifth force, that does not vanish in a local frame, but without necessarily violating the 
weak equivalence principle and the universality of free fall. Although each modified gravity 
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theory may generate a fifth force through dissimilar means, the necessary violation of the 
equivalence principle may be quantified by measuring the dynamical mass factor as the 
ratio of the dynamic mass to the observed baryonic mass in each of the theories, depending 
on its distance from the astrophysical center of the system. 

Motivations for each perspective are presented in greater detail, according to whether 
there exists physical dark matter candidates corresponding to as yet undetected massive 
particle fields, as in §1.1.1, or the hypothesis that modified gravity is responsible for the 
phantom of dark matter, as in §1.1.2, according to some violation of the strong equivalence, 
or relativity principle. The full hypothesis is detailed in §1.1.3, defining the candidate 
theories and the scope of astronomical case studies. The consequences for astrophysics 
and cosmology are identified in §1.1.4 and §1.1.5, respectively. 

Whereas each of the main ideas are theoretically developed in Part I, the specific 
objective of this thesis, documented in Part II, is to provide the results of a computational 
astrophysical data analysis, across a range of scales, in order to determine whether dark 
matter is real, or a phantom of a modified gravity theory, according to the method: 

1. To measure the visible baryon distributions in galaxies and clusters of galaxies within 
Newton's theory, Milgrom's modified dynamics (MOND) and Moffat's modified grav- 
ity theories (MOG), and to measure the dark matter distributions in a sample of 
galaxies and clusters of galaxies. 

2. To compute the dynamic mass factors in each of the sample's galaxies and clusters 
of galaxies, which provide the extra gravity and phantom dark matter in MOND and 
MOG; and to compare the dark matter factor with cosmological values. 

3. To observe whether the results may be universally understood across the range of 
kiloparsec to megaparsec scales within the sample. 

Conclusions are supplied in Part III. 

1.1.1 Non-Baryonic Dark Matter 

Despite over 20 years of focussed experimental effort, no direct evidence of dark matter 
particles has ever been found, and no annihilation radiation from any non-baryonic, dark 
matter candidate has ever been detected. In fact, no experiment has ever supported any 
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physics beyond the standard model of particle physics on which the dark matter hypothesis 
depends. One of the goals of CERN's new Large Hadron Collider (LHC) is to understand 
what constitutes dark matter. Capturing the notion that dark matter and dark energy 
are emergent gravitational phenomena due to terms beyond the Einstein-Hilbert action 
is a primary objective for this thesis. Just as the Michelson-Morley experiment gives a 
null result and falsifies the prediction of the luminiferous aether - the medium for the 
propagation of light as it was thought until the late 19th century - the LHC and the dark 
matter project (and future experiments) may never confirm any candidate for the dark 
matter particle. The scientific community was strongly resistant to accept the falsification 
of the luminiferous aether, and experiments continue to this day in search of this classically 
motivated, but expendable substance. The search for dark matter continues, as Kipling's 
poem IF reasons, regardless. 

The most popular candidates for non-baryonic dark matter, and the experimental status 
of their respective searches, are listed in Table 1.1. In all cases, a modification to the 
standard model of particle physics is required to explain why the A-CDM predicted density 
of dark matter dominates the baryon component by a factor of 5.7 ± 0.4 (Spergel et al, 
2007). 

The status of dark matter is entirely reversed in the astrophysics of galaxies and clusters 
of galaxies, where the baryon fraction in certain critical systems can be entirely neglected, 
and has led some astrophysicists to claim detection of dark galaxies devoid of stars (Minchin 
et al, 2007). The search for dark dwarf galaxies is an important prediction of the A-CDM 
cosmological model of structure formation. However, after years of intense search, all of the 
candidate dark dwarf galaxies have turned out to have been dwarf galaxies with a luminous 
stellar component, observed with high-powered optical telescopes. In fact, for decades, dark 
matter simulations have predicted many times more companion galaxies than are actually 
observed (Merrifield, 2005). Large galaxies like the Milky Way are accompanied by a 
local group including many small satellites. However, a robust prediction of dark matter 
simulations is that the Milky Way should be accompanied by ~ 500 satellites, whereas only 
35 have been observed (Moore et al : 1999a). Klypin et al (1999) used the circular velocity 
distribution of the galaxy satellites to conclude that unless a large fraction of the Local 
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Table 1.1: Standard model extensions for non-baryonic dark matter 

Some of the primary candidates for non-baryonic dark matter are 
listed (Munoz, 2004), (Taoso et a/., 2008). 



Particle 



Conjecture 



Neutral scalars associated with the spontaneous symmetry 
breaking of global £7(1) Peccei-Quinn symmetry as a mech- 
anism to solve the strong CP problem in QCD. These light 
- 10 /xeV - non-thermally produced species are viable CDM 
candidates with relic abundances matching the cosmological 
density. Reviewed in Sikivie (2005) and Raffelt (2007). 



Axion 



Weakly interacting massive particles whose annihilation rate 
fell below the cosmological expansion rate, freezing out at 
primordial density, are thermal relics. The required small 
annihilation cross-section introduces new physics at the 
weak scale. 



WIMP relics 



The lightest supersymmetric WIMP may be stable and may 
act as a heavy thermal relic. The near collisionless aspect 
of the neutralino makes it a prototype for cold dark matter. 
Supersymmetry is badly broken in nature, and requires new 
physics at the weak scale. 



Neutralino 



The first proposed dark matter candidate was a fourth gen- 
Neutrino eration heavy neutrino, but collider experiments have ruled 
this out up to 1 TeV. 

The supersymmetric partner to the neutrino has been ruled 
out by LEP up to 1 TeV, which is too massive to be the 

■at , • lightest stable supersymmetric particle. The left-handed 

sNeutrmo ° . % J . ^ . ... 

sNeutrmo has further been excluded m the minimal super- 
symmetric standard model (MSSM), but the right-handed 
sNeutrino is viable in extensions to the MSSM. 
The supersymmetric partner to the graviton may have a 
mass of order keV, and is a candidate for non-baryonic dark 
matter in the absence of the inflationary paradigm. How- 
ever, relic gravitinos at big bang nucleosynthesis would lead 
to a higher reheating temperature than permitted by ther- 
mal leptogenesis. 



Gravitino 
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Group satellites has been missed in observations, there is a dramatic discrepancy between 
observations and hierarchical models, regardless of the model parameters. Furthermore, 
Somerville et al (2004) showed that the A-CDM hierarchical scenario predicts that the 
largest most massive galaxies should form last, yet high redshift observations are beginning 
to indicate that some fraction of very luminous galaxies were present quite early in the 
process of structure formation. 

The dark matter paradigm has spawned an entire industry of computer simulations 
which attempts to model the formation of structure in the universe through the gravita- 
tional collapse of dark matter dominated clumps. As the resolution of very large computer 
simulations continues to improve, definite and robust results for the dark matter distribu- 
tion within individual galaxies show a central power-law cusp, p(r) oc r -7 . The steepness 
of the cusp is a topic of debate (Navarro et al, 2004), particularly in regard to the inclusion 
of baryons in the simulations. 

Bullock et al (2001) considered the angular momentum of dark matter halos, and pro- 
vided high-resolution N-body simulations of the ACDM cosmological structure formation 
of galactic halos, showing that the HI (and He) gas components cool at early times into 
small mass halos, leading to massive low- angular momentum cores in conflict with the 
observed exponential disks. The simulated ACDM galaxies have profiles which are too 
dense at small radii, and with tails extending too far. A possible solution is to associate 
the central excesses with bulge components and the outer regions with extended gaseous 
disks. 

In the cores of spiral galaxies that have been fitted, the power-law index, 1 < 7 < 1.5, 
and the general consensus is that at large radii, the profile steepens approaching a power- 
law index, 7 = 3. At intermediate distance, the fits reproduce the observed flat rotation 
curves with 7 = 2. However, the large amount of luminous stellar material in the core of 
spiral galaxies means that such cusps will have negligible effects on the rotation curves of 
spiral galaxies. Fortunately, a class of galaxies exists in which this issue is not the case. 
These low surface brightness galaxies contain a very low density of luminous material, even 
in the core, so that the observed dynamics should be dominated by the gravitational forces 
of the dark halo. Unlike in high surface brightness spiral galaxies, the small amount of 
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luminous matter should not be efficient in redistributing the dark matter, so the central 
cusp should remain, de Blok et al (2001b) showed that the dark matter fits to a sample 
of 30 low surface brightness galaxies do not fit the data, showing systematic deviations in 
the galaxy cores. Although the low surface galaxy rotation curves do tend to flatten off to 
the constant rotation velocity characteristic of the dark matter halo, they rise significantly 
more slowly out to several kiloparsecs than the best-fit cold dark matter prediction. The 
rotation curves of low surface brightness galaxies prefer dark matter distributions with 
constant density modified cores. 

Just as the luminiferous aether is considered an unnecessary addition to physics that 
violates the principle of Occam's razor, non-baryonic dark matter may be considered a 
superfluous component to be removed from our description of the cosmos. The application 
of Occam's razor here would be to decrease the computed masses of galaxies and clusters 
of galaxies by an order of magnitude, and several orders of magnitude in the case of the 
critical systems. Although the evidence we currently have for dark matter is concordant, 
it is unconfirmed, so it is vital to investigate the alternatives. In some cases, modified 
gravity may seem to fit the observations better than dark matter. 

The conflict between the cuspy dark matter halos predicted by hydrodynamical N-body 
simulations and the constant density cores preferred by dwarf and low surface brightness 
galaxies made it impossible to x 2 ~fit some of the galaxies in Part II with the NFW fitting 
formula of Navarro et al (1996, 1997) with a nonvanishing stellar mass-to-light ratio. 
Regardless, a cure to the cusp problem was found in Chapter 4 by implementing a core- 
modified fitting formula of Equation (2.12): 

p(r) = -rz^ (i-i) 

r 6 + r 6 s 

which provided excellent best-fits across the sample of high and low surface brightness 
galaxies including dwarf galaxies; and supports the existing A-CDM cosmology based on 
logarithmically divergent dark matter halos. For each rotation curve in the Ursa Major 
sample of §4.2, the power-law index was found to asymptotically rise (7 — ► 3) at large 
radii, whereas at small radii, the baryon-dominated constant density dark matter core 
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(7 — > 0) appeared, as shown in Figure 4.5. Including the baryons, the total power-law 
indices oscillate around 7 —> 1 at the small radii and asymptotically rise to 7 —> 3 at large 
radii, where the dark matter core dominates and the baryonic data runs out. 

Although cuspy dark matter halos are robust predictions from hydrodynamical N-body 
simulations of clusters of galaxies which ignore the baryon components, it is a feature which 
has not been observed. Whereas it was impossible to x 2_ fit any of the clusters of galaxies 
in Part II with the NFW fitting formula of Navarro et al. (1996, 1997), the core-modified 
profile of Equations (1.1) and (2.12) that solved the cusp problem for galaxy rotation curves 
was successfully applied to clusters of galaxies in Chapter 5, providing excellent best-fits 
across the sample from the smallest X-ray cluster in nearby Virgo to the largest and most 
radiant Bullet Cluster 1E0657-558 main cluster. The core-modified dark matter halos 
described in §5.2.1 do not dispute the existing A-CDM cosmology based on logarithmically 
divergent dark matter halos, but challenge dark matter simulations to reinsert the baryons, 
and make predictions that are consistent with astrophysical observations. 

1.1.2 Modified gravity 

Aristotle's notion of the motion of bodies was that a constant force maintains a body in 
uniform motion, and that force could only be applied by contact so action at a distance 
was considered impossible. Aristotle produced a number of arguments why the heavens 
revolved around the Earth, and denied the possibility that the Earth rotated on its axis. 
This geocentric model stood the test of time for over eighteen hundred years until Coper- 
nicus postulated that the sun was the center of the universe and the earth revolved around 
it at a distance related to the size of the orbit. This idea, although controversial, initi- 
ated a scientific revolution that allowed a mathematical description of the force of gravity. 
Kepler's empirical three laws of motion determined the elliptical orbits of planets orbiting 
the sun and opened an era of precision astronomy. Galileo continued the revolution with 
a series of experiments on projectile motion, including the legendary Tower of Pisa ex- 
periment on the universality of free fall, and developed the mathematical theory of falling 
bodies. Newton (1687) set the foundation for classical mechanics and introduced the law 
of universal gravitation and a derivation of Kepler's laws of planetary motion. Newton's 
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gravity theory describes action at a distance through Poisson's equation, in which any 
change in the matter distribution is instantaneously communicated to bodies in motion 
through the gravitational potential field. 

Scalar modified gravity 

The idea that Newton's constant varies from one point in space-time to another was first 
considered by Jordan (1959) and first implemented in Brans and Dicke (1961), which 
attempted to modify general relativity to be compatible with Mach's principle. A non- 
geometric, scalar field, 0(x), was coupled to the Ricci curvature scalar in the action for 
gravity. This field is massless, but self-interacting as it couples to its own kinetic term 
through the coupling constant, ujed, which must be determined experimentally. This self- 
consistent and energetically stable modification to general relativity leads to the result 
that the locally measured value of Newton's constant varies spatially, and depends on 
the expectation value of the inverse of the Jordan-Brans-Dicke field, G{x) = (l/0(x)). 
However, this scalar-tensor modification to general relativity violates the strong equivalence 
principle and leads to a variation in Kepler's third law - which is locally measurable through 
the 7 parameter of the parametrized post-Newtonian formalism. Solar system tests and to 
a lesser degree, PSR1913+16, have constrained Jordan-Brans-Dicke theory (Will, 2006). 

Renormalized gravity 

Even without the presence of auxiliary scalar fields which carry energy and momentum, the 
evidence from quantum field theory suggests that coupling constants like the fine structure 
constant, a ) in quantum electrodynamics (QED) are not really universal, but are scale 
dependent "running" quantities. The value measured in the lab for the "fundamental" 
charge of the electron depends on the renormalization scale, k. The physical mechanism 
behind the running fine structure constant is the appearance of a sea of virtual electron- 
positron pairs which are in a constant state of creation/annihilation. Offshell photons 
surround the test charge, and contribute to the polarization of the bare charge, screening it 
at large distances. Experimental measurements of the fine structure constant, at sufficiently 
large renormalization scale k ) show values a(k) 3> Just as the Coulomb force law in 
classical electromagnetism is modified by the quantum corrected Ueling potential, the 
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Newtonian force law may be subject to a running Newton's constant, leading to quantum 
corrected modified gravity (Reuter, 1998; Reuter, 2000). 

1.1.3 Hypothesis 

The solution to the missing mass problem in galaxy rotation curves and clusters of galaxies 
may be one of the candidates: 

1. Cold non-baryonic dark matter (CDM), 

2. Milgrom's modified Newtonian dynamics (MOND), 

3. Moffat's metric skew-tensor gravity theory (MSTG), 

4. Moffat's scalar tensor vector gravity theory (STVG). 

For the case of CDM, the total mass interior to a spherical region divided by the 
integrated baryonic mass of the combined, visible, components is a measure of how much 
dark matter is required - this is the dark matter factor. For the case of MOND, the dynamic 
mass interior to a region divided by the integrated baryonic mass on the same region is a 
measure of how much MOND is required - this dynamic mass factor is the inverse of the 
MOND interpolating function, /i. For the case of Moffat's modified gravity (MOG), the 
dynamic mass factor is the ratio of a running Newtonian constant, G{r) divided by the 
bare Newtonian constant, Gn = 6.67428 x 10 -11 m 3 kg _1 s 2 . 

The hypothesis will be tested, performing a best-fit to the data for each of the objects 
in Table 1.2 and measuring the dynamic mass factor for each of the candidate theories, 
completely determining the predictions that are distinct from Newton's theory. The radial 
profile of these dynamic mass factors across the Ursa Major sample of §4.2 are shown 
in Figure 4.4. The radial profile of the stellar mass-to-light ratios, for the same sample, 
are shown in Figure 4.6, providing the means for the best-fits to the galaxy rotation 
curves of Chapter 4. The high resolution properties of the solutions are used to study 
orphan features in §4.3.2, and the integrated results are used to study the theoretical Tully- 
Fisher relations, shown in Figure 4.9. In Chapter 5, the radial profile of these dynamic 
mass factors across a sample of best-fit clusters models, of §5.2, is shown in Figure 5.2, 
for each theory. Consistency with solar system experimental constraints are examined in 
Chapter 6, using data from §6.1 on the Pioneer 10/11 Anomaly, to set an upper bound on 
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dynamic mass factors, plotted in Figure 6.2. The planetary predictions and observational 



Table 1.2: Catalogue of astronomical case studies 



Case study 



Astronomical objects 



Chapter 4 Galaxy rotation curves 
High surface brightness galaxies 



Low surface brightness galaxies 



§4.2 Ursa Major filament of galaxies 

NGC 3726, NGC 3769, NGC 3877, NGC 3893, 

NGC 3949, NGC 3953, NGC 3972, NGC 3992, 

NGC 4013, NGC 4051, NGC 4085, NGC 4088, 

NGC 4100, NGC 4138, NGC 4157, NGC 4217, 

NGC 4389, UGC 6399, UGC 6973 

NGC 3917, NGC 4010, NGC 4183, UGC 6446, 

UGC 6667, UGC 6818, UGC 6917, UGC 6923, 

UGC 6983, UGC 7089 



Chapter 5 Clusters of galaxies 



§5.2 Best-fit cluster models 

Bullet 1E0657-558, Abell 2142, Coma, Abell 

2255, Perseus, Norma, Hydra-A, Centaurus, 

Abell 400, Fornax, Messier 49 

55.3 Bullet Cluster 1E0657-558 



limits are provided in Table 6.1 and the mean ephemerides of planetary orbits is provided 
in Table 6.2. Kepler's laws of motion are explored in §6.1.3, and the anomalous perihelion 
advance is studied in §6.1.5, with constraints listed in Table 6.3. 

As regards CDM halos, testing the hypothesis in Chapter 4 and Chapter 5 will 
explore the fine details of the cusp problem and too many dwarf problem, discussed in 
Chapter 2, and may indicate a consistent cure based upon the model of §2.1.3. Dynamic 
mass measurements, according to §2.2, are used to establish a fundamental connection 
to the Tully-Fisher relation, demonstrated in §4.3.4, in which the total mass vs. velocity 
relation, including dark matter, is plotted in Figure 4.8, and confirming the importance of 
including the baryons in curve-fitting dark matter halos. The core-modified dark matter 
model, as indicated in §5.2.1, also cures the cusp problem for the sample of clusters of 
galaxies studied in Chapter 5, whereas the well established NFW model of §2.1.1 does 
not allow x 2_ fits in these systems. 

As regards MOND, testing the hypothesis will provide a better measurement of the 
MOND acceleration, ag, and provide a more detailed understanding of the MOND in- 
terpolating function, /x, in the context of Milgrom's acceleration law presented in §3.2.1. 
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Dynamical mass measurements, according to §3.2.2, are used to establish the scales at 
which the theory fits observations; and where MOND potentially falls short, consider- 
ing both the possibilities of MOND without dark matter, in §5.2.2, and the possibility 
of neutrino halos in §5.3.5. Covariant theoretical foundations of MOND are explored in 



As regards Moffat's MOG theories, with MSTG presented in §3.3 and STVG presented 
in §3.3, the respective point source modified acceleration laws of Equations (3.71) and 
(3.139) are derived from an action principle in which a Yukawa fifth force combines with a 
Newtonian gravitational force, leading to modified Poisson equations, derived in §3.3.5 and 
§3.4.5, respectively. Whereas the MSTG theory has phenomenological parameters derived 
from the Tully-Fisher relation, the STVG theory has parameters which emerge as integra- 
tion constants from integrated field equations. In either case, testing the hypothesis will 
provide a better measurement of the MOG parameters. Dynamical mass measurements, 
according to §3.3.6 and §3.4.6, are used to provide support for the conjecture that the 
combination of a weak fifth force and a renormalization of Newton's constant will dynam- 
ically provide stronger gravity at astrophysical scales, providing a consistent solution to 
the missing mass problem without the addition of baryonic or non-baryonic dark matter. 

Although not specifically tested in the list of theories, Moffat's NGT is considered a 
candidate for halos of phantom dark matter, with overlapping predictions with MSTG and 
possibly STVG. It is natural that the source of the fifth force in either of the MOG theories 
is due to the fundamental NGT field excitations. 

1.1.4 Consequences for Astrophysics 

We know from general relativity that the Newtonian force is an emergent phenomenon due 
to the laws of physics in a curved space-time. Although we no longer treat gravity as a 
simple force, we do make such a simplification to perform astrophysics by means of the 
central potential, 



3.2.2. 




(1.2) 
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where Gn = 6.67428(67) x 10 -11 m 3 kg _1 s -2 is Newton's constant measured experimen- 
tally 1 . Newton's central potential, where $ is given by Equation (1.2), is an unshakeable 
foundation of modern physics, and appears upon identification of the constant of integra- 
tion in the static spherically symmetric Schwarzschild solution to general relativity, 

ds 2 = (l + ^) dt 2 - (l + ^) 1 (dx 2 + dy 2 + dz 2 ) . (1.3) 

Astrophysics assumes that the measured velocity dispersions and temperature profiles are 
determined to a good approximation by the Newtonian acceleration, 

a(t,x) = -V$(t,x), (1.4) 

which neglects relativistic effects. However, for realistic distributions of matter in galaxies, 
we have neither analytic, nor numerical solutions to general relativity from which orbits can 
be predicted. For realistic distributions of matter in clusters of galaxies, the high degree of 
symmetry improves the situation in general relativity, and we may well approximate the 
intracluster medium by the interior solution of the rotating, axially symmetric Kerr solution 
to general relativity, from which cluster masses and temperatures may be predicted with 
precision. However these interior solutions contain gravitomagnetic components due to the 
rotational energy of the system, which are not easily measured astronomically. Newton's 
universal law of motion will emerge from these other solutions with relativistic corrections. 
Unlike the simplest case of the Schwarzschild metric, the familiar distance-squared law may 
be modified, in effect. Any such modifications would have impact on the computations of 
the mass-to-light ratio in galaxy rotation curves, the temperature to mass relationship in 
X-ray clusters of galaxies, and gravitational lensing in galaxies and clusters of galaxies. 

1.1.5 Cosmological Models 

When matter sources are dominated by radiation, as in the early universe, the formation of 
local inhomogeneities is suppressed and the universe expands homogeneously and isotrop- 
ically, without the growth of structures such as stars, galaxies or clusters of galaxies. In 

iNIST 2006 CODATA value. 
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the idealized case that the dominant radiation density is constant, the Newtonian central 
potential vanishes from the solution to the Einstein equations, and the universe is perfectly 
described by the Friedmann Robertson Walker solution to general relativity, 



ds 2 = dt 2 - a 2 (t) 



-^— + r 2 (d6 2 + sin 2 6d(f ) 2 ) 
1 — kr 2 



(1.5) 



where k = {0, ±1} is the curvature index and a(t) is the expansion factor. The expansion 
in time, or redshift z, is governed by nonlinear Friedmann equations of motion, 

kc 2 

d{z)2 = 0(^1' (L6) 



where 



te G N p(z) 

[ ] 3 H{zf [ n 



is the cosmological density parameter, p(z) is the mean cosmological mass density. 

H(z) = ^ (1-8) 
a{z) 

is the Hubble parameter, measured experimentally as a fraction of today's value. 

Moffat and Toth (2007) explored the cosmological consequences of MOG, and found 
that it provides, using a minimal number of parameters, good fits to the data, including 
the cosmic microwave background temperature anisotropy, the galaxy power spectrum, 
and the supernova luminosity-distance observations without the necessity of dark matter. 

As astrophysics cannot simply do away with the Newtonian central potential, cosmology 
cannot simply do away with Newton's constant, but neither does gravity theory predict 
its value. However, in both astrophysics and cosmology, it can always be arranged that 
Newton's constant, Gjv, and the mass density, p, appear together in the combination GnP 
as in Equations (1.2) and (1.7). This argument is due to dimensional reasoning, and leads 
to an ambiguity between the necessity of dark matter, p m > and a running Newton's 
constant, G > G/v, or the existence of a MOND regime, p < 1. In the case of the dark 
matter paradigm, the density of matter p m exceeds the density of baryons p^ and the 
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combination G^Pm > G^Pb- In the case of a running Newton's constant, the combination 
Gpb > G^Pb) and the visible baryon distribution provides "extra gravity" without non- 
baryonic dark matter. However, the apparent degeneracy between dark matter and a 
running Newton's constant may be broken by calculations which involve a spatial integral 
or derivative of the combination Gp. Such is the case for galaxy and cluster lensing 
experiments and cosmological models. 

1.2 Citations to published results 

Large portions of Part II Astrophysics have been published: 

Brownstein and Moffat (2006a), "Galaxy rotation curves without non-baryonic 

dark matter", Astrophys. J. 636 721-741. arXiv:astro-ph/0506370 
Brownstein and Moffat (2006b), "Galaxy cluster masses without non-baryonic 

dark matter" , Mon. Not. Roy. Astron. Soc. 367 527-540. arXiv:astro-ph/0507222 
Brownstein and Moffat (2006c), "Gravitational solution to the Pioneer 10/11 

anomaly", Class. Quant. Gravity 23 3427-3436. arXiv:gr-qc/0511230 
Brownstein and Moffat (2007), "The Bullet Cluster 1E0657-558 evidence shows 

Modified Gravity in the absence of Dark Matter" , Mon. Not. Roy. Astron. Soc. 

382 29-47. arXiv:astro-ph/0702146. Roy. Astron. Soc. Press Note 07/44 

Some sections of Chapter 3, particularly §3.1 on violations of the strong equivalence 
principle, and §3.3 on the geometric origin of a fifth force, and Chapter 6, on Solar 
system tests, are motivated from my master's thesis, which is published: 

Moffat and Brownstein (1990), "Spinning test particles and the motion of a 
gyroscope in the nonsymmetric theory of gravitation" , Phys. Rev. D41 
3111 3117. 

1.3 Organization of the thesis 

Chapter 1 is an introduction to the problem of non-baryonic dark matter in astrophysics, 
and an overview of the solution in which a modified gravity phantom of dark matter appears 
at astrophysical distances, followed by a summary of motivations and objectives in §1.1, 
with citations to published results in §1.2. 

Part I of the thesis is divided into two chapters: Dark matter halo fitting formulae are 
provided in Chapter 2, and the derivations of the modified acceleration laws are provided 
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in Chapter 3. The two dark matter profiles used in curve-fitting, including baryons, are 
the Navarro- Frenk- White (NFW) profile, described in §2.1.1, and the core-modified halo, 
derived in §2.1.3. The three modified gravity theories used in curve-fitting, using only 
baryons, are Milgrom's modified Newtonian dynamics, described in §3.2, and Moffat's 
metric skew-tensor gravity (MSTG) and scalar-tensor- vector gravity (STVG), in §3.3 and 
§3.4, respectively, which produce a finite range, Yukawa-type, fifth force (Yukawa, 1935). 

The relativistic field theoretical versions of Milgrom's modified Newtonian dynamics 
(MOND), including Bekenstein's TEVES theory, and the general family of Einstein-aether 
gravity models that may provide a weak-field MOND-like phantom of dark matter are 
documented in §3.2.3. 

Part II covers a survey of astronomical observations across a tremendous range of astro- 
physical scales. The data used in the dissertation, and range of astrophysical phenomenon 
are organized in a catalogue of astronomical case studies, in Table 1.2. The investigation 
into the available data starts with galaxy rotation curves, in Chapter 4, and is concerned 
with dynamics in the weak field in §4.1, and uses the Ursa Major filament of galaxies, in 
§4.2, as the primary experimental observations between the 1 kiloparsec to 50 kiloparsec 
range. 

The best-fitting core-modified dark matter model of Equation (4.2) provide excellent 
fits, including the dwarf galaxies, consistent with the large distance power law behaviour 
of cold collisionless non-baryonic dark matter (CDM). All of the galaxy fits include the 
best-fitting Newtonian core model of §4.3.1, provided for comparison. The theoretical 
underpinning and the experimental status of the Tully-Fisher relation are reviewed in 
§4.3.4. 

Chapter 5 continues the investigation with X-ray clusters, in §5.1, as the primary ex- 
perimental observations between the 50 kiloparsec to 1000 kiloparsec range to the largest 
range of virialized matter, which compares the observed X-ray luminosities with the tem- 
perature profiles of the best-fit isothermal gas spheres in §5.2. The Bullet Cluster 1E0657- 
558 provides a laboratory to distinguish the direct gravitational lensing evidence for CDM 
with the modified gravity solution, in §5.3. 

The search for the phantom of dark matter within the solar system in Chapter 6, at 
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ranges between 1 AU to 50 AU, is primarily concerned with the Pioneer 10/11 Anomaly 
in §6.1, and experimental bounds. 

Conclusions are presented in Part III, which supplies a summary of contributions in 
Chapter 7, and a list of some possible future astrophysical tests in Chapter 8. Lessons 
learned from CDM halos and modified gravity theories are supplied in §7.1. Specific 
conclusions on galactic astrophysics are summarized in §7.2 with future tests in §8.1, and 
specific conclusions on cluster-scale astrophysics are summarized in §7.3 with future tests 
in 88.2. 



Part I 
Theory 



Chapter 2 



Non-baryonic dark matter 

"No great discovery was ever made without a bold guess. " 
Sir Isaac Newton 



For three centuries, Newton's theory has proven to be remarkably successful, but is limited 
to weak gravitational fields. As a classical nonrenormalizable effective theory, Einstein's 
theory has proven to be remarkably successful, and together with Newton's theory, is 
believed to fully describe the measurable gravitational physics in astrophysical systems and 
cosmology. The fact that these theories predict the necessary existence of non-baryonic 
dark matter which dominates the visible matter in the universe will create a new era for 
precision astrophysics - provided the dark matter candidate is identified and experimentally 
confirmed. Otherwise a modification of gravity, as in Chapter 3, may solve the missing 
mass question, provided there are gravitational degrees of freedom in nature that are not 
captured by Newton's or Einstein's theory. 

Whereas flat cosmological models with a mixture of radiation, ordinary baryonic mat- 
ter, cold collisionless dark matter and cosmological constant (or quintessence) and a nearly 
scale-invariant adiabatic spectrum of density fluctuations provide good fits to large scale 
(3> 1 Mpc) observations, there remains a large amount of data on galactic and sub-galactic 
scales (<C 100 kpc) which may be in conflict with the ACDM halo structure - or support a 
core-modified dark matter fitting formula which retains the large scale ACDM halo struc- 
ture (Zhao, 1996). 

2.1 Dark matter halos 

Based on three rotation curves, Roberts and Rots (1973) concluded that spiral galaxies 
must be larger than indicated by the usual photometric measurements, and suggested 
the existence of an unseen massive halo beyond the last measured point - to explain the 
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slower than Keplerian decline at large radii. This view challenged the notion of a constant 
mass-to-light ratio, with radius, and suggested a mass-to-light ratio which increases with 
distance from the center. 

Ostriker et al (1974) argued that the masses of ordinary galaxies - found by assuming 
a constant mass-to-light ratio - may have been underestimated by a factor of 10; but that 
the galaxy rotation curve in the inner region provides almost no information about the 
exterior halo mass. Upon application of a Newtonian force law, 

o(r) = _ G H M„(r) {u) 

one may obtain the Newtonian dynamic mass, M^v(r), which is the mass interior to the 
sphere of radius, r, needed to support the galaxy rotation curve. Ostriker et al (1974) 
observed that although the surface luminosity profiles, L(r), do appear to be convergent, 
Mjv(r), diverges with r either weakly (logarithmic) or strongly (linear) depending on the 
method of measurement, and concluded that within local giant spiral galaxies, 

M N (r) oc r for 20 kpc < r < 500 kpc. (2.2) 

This divergent mass-to-light ratio necessitates the existence of giant halos surrounding 
ordinary galaxies of dark matter - the implied density distribution similar to isothermal 
gas spheres in the outer parts, (Begeman et al, 1991) 

P(r) = J^~ 2 , (2-3) 

r z _|_ r z 

where r c is the core radius and p is the central dark matter density. In the limit of 
small r <C r c , the isothermal sphere model approaches a constant density core. Spherically 
integrating the constant density core model of Equation (2.3) one obtains a simple formula 
for the mass of dark matter, 

M(r) = 47rp r 3 c - tan" 1 (r/r c ) J , (2.4) 
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which diverges with the behaviour of Equation (2.2), for r 3> r c . 

Einasto et al. (1974) studied the distribution of missing mass, as it relates to galactic 
morphology, concluding that the distribution is suggestive of a corona (surrounding the lu- 
minous disk), increasing the total mass of the galaxy by an order of magnitude. Rubin et al. 
(1978) considered extended rotation curves of 10 high-luminosity galaxies, and reproduced 
observed velocities using mass distributions from disk or spherical models; and suggested 
that the flat rotation velocity, v outl is not correlated with luminosity or with radius, but 
with extended dark matter. However, the observations did not suggest whether spherical 
or disk models were favoured. Rubin et al. (1978) concluded that the total mass-to-light 
ratio is higher for early-type galaxies leading to a large intrinsic scatter in the Tully-Fisher 
relation. 

On larger than galaxy scales, Fillmore and Goldreich (1984) considered the self-similar 
gravitational collapse of collisionless dark matter in a perturbed Einstein-de Sitter universe, 
and suggested that spherically averaged solutions prefer similar halo mass profiles which 
may be approximated by a power-law in the distance from center of symmetry. In the 
case of structure evolving hierarchically from a scale-free Gaussian field of a given power 
spectrum, Hoffman and Shaham (1985); Hoffman (1988) suggested that the final virialized 
halo should have an asymptotic density profile given by 



where 7 = 2 assuming that the CDM power spectrum, on galactic scales, is effectively, 



2.1.1 Navarro- Frenk- White profile 

In search of a universal description of collisionless dark matter, Navarro et al. (1996, 1997) 
provided power-law fits to halo density profiles using N-body simulations, showing that 
halo profiles are shallower than r~ 2 near the center and steeper than r~ 2 near the virial 
radius. The NFW profile is then a simple fitting formula to Equation (2.5), with a radially 



p(r) oc r 



(2.5) 



P(k) oc fc neff , where n eff = 2. 



(2.6) 
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varying powerlaw 1 < 7(7*) < 3, to describe spherically averaged density profiles: 



p(r) 



(2.7) 



r(r + r s ) 2 



Spherically integrating the NFW profile of Equation (2.7) one obtains a simple formula for 
the mass of dark matter, 



which diverges logarithmically, for r 3> r s . In the limit of small r <C r 5 , the NFW fitting 
formula of Equation (2.7) approaches the power-law with 7 — ► 1; and in the limit of large 
r ^> r s approaches the power-law with 7 ^ 3 ~ which does not approximate isothermal 
spheres. Navarro et al (1996) reported that rotation curves from galaxies ranging in size 
from giant to dwarf, satellites and gaseous atmospheres are compatible with the NFW 
halo structure of Equation (2.7) provided the mass-to-light ratio increases with luminosity. 
Navarro et al (1996) determined that the central regions of the NFW distribution have 
densities comparable to the luminous parts of galaxies. 

Although the N-body problem can easily be defined, and numerically simulated in 
the world's best computers, the problem defies any rigorous analytic treatment. Zait 
et al (2008) reported that N-body numerical simulations are providing conflicting evidence 
regarding the asymptotic behaviours of the density slope, 7 of the profile at small radii (in 
the inner region of the halo). 

2.1.2 Generalized profile 

Burkert (1995) fitted a sample of several dark matter dominated dwarf galaxies employing 
a phenomenologically modified universal fitting formula, 



which, as in the case of the isothermal sphere of Equation (2.3), approximates a constant 
density core, 7 —> at r <C r s - instead of a divergent 7 = 1 core - but otherwise agrees 




(2.8) 



p(r) 



(2.9) 



(r + r s )(r 2 + r 2 s ) ' 
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with the NFW profile, with 7 —> 3 at r >> r s . Spherically integrating the Burkert model 
of Equation (2.9), one obtains an analytic formula for the mass of dark matter, 

M(r) = 7Tp r 3 s {ln(r 2 + r 2 s ) + 21n(r + r s ) - 41n(r 5 ) - 2 tan" 1 (r/r a )} , (2.10) 

which diverges logarithmically, for r » r s . 

Zhao (1996) hypothesized that the NFW fitting formula must be broadened to account 
for the basic observed features of galactic dynamics, including less cuspy cores: 

p(r) = —t — P ° r L m > ( 2 - n ) 

' V 7 rpC^rpd _|_ rpd j [O — C) / CL 1 V 7 

where (a, 6, c) are free parameters. The NFW fitting formula of Equations (2.7) and (2.8) 
correspond to Equation (2.11) with an inner cusp with logarithmic slope c = 1, an outer 
corona with logarithmic slope 6 = 3, and a "turnover" exponent of a = 1. Syer and White 
(1998) argued that the existence of a 7 C 1 core is inconsistent with the hierarchical for- 
mation scenario of dark halos, which are much more likely to result in cuspy central density 
distributions. The least cuspy fitting formula, the isothermal spheres of Equations (2.3) 
and (2.4) correspond to Equation (2.11) with a constant density inner core with logarith- 
mic slope c = 0, an outer corona with logarithmic slope 6 = 2, and a "turnover" exponent 
of a = 2. Although Burkert's fitting formula of Equation (2.9) cannot be expressed in the 
core-modified form of Equation (2.11), it does bridge the constant density, 7^0, core 
behaviour of the isothermal sphere with the 7^3 large r behaviour of the NFW profile. 

Kravtsov et al (1998) used the rotation curves of a sample of dark matter dominated 
dwarf and low surface brightness (LSB) galaxies, employing the modified universal fitting 
formula of Equation (2.11) with a shallow cusp, (a, 6, c) = (2,3,0.2), and computed that 
a dominant fraction (~ 95%) of the dynamical mass is due to dark matter at the last 
measured point in the rotation curve; but with 0.2 < 7 < 0.4 in the inner region, r <C r S) 
of every galaxy in the sample. 

McGaugh and de Blok (1998) enforced the claim that the severity of the mass dis- 
crepancy in spiral galaxies is strongly correlated with the central surface brightness of the 
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disk. Progressively lower surface brightness galaxies have ever larger mass discrepancies. 
No other parameter (luminosity, size, velocity, morphology) is so well correlated with the 
magnitude of the mass deficit. 

Regardless of the galactic and sub-galactic data, collisionless dark matter N-body sim- 
ulations continue to predict steep inner cusps (Moore et al, 1998). Moore et al (1999b) 
argued that a universe dominated by cold dark matter fails to reproduce the rotation 
curves of dark matter dominated dwarf and LSB galaxies; and instead provided fits em- 
ploying the modified universal fitting formula of Equation (2.11) with a steeper cusp, 
(a, 6, c) = (1.5, 3, 1.5). However, these fits purposely ignored the contribution from the HI 
gas and the stellar disk to maximize the dark matter halo in the core. In contrast, the 
stellar mass-to-light ratio, T, is critical to the study of galaxy rotation curves; and the 
requirement that the stellar mass-to-light ratio, T = 0, is too strong and therefore should 
be suspect as the reason for Moore et al (1999b) good core-modified best-fits. 

Navarro (1998) remarked that a subset of spiral galaxies have flat rotation curves, and 
suggested that disagreement with the rotation curves of a few dwarf galaxies may signal 
systematic departures from the NFW shape, and that the rotation curves for LSB galaxies 
may be better described by shallower central density profiles, than presumed in the NFW 
fitting formula, van den Bosch et al (2000) argued that the spatial resolution of LSB 
rotation curves is not sufficient to put any meaningful constraints on the dark matter 
density profiles, but conceded that the rotation curves of dark matter dominated dwarf 
galaxies are inconsistent with steeply cusped dark halos. Kleyna et al (2003) demonstrated 
that the most dark matter dominated dwarf galaxies in the Local Group have constant 
density 7^0 halo cores, and suggested that CDM disagrees with observations in that 
end of the galaxy mass spectrum. Based on a series of high-resolution N-body simulations 
designed to examine whether the density profiles of dark matter halos are universal, Jing 
and Suto (2000) found that the dark matter density depends on the total halo mass, making 
it difficult to link the inner slope with the primordial index of the fluctuation spectrum. 

The unexplained behaviour of the computed dark matter distribution in the core is 
known as the cusp problem and casts doubt on the choice of the NFW fitting formula which 
presupposes the core behaviour. These discrepancies at the galactic and sub-galactic scale 
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have stimulated a number of alternative proposals. Spergel and Steinhardt (2000); Ostriker 
and Steinhardt (2003) reviewed the situation for collisionless dark matter predictions - 
overly dense cores in the centers of galaxies and clusters and an overly large number of 
halos with the Local Group compared to actual observations - and suggested the alternative 
of self -interacting dark matter produces distinctive modifications on small scales that can 
be tested through improved astronomical observations. Stoehr et al (2002) commented 
that these self-interacting dark matter modifications either may fail to reproduce the large 
observed velocity dispersions in the Local Group dwarf galaxies; or may suffer from a fine- 
tuning problem. Modifying the microscopic physics of the dark matter particles may work 
to reduce the concentration in the central regions of galaxies and to reduce the abundance 
of halo substructure (unseen dwarf galaxies) . 

Binney and Evans (2001) claimed that the Milky Way has considerably less dark matter 
in the luminous disk than expected, particularly near the galactic center, and concluded 
that cuspy halos favoured by the cold dark matter cosmology (and its variants) are in- 
consistent with the observational data. Dave et al (2001) presented a comparison of halo 
properties in cosmological simulations, confirming that collisionless dark matter yields 
cuspy halos that are too centrally concentrated, as compared to observations, de Blok 
et al (2001a) found that, at small radii, the mass density distribution is dominated by a 
nearly constant density core with a core radius of a few kiloparsecs, and found no clear 
evidence for a cuspy halo in any of the low surface brightness galaxies. S waters et al 
(2003) presented a sample of 15 dwarf and low surface brightness galaxies, showing that 
most are equally well or better explained by constant density cores, and none require halos 
with steep cusps. Gentile et al (2004) confirmed that the distribution of dark matter in 
spiral galaxies is consistent with constant density cores. 

2.1.3 Core- modified profile 

Consider a fitting formula of the form of Equation (2.11) with a constant density inner 
core with logarithmic slope c = 0, an outer corona with logarithmic slope 6 = 3, and a 
"turnover" exponent of a = 3: 



24 



Chapter 2. NON-BARYONIC DARK MATTER 



which, as in the case of the isothermal sphere of Equation (2.3) and the Burkert model 
of Equation (2.9), approximates a constant density core, 7 — ► at r <C r s - instead of a 
divergent 7 = 1 core - but otherwise agrees with the NFW profile, with 7 —> 3 at r >> r s . 
Spherically integrating this core-modified model of Equation (2.12), one obtains a new 
analytic formula for the mass of dark matter, 

M(r) = ^ Ws 3 {ln(r 3 + r 3 ) - ln(r 3 )} , (2.13) 

which diverges logarithmically, for r » r s . 

Utilizing the form of the power-law of Equation (2.5), the power-law index of the profile 
of Equation (2.12) is minus the logarithm slope 

, x dlnp(r) 3r 3 . ^ A . 

i(r) = — = - — -. 2.14 

a In r r 6 + r 6 s 

The central density, p(0) = po 5 is finite and may be written in terms of the cosmological 
critical density, p c (^), and the concentration parameter, 5 C1 

Po = Pc{z)5 C) (2.15) 

where z is the redshift. Moreover, the dark matter density at r = r s is one-half the central 
density, 

P(r s ) = \po, (2.16) 
and the power-law index of Equation (2.14) is 

7W = 3/2, (2.17) 

which is the intermediate value between the inner core with logarithmic slope 7^0, and 
outer corona with logarithmic slope 7^3. This means that the halo's constant density 
core is limited to the region r < r s , where baryons dominate the galaxy, which is important 
for N-body simulations. 
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2.2 Dynamic mass 

Theoretical predictions of dynamical quantities such as galaxy rotation curves and cluster 
masses of galaxies, as in Chapter 4 and Chapter 5, respectively, are either difference 
calculations as in the case of dark matter, or divisive ones as in the case of the modified 
gravity models of Chapter 3, and the preferred frame gravity models of §3.2.3, including 
those with modified dynamics at small accelerations, as in §3.2.2. 

Each of the modified acceleration laws applied in the astrophysics computations of 
Part II determine the acceleration felt by a test particle, at a distance r from the center. 
This acceleration is proportional to the mass enclosed within a spherical region of radius, 
r, so that 

aWocM(r), where M(r) = f p(r f )dV, (2.18) 

Jo 

where p(r) is the density at the position, r, and dV denotes the spherical volume element. 
In this case, the cells of data can be related by a factor, 

M N (r) = M(r)r(r) (2.19) 

where Mjv(r) is the dynamical mass of the integrated cells of data within a spherical region 
or separation, r, and M(r) is the observed baryonic mass of the same region. T(r) is the 
dynamical mass factor and is related to the dark matter ratio, whereby 

T(r) = 1 + M dark 

matter 

(r)/M(r), (2.20) 

where Md ar k matter (r) is the integrated mass of dark matter inside the common spherical 
region. At cosmological scales, where r is course grained away, the dark matter factor of 
Spergel et al (2007) is 

T = ^matter/^baryon = 5.73 ± 0.40 (2.21) 

The dark matter fits to the Ursa Major sample of §4.2 confirm that T ^ 10 across the 
galaxies, therefore the dark matter factor is consistent with the ACDM scenario. 

Alternatively, the dynamic mass factor predicted by MSTG and STVG, as in §3.3 and 
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§3.4, respectively, is effectively due to a renormalized gravitation coupling of Equation 
(3.90) with 

T(r) = M N {v)/M{v) = G{r)/G N , (2.22) 

where G(r) is the best-fitted gravitational coupling to the dynamical data at the separation 
r, and Gn = 6.67428(67) x 10 -11 m 3 kg -1 s -2 is Newton's constant measured experimen- 
tally 1 . 

Dynamic mass factors, constructed from galaxy rotation curves in the Ursa Major 
filament of galaxies, are provided in Figure 4.4, and those constructed from a sample of 
clusters of galaxies are provided in Figure 5.2. Conclusions drawn from the astrophysics 
on CDM halos may be found in §7.1.1. 



iNIST 2006 CODATA value. 



Chapter 3 



Modified gravitation theory 

"A hundred times every day I remind myself that my inner and outer life depend 
on the labors of other men, living and dead, and that I must exert myself in 
order to give in the same measure as I have received and am still receiving. " 

Albert Einstein 



In the absence of a fifth force in nature, either the dark matter paradigm ensues at 
astrophysical and cosmological scales, or the relativity principle may come into question. 
Local SO (3,1) invariance is a foundation of relativistic gravity theories, and is made mani- 
fest by general covariance. However, the equivalence principle, as in §3.1, may be violated 
by fifth- force fields or preferred space-time frames. Milgrom's modified dynamics (MOND), 
as in §3.2, is phenomenologically derived from observations of galaxy rotation curves and 
the Tully-Fisher relation, relativistic theories with a preferred frame, as in §3.2.3, are man- 
ifestly covariant, but violate SO (3,1) Lorentz covariance by means of a constraint. Moffat's 
metric skew-tensor gravity, as in §3.3, is a relativistic metric gravity theory, with massive 
fifth-force fields, Moffat's scalar-tensor-vector gravity, as in §3.4, is furthermore without 
phenomenological input from the Tully-Fisher relation. 

3.1 Equivalence principle 
3.1.1 Local SO(3,l) theory 

At the turn of the last century, Lorentz conjectured that Newton's universal gravitation 
law needed to be modified so that changes in the gravitation field propagate with the 
speed of light. Days before the Einstein (1905) paper on special relativity, Poincare (1904, 
1905) suggested that all forces, including gravity, should transform according to Lorentz 
transformations. Einstein set himself the task of modifying Newton's gravity theory to 
accommodate the principles of special relativity, and proposed the equivalence principle 
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based on the empirically observed universality of free fall: 
Einstein equivalence principle 

As far as we know, the physical laws with respect to an accelerated system 
do not differ from those with respect to a system at rest; this is based on 
the fact that all bodies are equally accelerated in the gravitational field. At 
our present state of experience we have thus no reason to assume that the 
accelerating and inertial systems differ from each other in any respect, and in 
the discussion that follows, we shall therefore assume the complete physical 
equivalence of a gravitational field and a corresponding acceleration of the 
reference system (Einstein, 1907). 

Einstein (1916) formulated his gravity theory geometrically so that particles travel along 
geodesies in a curved space-time. Observables are invariant under local Lorentz trans- 
formations, generalizing the global Lorentz invariance of special relativity. The Newto- 
nian gravitational attraction is the effect outside of a test particle's rest frame, modulo 
relativistic corrections. It is the curvature of the pseudo-Riemannian manifold which is 
fundamental, and the space-time metric is a dynamical solution to the Einstein equations. 

3.1.2 Strong equivalence principle 

The demand that the laws of nature, in a sufficiently small region of a given space-time 
point, take the same form as they do in special relativity is stronger than the universality 
of free fall as it means that there are no fields unified with the metric. 

The is contrary to the case in which the unified field is associated with Maxwell's elec- 
tromagnetism - which obviously does not vanish locally. However, the strong equivalence 
principle holds for pure gravity, where the unified field is constrained to vanish by the 
metric-connection compatibility equations. Conversely, the dynamic nature of the connec- 
tion field does not imply a vanishing torsion trace, and there are two degrees of violation 
of the strong equivalence principle in the general hermitian theory. 

Even though Einstein's theory may be written formally as a gauge theory, with field 
variables suitably chosen, it does not predict the form of the Newtonian universal force 
law. Bianchi et al. (2006) considered the graviton propagator within background indepen- 
dent, nonperturbative quantum gravity, yielding results that are consistent with Newton's 
universal law, but the renormalized interaction remains to be calculated. 
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3.1.3 Violations of the strong equivalence principle 

The two mechanisms for potential violations of the equivalence principle are: 
Charge violations The possibility that fermions possess quantum numbers related to 
conserved fifth force charges leads to direct violations of the weak equivalence prin- 
ciple and severely constrain modified gravity theory (Will, 2006). The non-abelian 
gauge theory for gravity necessitates the consideration that the nonmetric degrees 
of freedom that are associated with the larger symmetry group carry their charged 
quantum numbers. This would cause violations of universality of free fall, although 
Cavendish and Eotvos lunar laser ranging experiments tightly constraint any charge 
associated with the hermitian theory. The possibility of weak equivalence principle 
violations due to the Earth's rotation have been tightly constrained (Moffat and 
Brownstein, 1990), and are not expected to be measurable by the Gravity Probe B 
in Earth's orbit. 

Scalar-vector-tensor violations Unlike the local SO(3,l) theory in §3.1.1, local grav- 
ity in the general theory cannot be removed due to the presence of dynamical 
scalar/vector fields which are not determined by the metric. This is also a general 
property of scalar- vector-tensor modifications, including Brans-Dicke gravity theory. 
These strong equivalence principle violations do not of themselves imply any vio- 
lation of the weak equivalence principle: Scalar- vector-tensor gravity preserves the 
universality of free fall. This is important for the consideration of astrophysical 
phenomena, for which the universality of free fall is assumed. 
Each of the modified gravity theories, including Modified Newtonian dynamics, in §3.2, 
Metric skew-tensor gravity, in §3.3, and Scalar tensor vector gravity, in §3.4, violate the 
strong equivalence principle, but maintain the universality of free fall for bodies in motion. 

3.2 Modified Newtonian dynamics 

Milgrom's modified Newtonian dynamics (MOND) is a nonrelativistic small acceleration 
model which softens the Newtonian 1/r 2 force law to the 1/r behaviour preferred by galaxy 
rotation curves, introduced by Milgrom ( 1983a, b). MOND violates the strong equivalence 
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principle since, at sufficiently low accelerations, the gravitational mass of a test particle 
exceeds the inertial mass. 

Bekenstein and Milgrom (1984, Appendix B) showed that a modified Newtonian poten- 
tial may emerge, in the case of spherical symmetry, from a covariant Lagrangian formalism 
in which a cosmological scalar field, sourced by ordinary baryons, is added to the Einstein- 
Hilbert action. This relativistic, metric-scalar gravity theory is a modification of Jordan 
(1959); Brans and Dicke (1961) theory and similarly leads to violations of the strong equiv- 
alence principle, as in §3.1.3. 



c 4 



16ttG 



(j)R-2K-uj- 



(3.1) 



where the Jordan-Brans-Dicke parameter, cj, is not treated as a universal constant, but 
instead is treated as a function of the magnitude of the scalar gradient (Sanders, 1986a): 

»-U + I) M, (3.2) 



where 



c 4 4>„<t» v 



X 4(2^0 + 4)^3' ( 3 - 3 ) 
and a is the Milgrom universal acceleration parameter, and 

»(x) = (3.4) 

is the MOND interpolating function.. Sanders (1986a) extended the Bekenstein-Milgrom 
modification to include a fixed Yukawa-type length scale (Yukawa, 1935), which fits the 
galaxy rotation curves studied by Sanders (1986b) so that at cosmic distances from the 
source, the gravitationally strong MOND force would vanish entirely. 

To address the hypothesis stated in §1.1.3, of fitting galaxy rotation curves and galaxy 
cluster masses without dominant dark matter, Milgrom's acceleration law is presented in 
§3.2.1, and the resulting modified dynamics are considered in §3.2.2. In addition, the 
notion of building a relativistic, metric-scalar version of MOND is presented, and theories 
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with dynamical preferred frames including Bekenstein's TEVES theory and the generalized 
Einstein-aether theory are presented in §3.2.3. 

3.2.1 Milgrom's acceleration law 

Milgrom (1983a) challenged the hidden mass hypothesis and introduced a nonrelativistic 
modification of Newtonian dynamics (MOND) at small accelerations, a < a , whereupon 
the gravitational acceleration of a test particle is given by 



ajji\ 



a 

a 0i 



<2n, 



(3.5) 



where ji{x) is a function that interpolates between the Newtonian regime, ji{x) = 1, when 
x»l and the MOND regime, ji{x) = x ) when x <C 1. Milgrom (1983b) introduced the 
interpolating function normally used for galaxy fitting, 



n(x) 



x 



(3.6) 



where 



x = x{r) 



V$(r) 




a(r) 


a 




a 



(3.7) 



and determined that the MOND acceleration was of the order a ~ cH /6, and proportional 
to the Hubble constant, implying a cosmological connection to the modified dynamics. 
Substituting Equations (3.6) and (3.7) into Equation (3.5) gives, 



a(r) 2 



a(r) 2 



al 



a N ( r )> 



which has the solution, 



(3. 



a(r) = a 



\ 



'a N {rY 



a 



+ 



1 



'a N (r)' 



a 



+ 



'a N (r)' 



a 



(3.9) 
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written in terms of the Newtonian acceleration of a test particle at a separation, r, 

G N M{r) 

a N {r) = — , (3.10) 

where M(r) is the baryonic mass integrated within a sphere of radius, r. 

Milgrom's acceleration law, given by Equation (3.9), is applied to galaxy rotation curves 
in Chapter 4, in Equation (4.6). The galaxy rotation curves, plotted in Figure 4.1, 
are one parameter best-fits by the stellar mass-to-light ratio, T, applying the MOND 
acceleration, ao of Equation (4.8), universally. Milgrom's acceleration law is applied to 
clusters of galaxies in Chapter 5, according to §5.2.2. In §5.2, the MOND mass is best- 
fitted to the X-ray gas mass of a sample of 11 clusters of galaxies, and plotted in Figure 5.1 
according to the best-fit cluster model parameters tabulated in Panel (b) of Table 5.2, for 
Milgrom's MOND. 

3.2.2 Modified dynamics at small acceleration 

Substituting Equation (3.10) into Equation (3.5), the MOND acceleration law can be 
written, 

= 1 G^MM 

and therefore MOND can be interpreted as gravity theory with a varying gravitational 
coupling 

a { r) = (3,2) 
C(r) - £ y (3,3) 

and G(r) ~ Gn in the Newtonian regime and G{r) > Gn in the MOND regime. It 
is important to note that MOND has a classical instability in the deep MOND regime 
corresponding to \i — > which leads to a divergent gravitational coupling of Equation 
(3.13), and that MOND violates the strong equivalence principle for all \i ^ 1. 

For gravity fields interior to galaxies and clusters of galaxies, the accelerations are 
sufficiently small that the MOND interpolating function, ji{x) 3> 1, so that the Newtonian 
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dynamic mass determined by MOND is much larger than the actual mass visible in the 
system. 

Angus et al (2007) clarified the central issue in regards to gravitational lensing and 
the modified dynamics at small acceleration, since the total mass of baryons enclosed in a 
sphere of radius, r, is determined from the divergence theorem, 

1 ; J 47iG(r) dr ' 1 } 

where $ is the modified gravitational potential, and G(r) is given by Equation (3.13). 
Therefore, the MOND dynamic mass factor is precisely the inverse of the MOND interpo- 
lating function, 

r(r) = ^ = lMx{r)), (3.15) 

plotted in Figure 4.4 for the Ursa Major filament of galaxies, and in Figure 5.2 for the 
sample of X-ray clusters of galaxies. 

Bekenstein and Magueijo (2006) considered the behaviour of the MOND interpolating 
function in the deep MOND regime signalled by the small gradient of the dynamical scalar 
field, (j) of §3.2.3, where ji(x) « x and Equation (3.7) implies 



/i(r) 



V$(r) 




a(r) 


a 




a 



(3.16) 



In this regime, say far outside a spherically symmetric point source of mass, M, the Milgrom 
(1983a) acceleration law, 

f a\ G N M 
a/i — = — 5—, (3.17) 
\a J r 2 

simplifies upon substitution of Equation (3.16): 



a G N M ^/a G N M 
< r ) = ^—^ = " » ( 3 - 18 ) 

and thus the modified dynamics, at small acceleration scales, yields the gravitational field 
as 1/r instead of the Newtonian 1/r 2 law. Substitution of Equation (3.18) into Equation 
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(3.16) gives 



Mr)~r-\/^, (3.19) 

V «0 



which is valid in the deep MOND regime. Substitution of Equation (3.19) into the dynamic 
mass factor of Equation (3.15) gives 



r <'')='VG^ <3 ' 20) 

in the deep MOND regime, which shows a linear dependence with r at large distances. 
Milgrom and Sanders (2008) defined a transition radius in MOND, 



G N M 

n = \ , 3.21) 

V «0 

so that the dynamic mass factor of Equation (3.20) in the deep MOND regime, can be 
written, 

T(r) = -, (3.22) 

n 

The dynamical mass factors plotted in Figure 4.4, in Chapter 4, do indeed show 
a monotonically near-linear increasing T(r) < 10, reaching the maximum value at the 
outermost observed data point, r out , where r ~ 10 • r t , typically. This may imply that the 
MOND interpolating function is bounded from below, fi > 0.1. Otherwise, as the gradient 
of the scalar field approaches zero, and the MOND interpolating function of Equation 
(3.16) approaches zero, the dynamic mass factor of Equation (3.15) grows without bound, 
T(r) — ► oc indicating a classical instability. 

All of the modified gravity models examined in this dissertation provide the needed 
phantom dark matter, which is quantified by the dynamic mass factor, T > 1. For MOND 
this corresponds to /i(a/a^) < 1, although it is not known if the MOND interpolating func- 
tion approaches 0, this would correspond to the ultra-deep MOND regime and if MOND's 
dynamic mass factor is not bounded, r —> 0, would effectively renormalize gravity's cou- 
pling G —> oc. Conversely, if the inverse of the MOND interpolating function approaches a 
finite value, so fi^ < fi < 1, then the instability of the theory is made regular (finite), and 
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instead gravity's coupling approaches an asymptotic value, G —> G^. This is consistent 
with the Ursa major sample of §4.2, from which it is clear that ji^ < 10. Such a cutoff 
applied to clusters of galaxies could potentially cure MOND's unfortunate prediction of 
requiring dominant dark matter to fit clusters of galaxies, as in Chapter 5. However, the 
final form of MOND's interpolating function should be dynamically determined from - or 
at least correlated with - the action of the covariant field theory from which it is derived. 

3.2.3 Dynamical preferred frames 

The antithesis of Einstein's theory of special relativity, with local Lorentz SO (3,1) invari- 
ance, is the aether theory in which the symmetry is broken. Named after the luminiferous 
aether - the medium for the propagation of light as it was thought until the late 19th 
century - the aether theory is a generally covariant extension of general relativity by the 
addition of a unit timelike vector field. The aether has a preferred rest frame, and thus 
breaks local Lorentz SO(3,l) invariance. In an address delivered on May 5, 1920, at the 
University of Leyden, Einstein commented, 

How does it come about that alongside of the idea of ponderable matter, which 
is derived by abstraction from everyday life, the physicists set the idea of the 
existence of another kind of matter, the aether? The explanation is probably 
to be sought in those phenomena which have given rise to the theory of action 
at a distance, and in the properties of light which have led to the undulatory 
theory. 

It is the space-time components of the Maxwell field, F Mzy , which undulate; whereas the 
aether vector field is not free to undulate because it is constrained to spacelike oscillations, 
and the vacuum cannot be empty of aether excitations. 

Bekenstein's TEVES theory 

Bekenstein (2004) introduced the tensor-vector-scalar (TEVES) theory as a relativistic 
implementation of Milgrom's modified Newtonian dynamics (MOND), as in §3.2, with an 
additional scalar field, 0, and also a non-dynamical scalar field, a. The vector field in 
TEVES, A^, has timelike unit norm, 



g^A,A v = -1, 



(3.23) 
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and dynamically selects a preferred reference frame, breaking local Lorentz SO(3,l) invari- 
ance. g^ u is the Einstein metric, with a well defined inverse, g Xfl , such that 

g x »g»v = (3.24) 

However, all types of matter see the same physical metric 

fiU = e~ 2( V " 2sinh(20)A^, (3.25) 

with a well defined inverse, 

g Xfi = e'^g^ + 2smh(2(f))g Xa g^A a A (3 , (3.26) 

so adding a preferred frame is not in conflict with the weak equivalence principle. However, 
because TEVES is a relativistic, bimetric theory, it permits the computation of geodesies in 
the presence of matter sources, and makes predictions for lensing convergences, time-delays 
and other metric effects (Zhao, 2006). 

In TEVES, the vector field action is taken to be that of a Maxwell vector field, A^ Ji) 
with an additional Lagrange multiplier, A, to enforce the timelike, unit norm constraint of 
the vector field of Equation (3.23). The action for the pair of scalar fields, (/) and a, is a 
generalization of the Bekenstein (1988) phase coupling gravity (PCG) theory including a 
vector-scalar interaction. The total action for TEVES is formed by combining the Einstein- 
Hilbert action of Equation (3.98) with the vector and scalar actions: 



S\g,A,(f>,a] = J d^xyj^g 



c 4 



WttG 



R-2A- * F^F*" + A (A^A" + 1) 



a 2 (g^ - A^A V ) + -G£~ 2 a 4 T(kGa 2 ) 

Zj 



(3.27) 



where F^ = A v ^ — is the Maxwell vector field strength, K and k are dimensionless 
couplings, £ is a positive constant with units of length, and T is a free dimensionless 
function, similar to the PCG potential, whose behaviour is determined phenomenologically 
by requiring that the dynamics at slow accelerations correspond to MOND. 
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Bekenstein and Magueijo (2006) predicted a universal acceleration scale, in terms of 
the positive coupling constant, k ) and the length scale, I, of the TEVES action of Equation 
(3.27) 



a = « 10" 8 cm s" 2 , (3.28) 

4:7T£ 

consistent with the MOND acceleration of Equation (4.8). 

Zhao and Famaey (2006) argued that the Bekenstein (2004) model produces a MOND 
interpolating function with the wrong behaviour to accurately fit galaxy rotation curves; 
and suggested a refinement to the TEVES Lagrangian to accommodate the standard 
MOND interpolating function of Equation (3.6). 

Einstein-aether theory 

Jacobson and Mattingly (2001) proposed a generally covariant model in which local Lorentz 
invariance is broken by a dynamical unit timelike vector field, which is nowhere van- 
ishing. The Einstein-aether theory leads to gravity with a dynamical preferred frame, via 
the Jacobson and Mattingly (2004) action 



S[g, A] = S[g] - J d'x^g ^ (K a ^V a A»V pA v + X(A a A a - 1) 

where S[g] is the Einstein-Hilbert action of Equation (3.37), and 

K a \ u = c l9 af3 g^ + c 2 5 a ^ u + c&tf W + c 4 g^A a A^ (3.30) 

is written in terms of four dimensionless coefficients, c i) and A is a lagrange multiplier 
which enforces the unit timelike nature of the vector field. 

Jacobson (2007) reviewed the theory, phenomenology, and observational constraints on 
the coupling parameters of Einstein-aether gravity, showing that the unit timelike vector 
field, which breaks the local Lorentz invariance, must be dynamical; and the preferred 
frame must therefore be dynamical. 

Jacobson and Mattingly (2001) showed that such a field carries a nonlinear representa- 
tion of the local Lorentz group since the field does not take values in a vector space on the 



(3.29) 
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tangent space, but on the unit hyperboloid. Jacobson and Mattingly (2004) developed the 
linearized Einstein-aether theory, finding the speeds and polarizations of the wave modes, 
determining in addition to the usual two transverse traceless metric modes, three coupled 
aether-metric modes. Eling and Jacobson (2006) claimed that regular perfect fluid star so- 
lutions exist with static aether exteriors, with the aether field pointing in the direction of a 
timelike Killing vector, but there are no spherically-symmetric solutions constructed purely 
from the aether without naked singularities. Seifert (2007) applied the action of Equation 

(3.29) , and found that the flat space solution and the static vacuum aether solution of Eling 
and Jacobson (2006) is stable to linear perturbations, provided the coefficients q satisfy 
an auxiliary inequality relation. 

Clayton (2001) showed that Einstein-aether theories of the type of Equation (3.29) 
are energetically unstable, having a Hamiltonian, in Minkowski flat space-time, that is 
unbounded from below; and the linearized analyses about configurations with a vanishing 
aether vector field are singular. Jacobson (2007) pointed out that Clayton (2001) considered 
the question of energy positivity, but examined a limited Maxwellian subclass of Equation 

(3.30) in which c 3 = —C\ and c 2 = C4 = 0, and restricted to the case where the coupling 
to gravity is neglected. Seifert (2007) confirmed that the subclass investigated by Clayton 
(2001) has spherically symmetric static solutions which are unstable, likely related to the 
unbounded Hamiltonian, even though kinetic terms in the unit timelike vector for a range 
of coefficients, c i) that were ignored by Clayton (2001), may stabilize the theory. Jacobson 
(2007) suggested that the linear perturbations all have positive energy for coefficients, q, 
within a particular range, but the total nonlinear energy has not been shown to be positive 
in this range. It is an unsatisfactory situation that the theory requires special values of 
the q. This places too great a burden on phenomenology, limiting the theory's ability 
to make testable and falsifiable predictions, but there have been no successful attempts 
to identify a principle of symmetry to restrict the action. Carroll et al. (2009a, b) found 
that a timelike vector field leads to an unbounded Hamiltonian, and generates instability, 
except provided the kinetic term in the action is in the form of a a-model, and introduced 
a a-model aether modified gravity theory, with a timelike vector field. 

Arkani-Hamed et al. (2005) studied the effects of direct couplings between the Goldstone 
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boson (which appear due to the broken time diffeomorphism symmetry), and standard 
model fermions, which necessarily accompany Lorentz- violating terms in the theory, finding 
that the aether field couples to spin in the non-relativistic limit. A spin moving relative to 
the aether rest frame will emit Goldstone-Cerenkov radiation. The Goldstone boson also 
induces a long-range inverse-square law force between spin sources. 

Generalized Einstein-aether theory 

Zlosnik et al (2006) interpreted TEVES as a special case of the Einstein-aether theory of 
§3.2.3 with non-canonical kinetic terms, and showed that there exists a tensor-vector-scalar 
theory equivalent to TEVES, without the additional scalar field, 0, but retains the non- 
dynamical scalar field, a. The equivalent theory is cast as a single-metric theory, because 
the Einstein metric which satisfies the Einstein-Hilbert action couples minimally to the 
matter fields. However, there would be modifications to gravity resulting from the metric 
coupling to the vector field as a direct consequence of the Lorentz violating, dynamical 
aether. Zlosnik et al (2007) generalized the Einstein-aether theory of Eling and Jacobson 
(2004), replacing 

K = M- 2 K a ^ v V a A»VpA\ (3.31) 

by ^(/C), where is given by Equation (3.30), but restricted to a class of theories 

spanned by the first three coefficients q, i = 1,2,3, and M has the dimension of mass 
in order to make Equation (3.31) dimensionless. Although the generalized Einstein-aether 
theory of Zlosnik et al (2007) does not include the TEVES equivalent theory of Zlosnik 
et al (2006), each of these theories are reducible to MOND in the weak-field limit due to 
the never vanishing vector field. Carroll et al (2009a,b) found that because the Lorentz 
violating timelike vector field has kinetic terms in the action that are not in the form of 
a a-model, the theory leads to an unbounded Hamiltonian, and is not stable, whereas 
a-model aether modified gravity probably does not have a low acceleration MOND limit. 

Seifert (2007) considered the stability of spherically symmetric solutions in TEVES, 
without matter fields, finding that the perturbational Hamiltonian arising from the varia- 
tional principle has an indefinite kinetic term. In the absence of a well-defined variational 
principle with a sensible inner product, Seifert (2007) applied a WKB analysis to measure 
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the instability of the spherically symmetric vacuum solution, and predicted a timescale of 
10 6 seconds - two weeks - before a solar mass object would collapse under the weight of 
the nonvanishing vector-scalar fields. 

Contaldi et al (2008) confirmed that TEVES is a fully causal theory for positive values 
of the scalar field, and represents a relativistic modification of gravity which may depend 
on acceleration (since one must have a reference frame to measure the acceleration), but 
develops classical singularities which may prevent the weak acceleration limit from resem- 
bling MOND; and argued that caustic singularities are symptomatic of Einstein- aether 
theory, in general. However, Contaldi et al (2008) speculated that problems with the 
vector field dynamics may be rectified by choosing more general kinetic terms, which may 
also include MOND in the nonrelativistic limit. 

3.3 Metric skew-tensor gravity 

van Nieuwenhuizen (1973) found that the only massive antisymmetric tensor fields free of 
ghosts, tachyons and higher-order poles in the propagator for linearized gravitation are the 
massive spin-1 Maxwell-Proca fields. Isenberg and Nester (1977) performed a Hamilton- 
Dirac analysis of vector fields, determining that only Maxwell fields, Proca-Maxwell fields, 
and purely longitudinal vector fields are free of instability when minimally coupled to 
gravity. 

In light of the difficulty in obtaining physically consistent modified gravity theories, it 
is instructive to study the emergent Kalb-Ramond-Proca field, as in §3.3.1. The action 
in §3.3.2 for the metric skew-tensor gravity (MSTG) theory, given by Equation (3.36), 
couples an Einstein metric background (the metric sector) to the Kalb-Ramond-Proca 
field (massive skew sector). To address the hypothesis stated in §1.1.3, of fitting galaxy 
rotation curves and galaxy cluster masses without dominant dark matter, it is sufficient 
to work in the weak-field spherically symmetric limit of MSTG, where the test particle 
equations of motion, calculated in §3.3.3, are used to derive the point source acceleration 
law in §3.3.4 and effective Poisson equations are deduced in §3.3.5 for distributions of 
matter. The quadratic equations for the MSTG dynamic mass are solved exactly in §3.3.6 
by Equations (3.95) and (3.96). 
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3.3.1 Kalb-Ramond-Proca field 

Clayton (1996) showed that the massive nonsymmetric gravity theory (NGT) becomes 
identical to a Kalb-Ramond-Proca field with an additional curvature coupling term when 
considered as a perturbation about a Ricci-flat background. Since the Kalb-Ramond-Proca 
theory does not require a conserved current and yet has no negative energy ghost modes, 
higher order poles or tachyons, the additional terms in the action for massive NGT allow 
the linearized field equations to take on this form in the antisymmetric sector. 

Moffat (1995b); Moffat and Sokolov (1996) determined that in the weak-field approxi- 
mation relevant to galaxy dynamics, a range dependent Yukawa-type, fifth force (Yukawa, 
1935) emerges in addition to the Newtonian 1/r 2 central force due to the exchange of the 
spin-l + skewons between fermions; and asserted that this additional potential due to the 
interaction of the field structure with matter in the halos of galaxies can explain galaxy 
rotation curves, as in Chapter 4, and is a candidate for phantom dark matter. This 
hypothesis is studied in the dynamics of the weak-field, as in §4.1, using the modified 
gravity theory of §4.1.3, and extended to clusters of galaxies, in Chapter 5, with running 
gravitational couplings, as in §5.2.3. 

Geodesic and path motion in the nonsymmetric gravitational theory (NGT) were shown 
in Moffat (1995a) to have similar weak-field limits. The correction to the weak-field grav- 
itational force was found to be due to a Yukawa potential, resulting in a renormalized 
gravitational coupling. The Yukawa interaction, considered as an alternative to dominant 
dark matter, must account for the majority of astrophysical forces and meanwhile be com- 
pletely undetected at terrestrial scales. It is remarkable that the astrophysical studies in 
Part II show that the dark matter to baryon ratio can be consistently explained using the 
same Yukawa meson theory, from the smallest dwarf galaxies to the clusters of galaxies. 
Measurements in the weak-field, according to §4.1, provide support of the hypothesis that 
dark matter is a phantom of MSTG, with galaxy specific density distributions. Whereas 
the best-fitting dark matter theory, according to Equation (1.1) of §1.1.1, requires at least 
two additional dark matter parameters, po,^s, per galaxy, MSTG provides low reduced- 
X 2 best-fits with universal mean parameters across galaxy scales. Clusters of galaxies, 
however, show significantly improved x 2 best-fits with variable parameters. 
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Violations of the strong equivalence principle, described in §3.1.3, are the means by 
which scalar-vector-tensor modifications to the action for gravity result in a fifth force 
which preserves the universality of free fall. The effect due to the scalar- vector-tensor 
fields on the motion of a test particle requires careful approximation, such as the weak- 
field limit of a static, spherically symmetric space-time. At astrophysical scales, we neglect 
any contribution to the fifth force due to baryons with charged quantum numbers. Al- 
ternatively, we seek gravitationally strong contributions to the fifth force from a Yukawa 
(range dependent) meson emerging from the spin-l + massive vector skewon of the Kalb- 
Ramond-Proca field, as in §3.3.2. 

The first measurable predictions for galaxy dynamics in the NGT appeared in Moffat 
(1994) and Moffat and Sokolov (1996), where the appearance of a Yukawa-like potential 
produced by a new spin 1 + boson interacting with fermions emerged. In Legare and Moffat 
(1996), the effects of three new interactions were identified, and possible modifications to 
the geodesic and path motion were calculated in the weak-field limit. It was recognized 
by these attempts to provide an alternate explanation to the dark matter paradigm that 
the static, weak-field, slow speed, spherically symmetric limit of NGT may provide an 
adequate solution to the missing mass problem through the nonvanishing skewon mass 
and the coupling to baryons. In the static, spherically symmetric limit, the skewon field 
strength tensor, 

F\jau\] = d^g u x] (3.32) 

has only one independent, non-zero component, ify^j. The modifications to the radial 
orbit equations of motion were explicitly computed, and the surviving Yukawa contribution 
- potentially attractive or repulsive - added a new phenomenology to the dynamics of 
astrophysical scale measurements. 

For the case of a repulsive Yukawa potential added to the attractive Newtonian poten- 
tial, Sanders (1984) provided a preliminary analysis of circular orbit velocities in which 
the combined potentials lead to a deviation from the inverse square-distance law and may 
produce rotation curves which are "nearly flat from 10 to 100 kpc". Sanders (1984) spec- 
ulated that "a very low mass vector boson carries an effective antigravity force which on 
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scales smaller than that of galaxies almost balances the normal attractive gravity force." 

In principle, for each astrophysical experiment, the Yukawa coupling constant and 
the mass of the vector boson (Yukawa range) provide two additional parameters which 
may be modelled through the mass-to-light ratio. Sanders (1986b) provided a best-fit 
to six galaxy rotation curves ranging in size from 5 to 40 kpc to determine whether the 
modification to gravity is associated with a fixed length scale. Using the overall best-fit 
Yukawa coupling and range parameters, Sanders (1986b) computed mass-to-light ratios 
between 1 and 3, which are considered reasonable, showing no systematic variation with 
the size of the galaxy. The observed infrared Tully-Fisher law is shown to be consistent 
with the predictions of the Yukawa modified gravity for large galaxies greater in size 
than 15 kpc, whereas the smaller galaxies under 10 kpc do not exhibit a maximum flat 
rotation velocity. Admittedly, the sample is too small to statistically determine whether 
the parameters are universal constants, although best-fitting universal constants for the 
finite length-scale Yukawa repulsive gravity does lead to agreement with the data without 
introducing mass discrepancies. 

The issue of whether the Yukawa meson coupling and range are universal is not certain 
in the weak-field limit of NGT, where the Yukawa potential is emergent. Moffat (1995c) 
derived the mismatch between the six degrees of freedom in the full nonlinear theory, and 
the three degrees of freedom that survive in the symmetry reduced, and linear, weak-field 
limit due to a Kalb and Ramond (1974) field, identified clearly as the skewon, to 
explain the effective, low energy coupling to the Yukawa meson. 

Whereas Moffat (2004b) developed the radial orbit equations of motion for the problem 
of galaxy rotation curves from the full NGT action, Moffat (2004a) derived the linear weak- 
field approximation, from which the Kalb and Ramond (1974) field emerges as the field 
strength of the massive skewon. The modified acceleration law corresponds to the low 
energy, low speed limit of NGT, effectively suppressing the high energy contributions of 
the full theory. Metric skew-tensor gravity (MSTG) is introduced in Moffat (2005), where 
the modified acceleration law results from coupling the massive skew symmetric F^ y \ field 
to Einstein's metric. At astrophysical scales, the emergent low energy Yukawa meson is 
the only feature of the full NGT left in MSTG to explain galaxy rotation curves. 
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3.3.2 Action 

Damour et al. (1993) analysed a class of physically consistent and ghost-free nonsymmetric 
gravity models with finite range massive spin-l + gauge boson described by a second rank 
skew symmetric tensor, A^, with an action in which the skewon's field strength tensor 
is coupled to a conserved fermion current vector with a dimensionless coupling constant. 
The similarity to Maxwell's field, but for a massive skewon (instead of a massless photon), 
is described by the massive Kalb and Ramond (1974) action, 

S F = f d 4 x^g(^F^ x F^ x - ^M^), (3.33) 

with 

F^ x = d^A„ x + d„A XfM + dxA^ (3.34) 

and \i is the mass of the Kalb-Ramond-Proca field, A^. The action is invariant under 
diffeomorphisms, and invariant under the U(l) local gauge transformation, 

SoApv = d^e u - d v e^ (3.35) 

only in the massless case. Therefore, the dependence of the action based on the Lagrangian 
of Equation (3.33) must be quadratic in the Proca field's strength tensor, F^ v \ of Equation 
(3.34). 

Damour et al. (1993) showed that, although the Kalb-Ramond-Proca field leads to 
minuscule - as yet unmeasured - deviations from Newtonian gravity at terrestrial scales 
consistent with stringent bounds on possible violations of the weak equivalence principle, 
the field may acquire gravitational strength at sufficiently large astrophysical scales because 
the coupling is unbounded as the range increases, and that the magnitude of the field is 
proportional to the coupling. 

Moffat (2005) introduced the metric skew-tensor gravity (MSTG) gravity theory by 
adding the Kalb-Ramond-Proca action, coupled to a conserved fermion current, to an 
Einstein-Hilbert action: 

S = Seh + Sp + Sfm + Sm-> (3.36) 
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where 



c 4 



J d 4 xv^(i?-2A), (3.37) 



is the Einstein-Hilbert action, and Sp is the Kalb-Ramond-Proca action of Equation (3.33). 
A possible action for the coupling between the Kalb-Ramond-Proca field and matter was 
suggested by Damour et al (1993) in regards to the NGT, and applied to MSTG by Moffat 
(2005), according to 

S FM = J d A xF Xllu r x ^ = -3 J d 4 xe a ^A aP d,J„, (3.38) 

where J M is a conserved vector current and J*^ uX = e^ uXa J a is the dual tensor current 
density. Sm is the ordinary matter action. 

Varying the action of Equation (3.36) with respect to the metric, 

1 S S 1 

(Tm/2u + Tpiw) , (3.39) 



gives the field equations, 



GW + = ^ (T Mlw + T FpLV ) , (3.40) 

where T M ^ and T FybV are the energy-momentum tensors for matter and the Kalb-Ramond- 
Proca field, Api,, respectively. The Bianchi identities satisfied by the Einstein tensor, 

Gyv = Rfjuy — ~j9\ivR (3-41) 

lead to the conservation laws 

V^(T M ^ + T F ^) = 0. (3.42) 
Varying the action of Equation (3.36) with respect to the field, A^, 

1 SS 1 



-g SAv -sf^g 



(3.43) 
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gives the Kalb-Ramond-Proca field equations, 



(3.44) 



where J„ v is the tensor density source for the A„ v field: 



7 c f)& jfl 



(3.45) 



3.3.3 Motion under the fifth force 

The equations of motion of a test particle are 

W + { a/3 } U ^ = 9 " af ™ U ^ (3 - 46) 

where r is the proper time along the path of the particle and u x = dx x /dr is the 4- velocity 
of the particle, and 

|^| = \g Xp ( W + 9pw ~ w) , (3.47) 

is the Christoffel connection, and f aiy is derived from the Euler-Lagrange equations for a 
test particle, of mass m, and fifth force charge Am, where A couples the skew field to the 
test particle and is assumed constant for the universality of free fall. Such a test particle 
has a point particle action (Legare and Moffat, 1994; Moffat, 2005), 



f i f e acrr i x 

S TP -m / dr^g af3 u a uP - Am / dT—=F av xg a i3U P . (3.48) 
j j \j g 



Variation of Equation (3.48) yields Equation (3.46) with 



fau = ^d [a — j =F U7] x9v\ l i • (3.49) 



For a spherically symmetric, static skew symmetric potential field there are two 
non- vanishing components, the "magnetic" field potential A 0r (r) = w(r) and the "electric" 
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potential field Ag^r) = f(r) s'm9. According to Moffat (2005), we may assume that there 
are no static magnetic poles, so only the electric field contribution /(r)sin0 is non-zero, 
Therefore, has only one non- vanishing component: 

Fet r = d r Aet = f sinO, (3.50) 

where the prime notation implies differentiation with respect to r, so f = df /dr. Substi- 
tuting Equation (3.50) into Equation (3.49), 

For a static spherically symmetric gravitational field the line element has the form 

ds 2 = ~/(r)dt 2 - a(r)dr 2 - r 2 (d9 2 + sin 2 ## 2 ), (3.52) 

and the equations of motion for a test particle are 

rfV a^fdr\ 2 _ r_ ( d0\ 2 _ f sm 2 9 \ /#\ 2 V /^\ 2 
dr 2 2a\dr) a\drj \ a J \dr J 2a\drJ 

= °' (3 ' 54) 

d 2 9 2 (dd\(dr\ . n Jd4>\ 2 
d 2 4> 2(d(f)\(dr\ n n (d(f)\(d6\ n , n . 

The motion of a test particle can be shown to lie in the plane, 9 = 7r/2, by an appropriate 
choice of axes. Integrating Equation (3.56) gives 



(3.57) 
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where J is the conserved orbital angular momentum. Integration of Equation (3.54) gives 

(3.58) 



dt _ 1 

dr 7 



'' ; + E 



where E > is the conserved orbital energy per unit mass (E = for the photon). 
Substituting Equation (3.58) into Equation (3.53), and using Equation (3.57), we obtain 



cPr | a' (dr\ 2 J 2 | j ( A 7 /' 



~*~ ) cry dr ( \/ayr^) ( \fcrjr i ) ^ ^ 



<ir 2 2a\dr J ar 3 2a r y 2 y^/ajr 4 
3.3.4 MSTG acceleration law 

Approximating the line element of Equation (3.52) by the Schwarzschild solution: 

, , r 2GM\~ 1 , , 1 2GM 
«(r) - 1 , 7(r) - 1 , (3.60) 



and making the approximations 2GM/r <C 1, A/'/r 2 « 1, //r 2 « 1 and the slow motion 
approximation dr /dt <C 1, Equation (3.59) becomes 

d 2 r J% GM d//'\ 

W-^ + — = X dr-{7 2 )> (3 - 61) 

where Jn is the Newtonian orbital angular momentum. 

Transforming the Kalb-Ramond-Proca field equations of Equation (3.44), to polar co- 
ordinates for the components = f(r) sinr9, 

/ _ 2GM\ „ _ 2 / _ + 8GM\ 

\ r J r \ r J \r z r J 

which has the solution (Moffat, 2005) 

f(r) = ^G 2 M 2 exp(-/xr)(l + /xr), (3.63) 



where s is a dimensionless constant. The skew field is therefore an excellent candidate 
for the phantom of dark matter due to the result of Equation (3.63), which leads to 
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gravitationally strong astrophysical effects, with a geometric originating fifth force, similar 
to the appearance of a fifth force charge, Q5, as an integration constant in Equation (3.132), 
described in §3.4, for STVG. 

Substituting Equation (3.63) into Equation (3.61) gives 

d 2 r Jlr GM aexp(-/xr), N /0 ^ 

dp ~ ^ = — ^~ + — t^C 1 + H, (3-64) 

where a is given by 

a = . (3.65) 

It is phenomenologically important for the modified acceleration law to be consistent 
with the observed Tully-Fisher relation, v 4 oc M (Tully and Fisher, 1977), so MSTG 
requires that the constant, 5, be of the form 



s = gM~ 3/2 , (3.66) 



so Equation (3.65) becomes: 



where 



a = JMM , (3.67) 



M = (If) , (3,8) 

is a parameter related to the strength of the coupling of the Kalb-Ramond-Proca field to 
matter, and 

r = (3.69) 
/i 

and the gravitational constant, G, in Equation (3.64) is taken to be: 



Goo = G*(l + y^j, (3.70) 

where Gn is Newton's gravitational constant. If the dependence of a on the source mass 
distribution is correctly modelled as per Equation (3.65), then the MSTG mass parameter, 
M , defined by Equation (3.68), will be a universal constant, /i denotes the effective mass 
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of the skewon, A^, with reciprocal identified as the MSTG length parameter, r , defined 
by Equation (3.69). 

Substituting Equations (3.67), (3.69) and (3.70) into Equation (3.64), and neglecting 
the Newtonian angular momentum, J/v, we obtain the MSTG acceleration law, 



. G N M M 



1 - exp(-r/r ) 1 + 



We can rewrite Equation (3.71) in the form 



(3.71) 



(3.72) 



where 



G(r) = G N \l + 



'Mo 
M 



1 - exp(-r/r ) 1 + 



(3.73) 



3.3.5 Poisson equations 

The experience of a test particle in the MSTG theory, moving in an extended matter 
distribution, is a combination of the force of gravity due to Einstein's metric gravity theory, 
and a coupled fifth force due to a Kalb-Ramond-Proca field. The weak-field, central 
potential for a static, spherically symmetric system can be split into two parts: 



$(r) = $jv(r) + $y(r), 



(3.74) 



where 



and 



$iv(r) 



$ y (r) 



Gjvcrexp(— jir) 



(3.75) 



(3.76) 



denote the Newtonian and Yukawa potentials, respectively, where M denotes the total 
constant mass of a point source. The point source gravitational coupling in Equation 
(3.75) is 

a ■ (3.77) 



Goo - G N (1 + — ), 
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where a, defined by Equation (3.65), is dependent on the source mass distribution through 
a power-law model, derived from the Tully-Fisher relation, leading to the phenomenolog- 
ical parametrizations of Equations (3.67) and (3.70). Since the Schwarzschild solution, 
according to Equation (3.60), was used in the derivation of Equation (3.77), Equations 
(3.67) and (3.70) may be generalized to static, spherically symmetric matter distributions, 
using the interior solution of the Schwarzschild metric, and we may set M to the active 
mass interior to a sphere of radius, r, 

M = J dVp(r'). (3.78) 

The Poisson equations for <$>n( y ) and $y(r) are given by 

V 2 $ N (r) = -Goop(r), (3.79) 

and 

(V 2 - ^ Y (r) = ^G N p(r), (3.80) 

respectively. For sufficiently weak fields, the Poisson Equations (3.79) and (3.80) are 
uncoupled and determine the potentials <&n( y ) and $y(r) for non-spherically symmetric 
systems, which can be solved analytically and numerically. The Green's function for the 
Yukawa Poisson equation is given by 

(V 2 -/i 2 )A y (r) = -5 3 (r). (3.81) 

The full solutions to the potentials are given by 

*"M = - G "M 1 + 5)SF^ < 3 ' 82 > 

and 

*y(*) = G N J d ^'^M-^-v'\)^-. (3.83) 
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The modified acceleration law is the gradient of the potential of Equation (3.74), 



(r) = -V$ = -(V$jv(r) + V$y(r)). 



(3.84) 



Combining Equations (3.82), (3.83) and (3.84), 



a r 



° N I ^ ^Ir-ff l 1 + li t 1 " eXP( "^ |r " r1)(1 + Hl * " r1) 



(3.85) 



Dividing the Tully-Fisher relation inspired phenomenological input of Equation (3.67) 
by M, 

(3.86) 



a_ = Mo 



and substituting Equations (3.69) and (3.86) into Equation (3.85), we obtain 



'« = -/ 



47r|r-r'| 3 ' 



(3.87) 



where 



G{t-t>) = G n {\ + 



M 



1/2 



Jd 3 r f p(r') 



1 — exp 



|r — r'l \ / ^ |r — r'l 



r 



1 + 



For a ^-function point source, 



(3.88) 



p(r) = M5 3 (r), 



(3.89) 



the modified acceleration law of Equations (3.87) and (3.88) reduces to the point source 
solution of Equations (3.72) and (3.73). 

For a static, spherically symmetric system, the effective modified acceleration law is: 



a{r) 



G{r)M(r) 



(3.90) 



where 



G(r) = G N {1 + 



M(r) 



1 - exp(-r/r ) ( 1 + 



r . 



(3.91) 
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We observe that G(r) — > G^ as r — > 0. 



3.3.6 Dynamical mass measurements 

Comparison of Equation (3.90) with the Newtonian acceleration law: 



a(r) 



G N M N (r) 



(3.92) 



allows the interpretation of the modified gravity dynamic mass as a scaled version of the 
Newtonian dynamic mass, 



Mmog(^) 



G N M N (r) 
G(r) 



(3.93) 



where the varying gravitation coupling, G(r), may take the form of Equation (3.91), derived 
in MSTG. The MSTG dynamic mass, 



M MSTG (r) = M N (r)\l + 



has the exact analytic solution: 



Mn 



M M stg(V) 



1 - exp(-r/r ) 1 + 



(3.94) 



M MST G(r) = M N (r) + M £(r) - ^M 2 f(r) 2 + 2M M N (r)^r), (3.95) 



£(r) = 



1 



1 - exp(-r/r ) 1 + 



(3.96) 



which is identified with the total baryonic mass within a separation, r from the center of 
the system. 

This MSTG acceleration law is applied to galaxy rotation curves in Chapter 4, in 
Equations (4.11) and (4.12). In §4.2 , in order to compute the overall best-fitting mean 
result, M and r are permitted to vary across the sample of 29 galaxies, as tabulated 
in Table 4.4. The galaxy rotation curves, plotted in Figure 4.1, are subsequently one 
parameter best-fits by the stellar mass-to-light ratio, T, applying the mean results of 
Equation (4.14) universally. The MSTG acceleration law is applied to clusters of galaxies 
in Chapter 5, according to §5.2.3: Metric skew-tensor gravity, in order to compute the 
scaling of the parameters, M and r . In §5.2, the MSTG mass is best-fitted to the X-ray 
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gas mass of a sample of 11 clusters of galaxies, and plotted in Figure 5.1 according to 
the best-fit cluster model parameters tabulated in Panel (c) of Table 5.2, for MSTG. A 
summary of lessons learned from the application of MSTG to the astrophysics of galaxies 
and clusters of galaxies is supplied in §7.1.2. 

3.4 Scalar-tensor-vector gravity 

Whereas the metric skew-tensor gravity theory, of §3.3, describes the effective, low energy 
Yukawa skewon as the Kalb-Ramond-Proca field, a separate solution is to model a gravity 
theory with a simpler Maxwell-Proca field, such as in the scalar-tensor-vector gravity 
(STVG) theory, which describes the low energy Yukawa phion, as a massive spin- 
1~ vector field, described in §3.4.1. The action in §3.4.2 for the STVG theory, given 
by Equation (3.97), includes an Einstein metric background (the gravity sector) to the 
Maxwell-Proca field in which the gravitational coupling, G, and the phion coupling, uj ) and 
the phion mass, /x, are treated as a triplet of scalar fields (scalar-tensor- vector sector). To 
address the hypothesis stated in §1.1.3, of fitting galaxy rotation curves and galaxy cluster 
masses without dominant dark matter, it is sufficient to work in the weak-field spherically 
symmetric limit of STVG, where the test particle equations of motion, calculated in §3.4.3, 
are used to derive the point source acceleration law in §3.4.4 and effective Poisson equations 
are deduced in §3.4.5 for distributions of matter. The STVG dynamic mass is provided 
in §3.4.6 by Equation (3.164) which is nonlinear through Equations (3.159) and (3.160), 
and requires a numerical solution unlike the exact analytic solution of Equations (3.95) 
and (3.96) for MSTG. Using the MSTG dynamic mass as the initial guess for the STVG 
numerical computation led to fast convergence in fewer than ten iterations. Lessons learned 
from the application of STVG to the astrophysics of galaxies and clusters of galaxies may 
be found in §7.1.2. 

3.4.1 Maxwell-Proca field 

van Nieuwenhuizen (1973) found that the only massive vector fields free of ghosts, tachyons 
and higher-order poles in the propagator for linearized gravitation are the massive spin-1 
Maxwell-Proca fields. 
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The expectation from the Yukawa skewon theory of §3.3.4 is that the gravitational 
coupling, G(r) of Equation (3.91) and the mass and range parameters, M and r , are 
scale dependent. Conversely, STVG theory models this astrophysical scale dependence 
with a renormalized triplet of self-interacting, cosmological, Klein-Gordon scalar fields. 

Moffat (2006b) introduced the scalar-tensor- vector gravity (STVG) theory by including 
a massive spin-1 - vector phion, which is the Maxwell-Proca field of §3.4.2, self-coupled and 
coupled to a matter current, to an Einstein-Hilbert action. Perhaps much simpler than the 
NGT, and possibly MSTG, the STVG effectively captures the fifth force due to a weak- 
field emergent Yukawa meson, simulating the predictions of NGT and MSTG, to a first 
order approximation. To address the hypothesis stated in §1.1.3, of fitting galaxy rotation 
curves and galaxy cluster masses without dominant dark matter, it is sufficient to work in 
the weak-field limit where the effective, low energy excitation is described by the Yukawa 
phion theory of §3.4.4. The cumulative renormalization of the phion mass, /i, self-coupling, 
a;, and the gravitational coupling, G, provide the gravitational strength. The central force 
law, for test particle motion in STVG, is derived in §3.4.5. 

3.4.2 Action 

The STVG action, with matter present, is based on the Lagrangian density, 



provides the general relativistic background, where A is the cosmological constant. The 
Maxwell-Proca spin-1 - vector phion, introduces the fifth force by a modification to 
gravity's action by the inclusion of the Lagrangian density, 



£ — £>eh + £0 + £>s + £m- 



(3.97) 



The Einstein-Hilbert Lagrangian density, 




(3.98) 




(3.99) 
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where fi is the phion mass, uj characterizes the coupling strength between the phion and 
matter, is the phion self-interaction potential, and the phion field strength tensor is 

B^ = d^-d^. (3.100) 



Isenberg and Nester (1977) showed that, when minimally coupled to gravity, both the 
Maxwell field photon, A a: and the Maxwell-Proca field phion, (j) a: where a = (0,a), have 
two constraints: The primary constraint sets the canonical momentum conjugate to A 
and 0o to zero, and therefore the longitudinal modes are non-propagating. A secondary 
constraint enforces the Gauss law, 

{0 photon 
(3.101) 
[i 2 (j)o phion, 

where E a is the canonical momentum conjugate to A a) or a , respectively. The Maxwell 
Hamiltonian has undetermined Lagrange multipliers which generate U(l) gauge invariance, 
but the Maxwell-Proca Hamiltonian is uniquely determined, since it is U(l) gauge non- 
invariant. Therefore, the Maxwell field has 4—1 — 1 = 2 degrees of freedom, whereas the 
Maxwell-Proca field has 4—1 = 3 degrees of freedom. Similar arguments apply to the 
MSTG massive spin-l + skewon, of §3.3.2, which is a Kalb-Ramond-Proca field. 

Moffat (2006b) confirmed that there are no pathological singularities in the Maxwell- 
Proca field coupled to gravity and promotes the three coupling constants of the theory, 
G, ji and a;, to scalar fields by introducing associated kinetic and potential terms in the 
Lagrangian density: 



£s = -g v I £2 + —^2 V„u;V„u; I + + + V u {u) 



G 



v-g, 

(3.102) 

where denotes covariant differentiation with respect to the local SO(3,l) invariant, 
symmetric metric g vu , while Vq, V^, and V w are the self- interaction potentials associated 
with the scalar fields. 
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The action principle for STVG in the presence of matter, 

SS = sj d 4 x (C EH + C^ + C S + C M ) = 0, (3.103) 

where Cm is the ordinary matter Lagrangian density, with Sm = / cPxCm- 
The total energy-momentum tensor takes the form, 

T^ u = Tm/j,v + T 4,^ + Tspv, (3.104) 

where 

2 $Sm __ rp 2 8S4 2 SS S _ , Qin ,^ 

Variation of the action with respect to g^ v yields the Einstein field equations in the 
presence of a massive vector phion: 

GW + + = (3.106) 
where is the Einstein tensor given by Equation (3.41), and 

- si { + ^ - v^) (« - I S «V) } 

+ S^j^ + ^ + w}. (3.107) 

A fifth force-matter current arises from extremizing the matter action under variations of 
the Maxwell-Proca phion field, 

j» = _1SSm_ (3108) 
Variation of the action with respect to <p u yields the Maxwell-Proca equations for the 
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massive vector phion: 



2 xv 



avM _ i 



r. 



(3.109) 



Variation of the action with respect to the gravitational coupling, G, the phion coupling, 
u>, and the phion mass, ji yields the scalar field equations: 



2 G ' 2 
-^(G) + G 



G_ 

4c 4 " 



G 

16tt 



2c 4 



+ 



c 21 



(R-2A) = 0, (3.110) 

0, (3.111) 
0. (3.112) 



3.4.3 Motion under the fifth force 



The equations of motion of a test particle are 



(3.113) 



where r is the proper time along the path of the particle and u x = dx x /dr is the 4- velocity 
of the particle, and 



9 p (9w,v + 



P) > 



(3.114) 



is the Christoffel connection. The acceleration, 0,5, is due to the fifth force derived from 
the Euler-Lagrange equations for a test particle, of mass m, and fifth force charge, 



<?5 «m, 



(3.115) 



where ft is a constant, independent of m. Such a test particle has a point particle ac- 
tion (Moffat, 2006b; Moffat and Toth, 2009), 



S TP = —m J dr^J g a pu a uP — q§ J druxj)^ . 



(3.116) 



§3.4. SCALAR-TENSOR- VECTOR GRAVITY 59 

Variation of Equation (3.116) yields the Euler-Lagrange equations corresponding to Equa- 
tion (3.113), where the velocity-dependent fifth force is given by 

ft = q, \uB\u v + ((/> a u a ) - V a uj (<pu a )] . (3.117) 

Dividing the fifth force of Equation (3.117) by the test particle mass, m, and using Equation 
(3.115), the mass m cancels, and the acceleration due to the fifth force becomes, 

a% = ^ = kuoB» v u v + kV^uj UaU a ) - nV a u , (3.118) 

m 

which is independent of the test particle mass, in exact agreement with the weak equiva- 
lence principle and the universality of free fall. Taking uj as constant, 

a% = kmB^vu". (3.119) 

For a static spherically symmetric gravitational field the line element has the form 

ds 2 = Bdt 2 - Adr 2 - r 2 (d9 2 + sin 2 0d<f> 2 ), (3.120) 

and the equations of motion for a test particle are 

rfV a^fdrV _ r_ ( d6\ 2 _ / sin 2 fl \ / Kwd<f> \ 2 dtV 
dr 2 2a\dr) a\drj \ a J \ dr J 2a\drJ 

d 2 6 2 (d6\(dr\ n „fd(f>\ 2 
d 2 (f) 2(d(j)\(dr\ n Jd<t>\(dB\ n , n , ni . 



The motion of a test particle can be shown to lie in the plane, = 7r/2, by an appropriate 
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choice of axes. Integrating Equation (3.124) gives 

r 2 ^ = J, (3.125) 
where J is the conserved orbital angular momentum. Integration of Equation (3.122) gives 



dt _ 1 

dr 7 



KW(/)q + E 



(3.126) 



where E > is the conserved orbital energy per unit mass {E — for the photon). 
Substituting Equation (3.126) into Equation (3.121), and using Equation (3.125), we obtain 



d 2 r a' ( dr\ 2 J 2 7' , _ 9 1 / dd \ , , _ , H . 

3.4.4 STVG acceleration law 

Approximating the line element of Equation (3.120) by the Schwarzschild solution: 

, , r 2GM\~ 1 . , 1 2C7M 

a(r)~(l , 7 r ~ 1 , 3.128) 

\ r J r 

and making the approximations 2GM/r <C 1, nuxfro <C 1, and the slow motion approxima- 
tion dr/dt <C 1, Equation (3.127) becomes 

d 2 r J% GM # ,„ 10Q , 

where Jjv is the Newtonian orbital angular momentum. 

In the limit of no phion self-interactions, V^(0) — > 0, and with constant, Equation 
(3.109) becomes 

V^ + ZiV = (3.130) 

In the weak-field, static spherically symmetric limit with J v = 0, the only nonpropagating 
mode, 0o, obeys the Maxwell-Proca equation 
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which has the Yukawa solution 

e -fir 

Mr) = -Qb • (3.132) 

r 

The constant Q$ emerges as a constant of integration, and should be interpreted as an 
effective Yukawa phion field strength, whereas the mass of the effective Yukawa phion, 
/x, should be interpreted as the range of the Yukawa interaction, A = l//i. Substituting 
Equation (3.132) into Equation (3.129), 

d 2 r J 2 N GM kujQ 5 exp(-/xr) 

dp ~ ^ = -—2- + ~2 (! + H- (3-133) 

Since the effective phion field strength is proportional to the source mass, M, with the 
same constant of proportionality as in Equation (3.115), (Moffat and Toth, 2009) 

Q 5 = kM, (3.134) 



we may write Equation (3.133) as 

d 2 r J% GM aG N Mexp(-/Lir) 
dt 2 r 3 r 2 r 2 

where 



+ (3.135) 



aG N = k 2 uj. (3.136) 

Demanding consistency with the observed Newtonian force law, for small r, when fir <C 1, 
the difference between Equation (3.135) and the Newtonian force law vanishes, 

G N M GM aMG N „ 

— 2 r + ^ = ' 3.137 

and the gravitational constant, G, in Equation (3.135) has the solution 

G OQ = G N (l + a). (3.138) 



Substituting Equation (3.138) into Equation (3.135) and neglecting the Newtonian angular 
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momentum, J N , we obtain the STVG acceleration law, 

exp(— fir) 



n G N M 
a(r) = —< l + a 



(1 +/ir) 



(3.139) 



Whereas the derivation of the MSTG acceleration law of Equation (3.71) relied upon the 
phenomenological input of Equations (3.66) and (3.67) leading to the MSTG parameters 
M and r , Moffat and Toth (2009) integrated the equations of motion in the weak-field, 
spherically symmetric limit, obtaining a and \i as functions of the mass M, 

a = ™ , 2 (§*-l). (3-140) 

D 

li = -=. 3.141 

y/M 



The parameters D and E are universal constants. We can rewrite Equation (3.139) in the 
form 

a(r) = (3.142) 

where 

G(r) = G N {l + a - ae~^(l + fir)} . (3.143) 
3.4.5 Poisson equations 

The experience of a test particle in the STVG theory, moving in an extended matter 
distribution, is a combination of the force of gravity due to Einstein's metric gravity theory, 
and a fifth force described by a triplet of scalar fields and a Maxwell-Proca field. The 
weak-field, central potential for a static, spherically symmetric system can be split into 
two parts: 

®(r) = $7v(r) + $y(r), (3.144) 

where 

Mr) = - — , (3-145) 
r 

and 

<Mr) = ^ 95eXP( "^ r) (3-146) 
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denote the Newtonian and Yukawa potentials, respectively, where M and Q$ denote the 
total constant mass and fifth force charge of a point source. is the gravitational 

coupling in Equation (3.145), and is given by Equation (3.138), and \± denotes the effective 
mass of the phion, in STVG. Since is proportional to M by Equation (3.134) and 
using Equation (3.136), Equation (3.146) becomes 

= cxG N M exp(-fir) ^ 
r 

The Poisson equations for $jv(r) and <&y(r) are given by 

V 2 <Mr) = -GooP(r), (3.148) 

and 

(V 2 - /x 2 )$ y (r) = <*G nP {t\ (3.149) 

respectively. For sufficiently weak fields, the Poisson Equations (3.148) and (3.149) are 
uncoupled and determine the potentials $w(r) and $y(r) for non-spherically symmetric 
systems, which can be solved analytically and numerically. The Green's function for the 
Yukawa Poisson equation is given by 

(V 2 -/i 2 )A y (r) = -8\v). (3.150) 

The full solutions to the potentials are given by 

* m = -G„f^^m ,3,51) 



4yr|r — r' 

and 

,ap(r') exp(— /i|r — r'|) 



$ y (r) = G N j d ^)^-^-r\) _ 
J 47r|r — r'| 

The modified acceleration law is the gradient of the potential of Equation (3.144), 



a(r) = -V$ = -(V$iv(r) + V$y(r)). 



(3.153) 
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Combining Equations (3.151), (3.152) and (3.153), 



a(r) = -G N I ^ r / ( r r )p( r ) ji + a - ae xp(-dr -r 7 |)(l + dr -r 7 |)l. (3.154) 
J 47T\r — r'| d [ J 



Therefore 

a M = - ( r1 3 v f - 

47r|r — r 



a ( r ) = " / ^Li/^W - A (3-155) 



where 

G(r - r') = (?7v|l + ol - aexp(-/x|r - r'|)(l + /i|r - r'|) j. (3.156) 
For a 5-function point source, 

p(r) = M5 3 (r), (3.157) 

the modified acceleration law of Equations (3.155) and (3.156) reduces to the point source 
solution of Equations (3.142) and (3.143). 

For a spherically symmetric system, the total baryonic mass within a separation, r, 
from the center of the system, is 

M(r) = f ^r' 2 dr'p(r'). (3.158) 
Jo 

Whereas the MSTG Poisson equations of Equations (3.79) and (3.80) relied upon the Tully- 
Fisher relation inspired phenomenological input of Equations (3.67) and (3.70) leading 
to the MSTG parameters M and r in Equation (3.91), a and /i can be obtained, by 
integrating the equations of motion in the weak-field, spherically symmetric limit (Moffat 
and Toth, 2009): 



(y/Mfrj + E)' 



D 

(i = , 3.160) 

lM{r) 
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For a static, spherically symmetric system, the effective modified acceleration law is: 

«M = -C^l, (3,61) 

where 

G(r) = G N {l + a - ae~^(l + fir)} . (3.162) 
where a and fi are given by Equations (3.159) and (3.160), respectively. 
3.4.6 Dynamical mass measurements 

Comparison of Equation (3.161) with the Newtonian acceleration law of Equation (3.92) 
allows the interpretation of the modified gravity dynamic mass as a scaled version of the 
Newtonian dynamic mass, 

M M oG{r) = — y^—, — , (3.163) 
G{r) 

where the varying gravitation coupling, G(r), may take the form of Equation (3.162), 
derived in STVG. The STVG dynamic mass, 

MsxvG(r) = M N {r) {l + a - ae~^(l + fir)}' 1 , (3.164) 

where a and fi are defined by Equations (3.159) and (3.160), respectively, is identified with 
the total baryonic mass within a separation, r from the center of the system. 

This STVG acceleration law is applied to galaxy rotation curves in Chapter 4, in 
Equations (4.11), (4.15), (4.16) and (4.17). In §4.2, in order to compute the overall best- 
fitting mean result, L>, E and are permitted to vary across the sample of 29 galaxies, as 
tabulated in Table 4.5. The galaxy rotation curves, plotted in Figure 4.1, are subsequently 
one parameter best-fits by the stellar mass-to-light ratio, T, applying the mean results of 
Equation (4.19) universally. The STVG acceleration law is applied to clusters of galaxies 
in Chapter 5, according to §5.2.3: Scalar-tensor- vector gravity, in order to compute the 
scaling of the asymptotic coupling, G^. In §5.2, the STVG mass of Equation (3.164) is 
fitted to the X-ray gas mass of a sample of 11 clusters, and plotted in Figure 5.1 according 
to the best-fit cluster model parameters tabulated in Panel (d) of Table 5.2, for STVG. 



Part II 
Astrophysics 



Chapter 4 



Galaxy rotation curves 

"To myself I am only a child playing on the beach, while vast oceans of truth 
lie undiscovered before me. " 

Sir Isaac Newton 



The creation of galaxy rotation curves from astrophysical observations is subject to 
model dependent assumptions. The road from photometry, in some observed electromag- 
netic band, to the mass profile, in some chosen gravity theory, takes its way through the 
dynamics of the galaxy, with the destination a rotational velocity profile. Spiral galaxies 
show a remarkable variation of the distribution and relative abundances of stellar material, 
distributed in bulges and disks, and the intergalactic medium, distributed in exponentially 
thin rings with vanishing amounts within galaxy cores. These are the three visible com- 
ponents - the sources of photometric data - that are used to reconstruct the dynamics of 
the galaxy, and the predicted galaxy rotation curves. 

The galaxy mass profiles are determined by a best-fit algorithm, within each gravity 
theory depending on dynamics in the weak-field, as in §4.1, for a sample from the Ursa 
Major filament of galaxies, in §4.2. Every galaxy studied, from the highest to lowest 
in surface brightness, from the most giant to the smallest dwarf, require some form of 
dark matter or some modification of gravity. Each of the candidates offer robust and 
distinguishable predictions for the mass luminosity relationship, as in §4.2.8. Although 
Milgrom's modified Newtonian dynamics and Moffat's modified gravity theories are sourced 
by ordinary baryons, there is evidence that each of these theories lead to measurable, 
and distinguishable halos of phantom dark matter, as described in §4.3. Dark matter 
distributions are sensitive to baryon distributions because x 2_ fitting algorithms recover 
the kinks and wiggles, repatriating the surface masses of orphan features, described in 
§4.3.2. If the best-fit stellar mass-to-light ratio for the model is near unity, T ~ 1, then 
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the Tully-Fisher relation, as in §4.3.4, follows from fundamental physics. 
4.1 Curve- fitting 

The observational data from galaxy rotation curves is compared to the predictions of cold 
non-baryonic dark matter (CDM) halos, Milgrom's modified Newtonian dynamics, and 
Moffat's modified gravity theories, in §4.1.1, §4.1.2, and §4.1.3, respectively. 

4.1.1 CDM halos 

In Chapter 2, the halo density power-law function of Equation (2.5), was shown to have 
fitting formulae for a power-law index between 1 < 7(7-) < 3, for the NFW formulae of 
Equations (2.7) and (2.8), 

p{r) = 



r(r + r s ) 2 ' (4.1) 
M(r) = 47rpo^ {ln(r + r 5 ) - ln(r 5 ) - 

and between < 7(7") < 3, for the core-modified formulae of Equations (2.12) and (2.13), 

pin - 



r + r t (4.2) 
M{r) = §7rp H? {ln(r 3 + rf) - ln(rf)} . 

Each profile is self-similar and describes the tremendous variation of galaxy-scale halos, 
without any further parameters. In addition, each profile has a simple analytic expression 
for the integrated mass function, M(r), relevant for curve-fitting. The two parameters 
that must be varied in both Equations (4.1) and (4.2), are po and r s . In the core-modified 
model these can be interpreted as the dark matter central density, and the radius at which 
the density is one-half the central density, respectively. Furthermore, the core-modified 
7^0 behaviour occurs in the baryon dominated galactic core, decreasing the dark matter 
density where the cusp problem prevents better fits using the NFW profile. 

Best-fits to the mass profiles of the dark matter halos, neglecting the stellar galactic 
disk (T = 0) were poor to gross for both the NFW profile and the core-modified model, 
whereas simultaneously best-fit parameters, po, T, produced low to very low x 1 ) as 
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shown in Table 4.2. The very low x 2 best-fits repatriated many of the orphan features seen 
within the galaxy rotation curves, as shown in Figure 4.1. Moreover, the predicted surface 
mass profile, S(r), extends gradually into the galaxy with a much broader center than the 
predictions of the modified gravity S-maps, as shown in Figure 4.2. 

Observations 

The CDM computations, using HI and K-band photometric surface luminosity data, de- 
tailed in §4.2, with galaxy rotation curves plotted in Figure 4.1, indicate the following: 

1. The sample may be universally fit with a common NFW profile given by Equation 
(4.1), where the NFW parameters are varied in order to best-fit the rotation curve - 
either with or without baryons. The fits without baryons lead to gross best-fits of the 
galaxy rotation curves, with very poor x 2 - Including the visible HI (and He) gaseous 
disk and the available luminous stellar disk with a variable stellar mass-to-light ratio 
T provides excellent fits to the large galaxies, but suppresses the best-fit T « 1, 
particularly in the case of the dwarf galaxies. The worst of these dwarf galaxies 
cannot be fitted using the NFW profile with any nonzero value of the stellar mass-to- 
light ratio. This confirms the cusp problem due to the singular NFW fitting formula 
and reinforces the importance of correctly incorporating the baryonic components 
into the galaxy models. 

2. Every galaxy in the UMa sample, from the highest to lowest in surface brightness 
may be universally fit with a common core-modified profile given by Equation (4.2) 
- with no extra parameters beyond those of the NFW parameters - provided the 
visible HI (and He) gaseous disk and the available luminous stellar disk are included. 
This model provides superior fits, with the lowest reduced x 2 statistic over all of 
the gravity theories considered, including all of the dwarf galaxies, and yields values 
of T ~ 1 as tabulated in Table 4.2. Moreover, the dark matter to baryon ratio 
at the outermost radial point, tabulated in Table 4.6 with mean values provided by 
Equation (4.53), is below the upper limit set by Spergel et al (2007) in the Wilkinson 
microwave anisotropy probe (WMAP) third year results. 

3. Every galaxy in the UMa sample, from the highest to lowest in surface brightness 
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has a central disk where the dark matter density differs strongly from a single power- 
law density profile, and 7(7-) of Equation (2.5) increases with radii, r, as shown in 
Figure 4.5. This solution to the dark matter cusp problem is studied in §4.2.7. 

4. The UMa sample can be fit by Newton's theory alone - without dark matter - using 
the visible HI (and He) gaseous disk and the available luminous stellar disk, within 
a Newtonian core up to some radius which varies across the galaxy sample. This 
maximizes the stellar mass-to-light ratio, T, and leads to systematically bad fits 
beyond the Newtonian core, indicating that the missing mass problem increases with 
radius. This point is elaborated in §4.3.1. The Newtonian core radii are plotted in 
Figure 4.6, and the galaxy rotation curves derived from this best Newtonian core 
model are plotted in Figure 4.1. 

5. The total mass and the shape of the dark matter halo varies in all galaxies, indepen- 
dent of the total mass of the visible HI (and He) gaseous and stellar disks. The dark 
matter parameters are neither correlated with galactic mass, nor the flat rotation 
velocity, v out , nor with the extent of the galaxy rotation curve, r out) as listed in Ta- 
ble 4.1. Sub kiloparsec, high-resolution E-map predictions are provided in Figure 4.2. 

6. Orphan features become traceable to a parent in either the gaseous disk, or the 
luminous stellar disk, for r < r s , but become increasingly orphaned for r ^> r s 
where the dark matter halo dominates. This provides the most obvious improvement 
between the quality of the fits, as compared to those of the NFW profile. 

Conclusions drawn upon identification of the missing mass as CDM is presented in the 
summary §7.1.1. 

4.1.2 Milgrom's acceleration law 

In Chapter 3, Milgrom's acceleration law of Equation (3.5), 



aji(x) <2n 



(4.3) 



with the interpolating function, 



x 



(4.4) 
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where 



x = x(r) = 



a(r) 
a 



was shown in Equation (3.9) to have the solution, 



(4.5) 



a(r) a 



a 



1 / ajv(r)' 



^4 



a 



+ 



'ajv(r)' 



a 



(4.6) 



written in terms of the Newtonian acceleration of a test particle at a separation, r, 



ajv(r) 



G N M(r) 



(4.7) 



where M(r) is the baryonic mass integrated within a sphere of radius, r. Each of the galaxy 
rotation curves in §4.2.3 are fitted in MOND by substituting the MOND acceleration law 
of Equation (4.6) into Equation (4.25) for the orbital velocity. 

Using the interpolating function of Equation (4.4), Sanders and McGaugh (2002) sug- 
gested that using the fits to the rotation curves of §4.2, using a revised cluster distance 
of 18.6 Mpc to Ursa Major, from the Cepheid-based re-calibrated Tully-Fisher relation of 
Sakai et al (2000), would imply that the MOND universal acceleration should be reduced 
to 



a 1.0 x 10 cm s" 



(4.8) 



In §4.2, ao is permitted to vary across the sample of 29 galaxies, in order to compute the 
MOND universal acceleration parameter, in Table 4.3, with the best-fitting results, 



a = < 



(1.34 ± 0.66) x lCT 8 cm s" 2 HSB galaxy subsample, 
(1.02 ± 0.78) x 10 -8 cm s" 2 LSB galaxy subsample. 



(4.9) 



Because of the gross uncertainty in the mean results, the galaxy rotation curves of Fig- 
ure 4.1 are one parameter best-fits by the stellar mass-to-light ratio, T, applying Equa- 
tion (4.8) universally. 
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For galaxies of sufficiently high surface brightness, the asymptotic circular velocity, 

v 4 out = a G N M, (4.10) 

satisfies the empirical Tully-Fisher relation L oc v^ ut provided one uses a luminosity param- 
eter which is proportional to the observed mass. This is shown more precisely in §4.2.8, 
where the stellar mass-to-light ratio, T, is the single free parameter, treated as a constant 
within a galaxy (recovering the empirical Tully-Fisher relation), but varying from galaxy 
to galaxy depending on the best-fit, with results tabulated in Table 4.3. 

Observations 

The MOND computations using HI and K-band photometric surface luminosity data, 
detailed in §4.2, with galaxy rotation curves plotted in Figure 4.1, indicate the following: 

1. The sample may be universally fit with a single MOND interpolating function and 
MOND acceleration constant given by Equations (4.4) and (4.8) yielding a best-fit 
stellar mass-to-light ratio T providing excellent to poor fits. 

2. The sample may be fit with a single MOND interpolating function but a varying, 
best-fit MOND acceleration parameter, tabulated in Table 4.3, yielding stellar mass- 
to-light ratios closer to T ~ 1, and correcting those poor fits with the universal 
MOND acceleration, but providing minor correction to those fits that were already 
good. The best-fit MOND acceleration parameter is not correlated with the galactic 
surface brightness. 

3. Every galaxy from the highest to lowest in surface brightness has a central disk that 
is dominated by the Newtonian potential, where the MOND interpolating function 
remains in the Newtonian core, \± ~ 1, and a MOND regime where ji > 1 outside of 
the core. 

4. Once within the MOND regime, the dynamics within the galactic disk continue to 
dominate, rising monotonically with orbital distance, as shown in Figure 4.4 which 
plot T(r) = vs. r. Unless \i is bounded from below, the dynamical mass factor, 
r(r) — > oc suggesting a classical instability. Conversely, since for every galaxy in the 
sample r(r out < 10, there is evidence that ji > ji^ ~ 0.1 is bounded by a cosmological 
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lower limit. 

5. The best-fit stellar mass-to-light ratio, T, is generally too large in the Newtonian 
core and too small in the MOND regime for those HSB galaxies that show poor fits, 
however the trend is reversed in those LSB galaxies that MOND does not fit well, 
as evident in Figure 4.6. For the lowest surface brightness galaxies, the increased 
uncertainty in the stellar mass-to-light ratio in the MOND regime lead to dramatic 
increases in the uncertainty in the total galaxy mass and relatively weak fits in the 
Newtonian core. 

6. Orphan features are traceable to a parent in either the gaseous or luminous stellar 
disks and are generally independent of choice of either the best-fit or universal accel- 
eration parameter, but become increasingly pronounced toward the outermost radial 
data point in the velocity rotation curve. 

4.1.3 Moffat's modified gravity 

Moffat's modified gravity theory predicts that galaxy rotation curves are explained by the 
radial acceleration law, 

o(r) _ _G(rj£M (4 . u) 

where G(r) is the effective gravitational constant, and varies through the galaxy. In the 
cores of each galaxy, where Newtonian gravity dominates the dynamics, G(r) ~ G/v, the 
value of Newton's constant. However, within a few kiloparsecs away from the core, the 
repulsive Yukawa forces becomes appreciable, G(r) > Gjy. For the analysis of galaxy 
rotation curves, we will consider the effective gravitational constant given in §3.3.5 by 
Equation (3.91), derived from MSTG (Moffat, 2005): 



G(r) = G N {1 + 



M(r) 



T 

1 - exp(-r/r ) ( H 

r 



(4.12) 



Brownstein and Moffat (2006a) applied the MSTG acceleration law of Equations (4.11) 
and (4.12) to a large sample of LSB and HSB galaxy rotation curves, obtaining satisfactory 
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fits with the parameters 

M = 96.0 x 1O 1O M , r = 13.92 kpc. (4.13) 

However, the dwarf LSB and HSB galaxy rotation curves were better fit with smaller values 
for these parameters. In §4.2, M and r are permitted to vary across the sample of 29 
galaxies, in order to compute the MSTG mean-universal parameters, in Table 4.4, with 
the overall best-fitting result, 

M = (98.6 ± 21.6) x 1O 1O M , r = (16.4 ± 6.1) kpc. (4.14) 

Figure 4.1 are subsequently one parameter best-fits by the stellar mass-to-light ratio, T, 
applying the mean results of Equation (4.14) universally. 

Moffat (2006b) introduced the scalar-tensor-vector gravity theory (STVG) where the 
weak-field, massive skew symmetric sector of NGT and MSTG are reduced to the simplest 
representation of the Yukawa meson - a massive Maxwell-Proca spin-1 - vector field, M 
and a triplet of scalar fields, The STVG modified acceleration law results from 

coupling the additional degrees of freedom to the Einstein metric, where Newton's constant 
and the Yukawa meson's coupling and range are dynamical scalar fields. 

The predictions of STVG mimics those of MSTG and NGT at astrophysical scales, 
but since the basic excitations of the three theories are qualitatively different, fits to astro- 
physical phenomena may constrain the phenomenological parameter space. Moreover, since 
STVG is a relatively simple gauge theory of gravitation, the static, spherically symmetric 
solution has been calculated exactly and resembles the Reissner-Nordstrom solution, but 
with the source electromagnetic charge replaced by the source "fifth force" charge between 
fermions and the massive Maxwell-Proca spin-1 - vector field. 

A derivation of a new acceleration law in STVG - from the action principle, but without 
necessary ad-hoc phenomenological input - provided a modified acceleration law of the 
form of Equation (4.11), where the effective gravitational constant is determined from the 
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modified central force law, given in §3.4.5 by Equations (3.159), (3.160) and (3.162): 

G(r) = G N {l + a-ae-^ r (l + fir)} , (4.15) 

(yjM(r) + E) VGw 7 

/i(r) = -j§= (4.17) 
^M(r) 

obtaining satisfactory fits with universal parameters, (Moffat and Toth, 2009) 

D = 6.25Mq /2 /pc, E = 25000Mq /2 , = 20G N . (4.18) 

In §4.2, D, E and are permitted to vary across the sample of 29 galaxies, in order to 
compute the STVG mean-universal parameters, in Table 4.5, with the overall best-fitting 
result, 



D = (6.44 ± O.2O)0W©pc-\ E = (28.4 ± 7.9) x 10 3 Mq /2 , = (24.4 ± 18.0)G N . 

(4.19) 

The galaxy rotation curves of Figure 4.1 are subsequently one parameter best-fits by the 
stellar mass-to-light ratio, T, applying the mean results of Equation (4.19) universally. 

Observations 

The MSTG and STVG computations using HI and K-band photometric surface luminosity 
data, detailed in §4.2, with solutions plotted in Figure 4.1, indicate the following: 

1. The sample may be universally fit with the MOG acceleration law given by Equations 
(4.11) and (4.12) in MSTG, or Equations (4.15), (4.16) and (4.17) in STVG, and 
either universal MSTG or universal STVG parameters of Equations (4.14) and (4.19), 
respectively, yielding best-fit stellar mass-to-light ratios, T, providing fits and galaxy 
masses comparable to MOND. 

2. The sample may be fit with varying, best-fit MSTG parameters of Equation (4.12), 
tabulated in Table 4.4, yielding stellar mass-to-light ratios closer to T ~ 1 for the HSB 
galaxies and farther from unity for the LSB galaxies T > 1, and correcting the fits 
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using the mean-universal MSTG parameters. The best-fit MSTG scale parameters 
are correlated with the galactic surface brightness, with larger values of M and 
smaller values of r preferred by HSB galaxies, with the reverse trend in LSB galaxies. 

3. The sample may be fit with varying, best-fit STVG parameters of Equation (4.15), 
tabulated in Table 4.5, yielding the most robust stellar mass-to-light ratios T ~ 1 of 
all the gravity theories tested, and correcting the fits with the mean-universal STVG 
parameters. The best-fit STVG parameters are strongly correlated with the galactic 
surface brightness, with larger values preferred by HSB galaxies, and smaller values 
preferred by LSB galaxies, with the greatest uncertainty in the STVG cosmological 
parameter, G^. 

4. Every galaxy from the highest to lowest in surface brightness has a central disk that 
is dominated by the Newtonian potential, where the MOG running gravitational 
coupling remains in the Newtonian core, G(r) ~ G^, and a MOG regime where 
G(r) > Gn outside of the core. 

5. Once within the MOG regime, the dynamics within the galactic disk continue to 
dominate, rising monotonically with orbital distance, as shown in Figure 4.4 which 
plot T(r) = G(t)/Gn vs. r. Since G(r) is bounded from above, the dynamical mass 
factor, r(r) — > confirming the asymptotic stability of MOG. This is the primary 
difference between the MOG and MOND predictions, in which the former return 
to a Keplerian behaviour (with larger than Newton gravitational coupling) whereas 
MOND favours asymptotically flat galaxy rotation curves until the dynamics are 
correlated with other systems. 

Lessons learned from modified gravity theory are presented in §7.1.2 of the summary 
chapter of Part III. 

4.2 Ursa Major filament of galaxies 

Surrounding the local group of the Milky Way, the Coma-Sculptor cloud is our home in 
the Virgo supercluster. Tully et al (1996) identified 79 galaxies of the Ursa Major (UMa) 
filament, in the first of a trilogy of works, and provided surface brightness measurements 
in the blue, red and infrared bands. UMa lies in the plane of the Virgo supercluster at 
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the junction of filamentary structures, beyond the long axis of the filament of the Coma- 
Sculptor cloud. 

In a sequel, Tully and Verheijen (1997) explored the differences in surface brightness 
amongst the 62 galaxies of the complete sample, and identified two distinct radial config- 
urations of spiral galaxies of varying size, all of which are unevolved and rich in HI gas 
consistent with observations in Virgo and Fornax. These are the high surface brightness 
(HSB) galaxies and the low surface brightness (LSB) galaxies. Remarkably, even though 
the configurations differ, both types of HSB and LSB galaxies exhibit a common expo- 
nential disk for the central surface brightness profiles. However, the mass-to-light ratio in 
LSB galaxies is difficult to explain using Newtonian gravity without dark matter dominated 
cores, whereas the mass-to-light ratio in HSB galaxies is difficult to explain using Newto- 
nian gravity without baryon dominated cores and extended dark matter halos. The UMa 
sample suggests that structure formation avoids the region of parameter space between 
LSB and HSB galaxies, possibly due to different angular momentum regimes. Passing 
from high to low specific angular momentum, there is first the transition from LSB to 
HSB regimes, and at very low specific angular momentum, there is another transition from 
HSB galaxies which are exponential disk dominated to disk and bulge dominated. Since 
the first transition can be modelled by a single parameter - the mass-to-light ratio of the 
stellar exponential disk, whereas the second transition requires a second parameter - the 
mass-to-light ratio of the bulge, this chapter will focus exclusively on the subsample where 
the bulge can be neglected; and one parameter fits are possible for both HSB and LSB 
galaxies. 

Sanders and Verheijen (1998) presented a third paper in the series on the rotation curves 
of UMa galaxies which focussed on the near-infrared band, because it is relatively free of 
the effects of dust absorption and less sensitive to recent star formation. In this work, it 
was observed that the exercise of fitting dark matter halos to galaxy rotation curve data 
required at least three free parameters per galaxy (stellar disk mass-to-light ratio, halo 
core radius and density normalization) and essentially any observed rotation curve can be 
reproduced. However, dark matter gained more predictive power when the density law 
was parametrized by singular r -7 halos, with 1 < 7 < 3. However, although these singular 
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halo models such as Navarro et al (1996) produced acceptable fits to HSB galaxies, they 
generally failed for LSB galaxies. McGaugh and de Blok (1998) tested the dark matter 
hypothesis with LSB galaxies, finding that progressively lower surface brightness galaxies 
have progressively larger mass deficits, requiring high concentrations of dark matter deep 
in the galaxy core, rendering the visible components insignificant to the galaxy dynamics 
and leading to fine-tuning problems. In comparison, MOND was shown to fit all of the 
galaxies in the sample with only a single free parameter (disk mass-to-light ratio) although 
MOND itself has the Milgrom acceleration parameter, ag and the best-fit mass-to-light 
ratio also depends on the choice of a universal interpolating function fi(x). The notion 
of fitting galaxy rotation curves without dark matter was further explored in Brownstein 
and Moffat (2006a) utilizing a larger sample (including UMa) where it was confirmed that 
MOND provided good one parameter fits to the sample's galaxy rotation curves with a 
universal choice of a and /i(x). 

4.2.1 Photometry 

According to Sanders and Verheijen (1998), the existence of X-band surface photometry 
is a great advantage since the near-infrared emission, being relatively free of the effects of 
dust absorption and less sensitive to recent star formation, is a more precise tracer of the 
mean radial distribution of the dominant stellar population. The principal advantages of 
using infrared luminosities is that stellar mass-to-light ratios are less affected by population 
differences and extinction corrections are minimal (Verheijen, 2001). 

The galaxy rotation curves of §4.2 are divided into high and low surface brightness 
galaxies , as in §4.2.3. The component velocities plotted in Figure 4.1 are based on the 
surface photometric data of the gaseous disk (HI plus He) component and luminous stellar 
disk component. The method of generating the rotation curves closely followed Sanders 
and Verheijen (1998) and (Verheijen, 2001). The ROTMOD task of that group's Gronin- 
gen Image Processing System (GIPSY) 1 was used to analyse the HI and X-band surface 
photometry data to produce the velocity profiles of the gaseous disk (HI plus He) distri- 
bution and luminous stellar disks, accounting for the Verheijen and Sancisi (2001) revised 

1 http : //www . astro . rug . nl/~gipsy/ 
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distance estimate to UMa from D = 15.5 Mpc to D = 18.6 Mpc. 

For each galaxy in the UMa sample, the photometric data is best-fit to the galaxy 
rotation curve data through a non-linear least squares fitting algorithm, which minimizes 
the weighted sum of squares of deviations between the fit and the data. The sum of 
squares of deviations is characterized by the estimated variance of the fit. The reduced 
X 2 jv statistic is computed as the value of x 2 per degree of freedom, 



where the number of degrees of freedom, z/, is the difference between the number of data 
points in the galaxy rotation curve, N ) and the number of free parameters, p. 

4.2.2 Surface mass computation 

The gaseous disk is modelled as an infinitely thin, uniform disk and the surface mass density 
profile is derived numerically by means of a computation allowing for high resolution 
sampling of the HI gas data. The UMa sample was resolved at the sub-kiloparsec scale, 
equivalent to a resolution of 



where r out is the outermost observed radial position, measured in kiloparsecs, of the rotation 
velocity data, listed for each galaxy in Column (4) of Table 4.1; and Column (5) is the 
observed velocity at the outermost observed radial position. 
The surface brightness computation, 



derives the absolute surface mass density of the HI gas, where / Sdv is the integrated 
HI flux density in units of Jy km/s as measured from the global HI profile - taken from 
Column (15) of Table 2 of Verheijen (2001), and D is the distance in Mpc. The computation 
results in a radial surface mass profile, at the resolution, g ) of Equation (4.21), and a total 
(integrated) result which is an absolute measurement of the HI disk mass. Although the 



reduced — x? = X 2 /^> where v = N — p ) 



(4.20) 




(4.22) 
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computation is free of unspecified parameters, it is systematically affected by changes in 
distance estimates. 

However, since all of the UMa galaxies are located within a filament of the Coma- 
Sculptor cloud - and at similar redshift, listed in Column (3) of Table 4.1 - luminosity 
distances are common and the uncertainty in the mass-to-light ratio is greatly reduced. 
This improves the certainty in identifying orphan features, as in §4.3.2, which appear from 
emergent surface mass profiles. The total (integrated) HI gas masses, determined by the 
computation, are listed in Column (2) of Table 4.6, and the total mass of the HI (and He) 
gaseous disk is determined by scaling the HI gas mass by the Big bang nucleosynthesis 
Helium fraction, 

M gas = ^M HI) (4.23) 

where the 4/3 BBN He fraction is enforced across the UMa sample. This introduces a 
margin of uncertainty, in Equation (4.23), which increases radially due to evolutionary 
changes in the distribution of HI and He, since the formation of the UMa filament of 
galaxies. Hoekstra et al (2001) showed that the BBN scale factor would have to increase 
by a factor of ~ 7 to fit a sample of 24 spiral galaxies without dark matter, obtaining 
good fits for most galaxies, but not for those galaxies which show a rapid decline of the HI 
surface density in the outermost regions. 

The luminous stellar disk was assumed to be described by the Van der Kruit and Searle 
law, where the disk density distribution as a function of z (vertical height from the plane 
of the disk) is given by 

£(*) = sech 2 (z/z )/z , (4.24) 

where z is the vertical scale height of the luminous stellar disk, and was assumed to 
be 20% of the near infrared exponential disk scale length according to Column (13) of 
Table 2 of Verheijen (2001), as listed in Column (4) of Table 4.1. The surface brightness 
computations using Equation (4.24) return the surface mass density of the stellar disk to 
within an overall multiplicative factor, T, which is strictly set to unity in the computation, 
as listed in Column (3) of Table 4.6. The idea of a varying the stellar mass-to-light ratio, 
T(r) throughout a galaxy - shown in Figure 4.6 - would lead to perfect fits for any gravity 
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Table 4.1: Galaxy properties of the Ursa Major sample 
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Notes. — Relevant galaxy properties of the UMa sample: Column (1) is the NGC/UGC galaxy 
number. Column (2) is the galaxy morphological type. Column (3) is the observed redshift 
from the NASA/IPAC Extragalactic Database. Column (4) is the K-band vertical scale height 
of the luminous stellar disk, and Column (5) is the K-band luminosity data converted from the 
2MASS K-band apparent magnitude via Equation (4.58). Column (6) is the velocity amplitude 
(maximum) of the rotation curve. Column (7) is the outermost observed radial position in the 
rotation velocity data; and Column (8) is the observed velocity at the outermost observed position. 
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theory, but the arbitrariness of such a solution would lead to a fine-tuning problem, and 
instead a best-fit T is computed by a nonlinear least-squares algorithm - shown in the 
same figure - with results provided for the best-fit NFW and core-modified dark matter 
profiles in Table 4.2, the best-fit MOND universal acceleration in Table 4.3, and the best-fit 
MSTG and STVG parameters in Table 4.4 and Table 4.5, respectively. 

4.2.3 High and low surface brightness galaxies 

The galaxy rotation curves, in Figure 4.1, plot the rotation velocity profiles, 



in km s x , vs. r in kpc, where the acceleration law, a(r), is given by Equations (2.1) and 
(4.7) for Newton's theory (with and without dark matter), Equation (4.6) for Milgrom's 
MOND, and Equation (4.11) for Moffat's MOG. 

Shown for each galaxy are the mean-universal best-fits according to Moffat's STVG and 
MSTG theories and Milgrom's MOND; and the best-fitting core-modified dark matter - 
and the corresponding core-modified dark matter halo components. The best-fit Newtonian 
results (visible baryons only) are plotted for comparison. 

Newtonian core 

Each of the gravity theories which fit high and low surface brightness galaxies disagree 
with Newton's theory without dark matter. The disagreement is small in the core of each 
galaxy, but increases with separation from the center. In each galaxy in the sample, there 
is a Newtonian core where neither modified gravity nor dark matter is required to fit the 
galaxy rotation curves. The computation of the radius of the Newtonian core weights 
the velocity points inside the core, and discards the velocity points outside the core, and 
yields a single parameter, the best-fit stellar mass-to-light ratio, T, for each galaxy, as 
detailed in §4.2.8. For every galaxy in the sample, the best-fit Newtonian core model, 
plotted in brown dot-dotted lines, shows the characteristic Keplerian behaviour outside 
the Newtonian core, which disagrees with the galaxy rotation curves, but the model shows 
reduced x 2 jv comparable to the modified gravity theories within the Newtonian core. 




(4.25) 
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Figure 4.1: UMa — Rotation curves. 



The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc for 19 
HSB and 10 LSB galaxies. The dynamic data consist of the measured orbital velocities. 
The photometric data sets consist of the actual HI gas component and the stellar disk 
component, with a normalized stellar mass-to-light ratio, T = 1. The computed best-fit 
results by varying the stellar mass-to-light ratio, T, are plotted for Moffat's STVG 
and MSTG theories and Milgrom's MOND theory with mean-universal parameters. 
Results are plotted for the best-fit core-modified dark matter theory including visible 
baryons, and the corresponding dark matter halo component. The best-fit Newtonian 
core model (visible baryons only) is plotted for comparison. The figure is continued. 

Best-fit and universal parameters 

Other than Newton's theory without dark matter, each of the gravity theories which fits 
the sample of high and low surface brightness galaxies contain either one, two or three 
parameters unique to their respective acceleration laws. In the case of dark matter, as 
discussed in §4.1.1, either the NFW profile of Equation (4.1) or the core-modified profile 
of Equation (4.2) contain two parameters - the central dark matter density, p , and the 
scale radius, r s - which are taken to vary from galaxy to galaxy. 

Conversely, the hope for the modified gravity theories is that there exist universal 
parameters which simultaneously fit all of the data for high and low surface brightness 
galaxies. To research this possibility, each of the parameters unique to each of the modified 
gravity theories and the stellar mass-to-light ratio, T, were first allowed to vary and a table 
of best-fit values was constructed. Then each parameter was averaged across the 
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Figure 4.1 continued: UMa — Rotation curves. 
The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc. 
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Figure 4.1 continued: UMa — Rotation curves. 
The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc. 
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Figure 4.1 continued: UMa — Rotation curves. 
The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc. 
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Figure 4.1 continued: UMa — Rotation curves. 
The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc. 
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Figure 4.1 continued: UMa — Rotation curves. 

The rotation velocity profile, v(r) in km s _1 , vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The dynamic data consist of the measured orbital velocities. The photometric 
data sets consist of the actual HI gas component and the stellar disk component, with a 
normalized stellar mass-to-light ratio, T = 1. The computed best-fit results by varying 
the stellar mass-to-light ratio, T, are plotted for Moffat's STVG and MSTG theories and 
Milgrom's MOND theory with mean-universal parameters. Results are plotted for the 
best-fit core-modified dark matter theory including visible baryons, and the corresponding 
dark matter halo component. The best-fit Newtonian core model (visible baryons only) is 
plotted for comparison. 
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subsamples of HSB and LSB galaxies, and across the complete sample, and mean-universal 
parameters were chosen. Next the galaxy rotation curves were refitted and replotted 
with one free parameter, the stellar mass-to-light ratio, T, using the tabulated universal 
parameters. Therefore all of the dark matter fits have three free parameters, whereas all of 
the modified gravity theories have only one free parameter, and should be compared with 
this in mind. 

Overall, the core-modified dark matter model shows the lowest reduced xV^ statistic, 
but the model is less predictive than MOND, MSTG, or STVG due to the variation across 
the sample in the p and r s parameters in the dark matter fitting formulae. The only 
theory that fails to produce x 2 best-fits for some dwarf galaxies was the NFW model for 
NGC 4389 and UGC 6818, as shown in Table 4.2. 

When MOND as in Table 4.3, MSTG as in Table 4.4, and STVG as in Table 4.5, 
use the best-fit parameters instead of the mean-universal parameters, the reduced \ 2 l v 
statistic of Equation (4.20) decreases, but there is no case where using mean-universal 
parameters leads to disagreement with the photometry. It is reasonable to conclude that 
MOND, MSTG and STVG provide acceptable fits to galaxy rotation curves with universal 
parameters and variable mass-to-light ratios, as shown in Figure 4.1. 

4.2.4 Surface mass density maps 

The importance of being able to determine the distribution of matter in astronomical 
objects is that it allows predictions for ongoing and future experiments, such as galaxy- 
galaxy lensing, which measures the line-of-sight surface mass density, 




(4.26) 



through the convergence, 



k{t) 



S(r) 



(4.27) 



where 



c 2 D s 



(4.28) 
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Figure 4.2: UMa — Surface mass densities. 



The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc for 19 HSB and 10 
LSB galaxies. The photometric data sets consist of the actual HI gas component and the 
stellar disk component, with a normalized stellar mass-to-light ratio, T = 1. The com- 
puted best-fit results by varying the stellar mass-to-light ratio, T, are plotted for Moffat's 
STVG and MSTG theories and Milgrom's MOND theory with mean-universal param- 
eters. Results are plotted for the best-fit core-modified dark matter theory including 
visible baryons, and the corresponding dark matter halo component. The best-fit Newto- 
nian core model (visible baryons only) is plotted for comparison. The figure is continued. 

is the Newtonian critical surface mass density (with vanishing shear), D s is the angular 
distance to the source, background galaxy, D\ is the angular distance to the lens, foreground 
galaxy. (Peacock, 2003, Chapter 4). 

The E-maps plotted in Figure 4.2 provide high resolution sub-kiloparsec predictions, 
whereas the current state of the art in galaxy-galaxy lensing yield only course grained obser- 
vations, with resolutions of < 10 kpc/pixel. The MOG predictions for future high resolution 



ft-maps must account for the modified acceleration law of Equation (4.11) Brownstein and 
Moffat (2007): 




(4.29) 



where 




(4.30) 
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Figure 4.2 continued: UMa — Surface mass densities. 

The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc. 
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Figure 4.2 continued'. UMa — Surface mass densities. 

The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc. 
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Figure 4.2 continued'. UMa — Surface mass densities. 

The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc. 
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Figure 4.2 continued: UMa — Surface mass densities. 

The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc. 
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Figure 4.2 continued: UMa — Surface mass densities. 

The surface mass density, S(r) in M /pc 2 , vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The photometric data sets consist of the actual HI gas component and the stellar 
disk component, with a normalized stellar mass-to-light ratio, T = 1. The computed best- 
fit results by varying the stellar mass-to-light ratio, T, are plotted for Moffat's STVG and 
MSTG theories and Milgrom's MOND theory with mean-universal parameters. Results 
are plotted for the best-fit core-modified dark matter theory including visible baryons, 
and the corresponding dark matter halo component. The best-fit Newtonian core model 
(visible baryons only) is plotted for comparison. 
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is the weighted surface mass density, and S c is the usual Newtonian critical surface mass 
density Equation (4.28). Moffat and Toth (2008) simplified Equations (4.29) and (4.30), 
in STVG, in the case that the lens may be treated as a point source, but not for extended 
mass distributions relevant for galaxy-galaxy lensing. 

The surface mass density due to the visible component is, 

^ bar yon (0 Sg as (r) + TS disk (r), (4.31) 

and therefore the E-map computed for each galaxy depends on the best-fitting stellar 
mass-to-light ratio, T, determined separately for each gravity theory: 



E(r) 



Sbaryon(^) Newtonian core, MOND 

S baryon (r) MSTG, STVG (4.32) 

Sbaryon(r) + S ha i (r) Dark matter. 



4.2.5 Radial mass profiles 

The missing mass problem is best visualized by solving Equation (4.25) for the Newtonian 
dynamic mass of Equation (4.7), 



M N (r) = 




(4.33) 



where the velocity points are dynamic variables determined from the galaxy rotation curves 
of Figure 4.1. 

The visible component mass profiles, plotted in Figure 4.3, are based on the surface 
photometric data of the gaseous disk (HI plus He) component and luminous stellar disk 
component, 

M gas (r) = ^ r 2^r'E gas (r0dr', (4.34) 
M disk (r) = f r 27rr f E dlsk (r f )dr f . (4.35) 



54.2. URSA MAJOR FILAMENT OF GALAXIES 



105 



M(r) 


red crosses with error bars 


Component 


lines 


Actual HI gas 


green dashed line 


Stellar disk [T = 1] 


magenta dotted line 


Dark matter halo 


orange fine-dotted line 


Theory- 


lines 


Newtonian core 


brown dot-dotted line 


STVG 


black solid line 


MSTG 


blue short dash-dotted line 


MOND 


cyan long dash-dotted line 


Dark matter 


orange solid line 




100 



HSB (a) NGC 3726 



Figure 4.3: UMa — Mass profiles. 



The radial mass profile, M(r) in M , vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The dynamic data consist of the Newtonian dynamic mass due to the measured 
orbital velocities. The photometric data sets consist of the actual HI gas component and 
the stellar disk component, with a normalized stellar mass-to-light ratio, T = 1. The com- 
puted best-fit results by varying the stellar mass-to-light ratio, T, are plotted for Moffat's 
STVG and MSTG theories and Milgrom's MOND theory with mean-universal param- 
eters. Results are plotted for the best-fit core-modified dark matter theory including 
visible baryons, and the corresponding dark matter halo component. The best-fit Newto- 
nian core model (visible baryons only) is plotted for comparison. The figure is continued. 

The integrated mass profile due to the visible component is therefore, 



M baryon (r) = M gas (r) + TM disk (r), 



(4.36) 



which depends on the best-fitting stellar mass-to-light ratio, T, for each galaxy, determined 
separately for each gravity theory. 



M(r) = 



^baryon(^) Newtonian core, MOND, MSTG, STVG 

M baryon (r) + M ha i (r) Dark matter, 



(4.37) 



where the dark matter halo may be computed according to either the NFW formula of 
Equation (4.1), or alternatively the core-modified formula of Equation (4.2). 
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Figure 4.3 continued: UMa — Mass profiles. 
The radial mass profile, M(r) in M , vs. orbital distance, r in kpc. 
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Figure 4.3 continued: UMa — Mass profiles. 
The radial mass profile, M(r) in M , vs. orbital distance, r in kpc. 
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Figure 4.3 continued: UMa — Mass profiles. 
The radial mass profile, M(r) in M , vs. orbital distance, r in kpc. 
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Figure 4.3 continued: UMa — Mass profiles. 
The radial mass profile, M(r) in Mq, vs. orbital distance, r in kpc. 
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Figure 4.3 continued: UMa — Mass profiles. 



The radial mass profile, M{r) in Mq, vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The dynamic data consist of the Newtonian dynamic mass due to the measured 
orbital velocities. The photometric data sets consist of the actual HI gas component and 
the stellar disk component, with a normalized stellar mass-to-light ratio, T = 1. The com- 
puted best-fit results by varying the stellar mass-to-light ratio, T, are plotted for Moffat's 
STVG and MSTG theories and Milgrom's MOND theory with mean-universal parame- 
ters. Results are plotted for the best-fit core-modified dark matter theory including visible 
baryons, and the corresponding dark matter halo component. The best-fit Newtonian core 
model (visible baryons only) is plotted for comparison. 
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The total mass of each galaxy due to the visible components is 

M baryon = M gas + TM disk , (4.38) 

and therefore the total mass of each galaxy depends on the best-fitting stellar mass-to-light 
ratio, T, determined separately for each gravity theory: 



Mi 



total 



M barvon Newtonian core, MOND, MSTG, STVG 

(4.39) 

M baryon + M ha i Dark matter, 



where the dark matter halo mass may be computed according to Equation (4.1) according 
to the NFW fitting formula, or alternatively according to the core-modified fitting formula 
of Equation (4.2). The final results for the total galaxy masses, according to the best-fitting 
stellar mass-to-light ratio, T, for each gravity theory are provided in Table 4.6. 

Every galaxy studied, from the highest to lowest in surface brightness, from the most 
giant to the smallest dwarf, exhibit Newtonian dynamic masses far in excess of the mass 
profiles due to the visible components, outside the Newtonian core. The situation within 
the Newtonian core depends on the particular gravity theory being applied. Milgrom's 
theory provides a region inside the MOND regime where accelerations are larger than 
a , where Moffat's MOG theories provide a region inside the MOG regime where the 
gravitational coupling G(r) ~ Gn- The core-modified dark matter halo is spherical, and 
does not dominate the visible disks until a critical radius is reached. In all cases, there is 
a transition region just outside the Newtonian core where either some form of dark matter 
is required, or some modification of gravity provides sufficient violations of the strong 
equivalence principle. 
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4.2.6 Dynamic mass factor 

Since the Newtonian dynamic mass greatly exceeds the baryonic mass outside the Newto- 
nian core, each gravity theory must make up the difference in order to fit the data. For 
dark matter, the difference is the halo component, 

M N (r) = M baryon + M halo . (4.40) 

The dynamic mass factor is defined as the Newtonian dynamic mass per unit baryonic 
mass 

r ( r )- Mjv(r) -1 i Mhalo(r) r 4 zm 

For Milgrom's MOND, the difference is due to the reciprocal factor of the smaller than 
unity MOND interpolating function, 

M N (r) = Mbar H (r) , (4-42) 
and the dynamic mass factor is denned as 

r — (r) s S = wr (4 ' 43) 

For Moffat's MOG, the difference is due to the multiplicative factor of the larger than 
Newton gravitational coupling, 

m„ W = gW^g-M (4 . 44) 

and the dynamic mass factor is defined as 

_ , , M N {r) G{r) lA A , 

r MOG (r) = - = (4.45) 

The dynamic mass factors, plotted in Figure 4.4 for each galaxy in the UMa sample, 



4.2. URSA MAJOR FILAMENT OF GALAXIES 



115 




r [kpc] 
HSB (a) NGC 3726 

Figure 4.4: UMa — Dynamical mass factors. 

The Dynamical mass factor, T(r), vs. orbital distance, r in kpc for 19 HSB and 10 
LSB galaxies. The dynamic data consist of the Newtonian dynamic mass due to 
the measured orbital velocities per unit baryonic mass per gravity theory, shown 
with error bars. The computed best-fit results by varying the stellar mass-to- 
light ratio, T, are plotted for Moffat's STVG and MSTG theories and Milgrom's 
MOND theory with mean-universal parameters. Results are plotted for the best-fit 
core-modified dark matter theory including visible baryons. The figure is continued. 

are shown to be a monotonically rising (near linear) functions, with similar properties: 

r(r) ~ 1 within Newtonian core, 

T(r) 3> 1 beyond Newtonian core, (4.46) 
T(r) < 10 within galaxy, 

and is the measure of the missing mass factor. Each theory may be judged by how well the 
Newtonian dynamic mass due to the measured orbital velocities per unit baryonic mass per 
gravity theory, shown with error bars, corresponds to the predictions of Equations (4.41), 
(4.43) and (4.45). The dynamic mass factor provides a unifying picture for dark matter 
and phantom dark matter and can be phenomenologically applied to constrain the choice 
of the MOND interpolating function - without ad hoc choices - and the form of Moffat's 
varying gravitational coupling. 
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Figure 4.4 continued: UMa — Dynamical mass factors. 

The Dynamical mass factor, I\r), vs. orbital distance, r in kpc. 
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Figure 4.4 continued: UMa — Dynamical mass factors. 

The Dynamical mass factor, I\r), vs. orbital distance, r in kpc. 
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Figure 4.4 continued'. UMa — Dynamical mass factors. 

The Dynamical mass factor, I\r), vs. orbital distance, r in kpc. 
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Figure 4.4 continued: UMa — Dynamical mass factors. 

The Dynamical mass factor, T(r), vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The dynamic data consist of the Newtonian dynamic mass due to the measured 
orbital velocities per unit baryonic mass per gravity theory, shown with error bars. The 
computed best-fit results by varying the stellar mass-to-light ratio, T, are plotted for 
Moffat's STVG and MSTG theories and Milgrom's MOND theory with mean-universal 
parameters. Results are plotted for the best-fit core-modified dark matter theory including 
visible baryons. 
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4.2.7 Core- modified dark matter halos 

The simple observation that galaxy rotation curves are approximately flat at large radii, 
where the orbital velocity of Equation (4.25) is constant, leads to the conclusion of Equation 
(2.2) that the Newtonian dynamic mass of Equation (4.33) grows linearly with radius, and 
therefore, since we are not neglecting baryons, 

P(r) = Pbaryon(r) + Phalo(0 OC r" 2 , (4.47) 

is valid where the galaxy rotation curves are approximately flat. However, the radial 
distribution of spherically averaged dark matter halos is unlike either of the baryonic 
components which accumulate in exponentially thin HI (and He) gaseous disks or luminous 
stellar disks of Equation (4.24), and 

Pbaryon(^) > Phaio(^) on the galactic plane, (4.48) 
Pbaryon(^) <C Phaio(^) of f the galactic plane, (4.49) 

4 I I 



3 - 
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Figure 4.5: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(r), vs. orbital distance, r in kpc for 19 HSB 
and 10 LSB galaxies. The logarithm slope, 7(7-) of the dark matter power- 
law, vs. orbital distance, r in kpc. The computed best-fit results are plot- 
ted for core-modified dark matter theory including visible baryons and 
the corresponding dark matter halo component. The figure is continued. 
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Therefore in order to generalize the power-law index of Equation (2.14) to include 
baryons, it is convenient to derive the relation using the spherically integrated power-law 
density of Equation (2.5), 

M N (r) oc jT ^dr', (4.50) 

and 

dM N {r) oc r 2 ~ l{r) dr. (4.51) 

Thus the spherically averaged power-law index may be defined in terms of the logarithm 
slope, 

dmr 

The power-law indices for the best-fit core-modified dark matter halo, given by Equation 
(2.14), and the Newtonian dynamic mass including baryons, according to Equation (4.52), 
are plotted in Figure 4.5. 

Since the virial radius of the halo naturally extends beyond the outermost radial point 
in the galaxy rotation curve, r out , the dark matter to baryon fraction can grow without 
bound until the cosmological limit is reached. Within each galaxy in the sample, the dark 
matter to baryon fraction is tabulated to the outermost radial point in the galaxy rotation 
curve - in Column (8) of Table 4.6 - with mean values: 



M ha i Q (r out ) 

-^baryon (^out ) 



2.4 ±2.1 HSB galaxies 

1.3 ±0.5 LSB galaxies (4.53) 
2.0 ±1.8 full sample. 



which are consistently below the upper limit set by Spergel et al (2007) in the Wilkinson 
microwave anisotropy probe (WMAP) third year results. 

Solution to the dark matter cusp problem 

The conflict between the cuspy dark matter halos predicted by N-body simulations and 
the constant density cores preferred by dwarf and low surface brightness galaxies may be 
resolved by a universal core-modified fitting formula with a constant density core, while 
including the visible baryons which are dominant in the galaxy core. 
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Figure 4.5 continued: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(r), vs. orbital distance, r in kpc. 
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Figure 4.5 continued: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(7*) , vs. orbital distance, r in kpc. 
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Figure 4.5 continued: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(7™) , vs. orbital distance, r in kpc. 
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Figure 4.5 continued: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(7*) , vs. orbital distance, r in kpc. 
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Figure 4.5 continued: UMa — Dark matter power-law logarithm slopes. 

The logarithm slope profile, 7(7™), vs. orbital distance, r in kpc for 19 HSB and 10 LSB galaxies. 
The logarithm slope, 7(7*) of the dark matter power-law, vs. orbital distance, r in kpc. 
The computed best-fit results are plotted for core-modified dark matter theory including 
visible baryons - and the corresponding dark matter halo component. 
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As the plot of the dark matter power-law proves in Figure 4.5, at large distances from 
the center of each galaxy in the sample, the density profile of the dark matter halo is well 
described by a steep power-law, with power-law index 7^3, whereas at distances toward 
the center of the galaxy an increasingly shallow power-law is observed. For distances less 
than the dark matter halo core radius, r < r s ,the total density profile including baryons 
shows a universal 7 — ► 1 power-law index, and the density profile of the dark matter 
component alone approaches a rarified, constant density core. 

A comparison of Table 4.2, show a statistically significant reduction of the x 2 l v test in 
^ 90% of the galaxies results from using the core-modified profile of Equation (4.2) instead 
of the NFW profile of Equation (4.1). Moreover, in those galaxies that the singular NFW 
profile fits well, the best-fit stellar mass-to-light ratio, T <C 1, which has prompted the dark 
matter community of physicists to disregard the baryonic component in their simulations. 
Most strikingly, in one HSB and one LSB galaxy, there does not exist a best-fit NFW profile 
with any nonzero stellar mass-to-light ratio, forcing T = for these galaxies. Alternatively, 
the core-modified profile prefers values for the best-fit stellar mass-to-light ratio, T ~ 1, 
which are physically acceptable for every galaxy in the sample. 

This core-modified dark matter galaxy model produces excellent fits to the galaxy 
rotation curves of Figure 4.1, and enables predictions of detailed surface mass density 
maps, as shown in Figure 4.2, and demonstrates excellent fits to the mass profiles of 
Figure 4.3, with dark matter to baryon fractions consistent with cosmologically observed 
values. 

4.2.8 The mass luminosity relationship 

Throughout this work, the stellar mass-to-light ratio, T, is treated as a free parameter, 
with results near unity considered reasonable. Each gravity theory which attempts to fit 
the galaxy rotation curve to the integrated surface mass densities of the HI (and He) gas 
and stellar disk components will select a best- fit stellar mass-to-light ratio, for each galaxy. 

The behaviour of the mass-to-light ratio within each galaxy and the change in the 
behaviour from galaxy to galaxy are important concerns of a good fit. For the sample of 
galaxies considered in §4.2, the basic computation is that of the surface mass computation 
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Figure 4.6: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc for 19 HSB and 10 
LSB galaxies. The stellar mass-to-light ratio, T, required to fit the galaxy rotation 
curve at each data point without dark matter is plotted for Moffat's STVG and MSTG 
theories, Milgrom's MOND theory, and Newton's theory. Within each theory, the 
best-fit value of the stellar mass-to-light ratio, T, is shown with a horizontal line, 
including the best-fit core-modified dark matter; and the best-fit Newtonian core 
model - the extent of the core shown with a vertical line. The figure is continued. 

of §4.2.2 of the individually detected components: the exponentially thin gaseous (HI and 
He) disk of Equations (4.22) and (4.23) and the luminous stellar disk of Equation (4.24), 
with the bulge neglected for the reasons stipulated at the start of this section in regards 
to Tully and Verheijen (1997). 

In order to calculate the total mass of each galaxy from the photometry, within the 
context of each gravity theory, the mass luminosity relationship is based on the best- 
fitting stellar mass-to-light ratio, T, according to Equation (4.38). The HI (and He) 
gaseous component is the only computation that is independent of assumptions on the 
mass-to-light ratio and is determined by big bang nucleosynthesis according to Equation 
(4.23). 

The best-fit stellar mass-to-light ratio, T, and the computed total galaxy mass is listed 
for each galaxy in Table 4.6. It is clear that the best-fit mass-to-light ratio varies from 
galaxy to galaxy; and none of the galaxy rotation curves in the UMa sample can be 
fit by a universal-mean stellar mass-to-light ratio. The possibility that the stellar mass- 
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Figure 4.6 continued: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc. 
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Figure 4.6 continued: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc. 



20 



T(r) 


Best fit 


with horizontal lines 


Mass to light profile 


with lines 


Newtonian core 


brown dot-dotted horizontal 


Newton's theory 


brown dot-dotted line 


STVG 


black solid horizontal 


STVG 


black solid line 


MSTG 


blue short dash-dotted horiz. 


MSTG 


blue short dash-dotted line 


MOND 


cyan long dash-dotted horiz. 


MOND 


cyan long dash-dotted line 


Dark matter 


orange solid horizontal 



132 



Chapter 4. GALAXY ROTATION CURVES 




r [kpc] r [kpc] 

HSB (r) UGC 6399 HSB (s) UGC 6973 



Figure 4.6 continued: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc. 



T(r) 


Best fit 


with horizontal lines 


Mass to light profile 


with lines 


Newtonian core 


brown dot-dotted horizontal 


Newton's theory 


brown dot-dotted line 


STVG 


black solid horizontal 


STVG 


black solid line 


MSTG 


blue short dash-dotted horiz. 


MSTG 


blue short dash-dotted line 


MOND 


cyan long dash-dotted horiz. 


MOND 


cyan long dash-dotted line 


Dark matter 


orange solid horizontal 



4.2. URSA MAJOR FILAMENT OF GALAXIES 



133 




5 10 
r [kpc] 
LSB (a) NGC 3917 



15 



5 10 
r [kpc] 
LSB (b) NGC 4010 




10 
r [kpc] 
LSB (c) NGC 4183 




5 10 
r [kpc] 
LSB (d) UGC 6446 



15 




r [kpc] 
LSB (e) UGC 6667 



6 
5 
4 

3 F- 
2 
1 





r [kpc] 
LSB (f) UGC 6818 



Figure 4.6 continued: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc. 
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Figure 4.6 continued: UMa — Stellar mass-to-light ratios. 

The stellar mass-to-light ratio, T(r), vs. orbital distance, r in kpc for 19 HSB and 10 LSB 
galaxies. The stellar mass-to-light ratio, T, required to fit the galaxy rotation curve at 
each data point without dark matter is plotted for Moffat's STVG and MSTG theories, 
Milgrom's MOND theory, and Newton's theory. Within each theory, the best-fit value of 
the stellar mass-to-light ratio, T, is shown with a horizontal line, including the best-fit 
core-modified dark matter; and the best-fit Newtonian core model - the extent of the core 
shown with a vertical line. 
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to-light ratio changes within a galaxy, T = T(r), may be explored by calculating T(r) 
independently at every point in the galaxy rotation curve, for each gravity theory. 

Figure 4.6 plots the stellar mass-to-light ratio, T(r) vs.r, showing that the variation 
in T(r) in every galaxy exceeds the variation in the best-fit stellar mass-to-light ratio, T, 
across the sample of galaxies. 

4.3 Halos of phantom dark matter 

Whether galaxy rotation curves are best described by cold non-baryonic dark matter 
(CDM) halos, as in §4.1.1, or Milgrom's modified Newtonian dynamics, as in §3.2, or 
Moffat's modified gravity, as in §4.1.3, or other sensible theories, there are certain regular- 
ities in galactic structure which may have theoretical underpinning. Each of the galaxies 
studied in Brownstein and Moffat (2006a) and the Ursa Major filament of galaxies, in §4.2, 
exhibit a core region where Newton's theory provides acceptable fits and the galaxy may be 
modelled by a Newtonian core model, as in §4.3.1. The observation that mass follows light 
and the appearance of orphan features beyond the Newtonian core, as in §4.3.2, confirm 
that the baryons are dynamically important. This result is natural in MOND and MOG 
theories which are sourced by baryons alone, and also supports the alternative model of 
core-modified dark matter, as described in §4.3.3, which fits the galaxy rotation curves of 
§4.2.3, including all of the dwarfs, with physically reasonable stellar mass-to-light ratios of 
T « 1. This is difficult to achieve using the cuspy NFW profile for some LSB and dwarf 
galaxies which prefer T ~ 0. Comparison of the dynamic mass distribution inferred from 
galaxy rotation curves and the visible baryon distribution derived from each gravity theory 
enable a fundamental explanation to the Tully-Fisher relation, as in §4.3.4. 

4.3.1 Newtonian cores 

The radius of the Newtonian core is easily measured by plotting the stellar mass-to-light 
ratio, as in Figure 4.6, which shows the variation in the Newtonian T(r) in brown dot- 
dotted lines. The profile is flat T ~ 1 in the Newtonian core, and then rises rapidly 
outside the Newtonian core, as shown on the figure for each galaxy. Thus the best-fitting 
stellar mass-to-light ratio for Newton's theory without dark matter may be computed by 
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weighting the region inside the core radius, for each galaxy, as shown by the horizontal 
brown-dot-dotted lines in the figure. As a result, the best-fit Newtonian core model predicts 
values of T larger than the other theories because there is less gravity due to the visible 
baryons, without dark matter. 

4.3.2 Orphan features 

Kent (1986) presented a sample of 37 Sb and Sc galaxies with photometric data, discovering 
that the component mass distributions admit decompositions into baryon and dark matter 
components, but could not simultaneously constrain the dark matter distribution and 
the stellar mass-to-light ratio for the baryons. At one extreme, the stellar mass-to-light 
ratio was set to the maximum value permitted by the rotation curves, and a modest halo 
component produced good fits, but most galaxies were also well fit by models at the other 
extreme, with much more massive dark matter halos and correspondingly reduced stellar 
mass-to-light ratios. 

Kent (1987) presented a sample of 16 spiral galaxies with photometric data and ex- 
tended HI gas, and provided least-squares fits to the rotation curves, finding that a halo 
component is needed in each galaxy, but is tightly coupled to the stellar mass-to-light 
ratio. Although Kent (1986, 1987) assumed the constant density core dark matter distri- 
bution of Equation (2.3), the uncertainty in the stellar mass-to-light ratio is a result of the 
uncertainty in the dark matter distribution. 

However, this fine-tuning problem, which is known as the disk-halo conspiracy, is re- 
solved by a correlation between the shape of the rotation curve and the shape of the 
baryonic luminosity measurements, first observed by Burstein et al. (1982), which suggests 
the presence of some features in the rotation curves at the transition from the baryon 
dominated core to the dominant dark matter halo. Salucci and Frenk (1989) showed that 
the fractional amount of mass from the luminous disk is an increasing function of the lu- 
minosity, and argued that the shape of the rotation curve near the edge of the optical disk 
should vary systematically with luminosity, leading to distinct features in galaxy rotation 
curves. 

Hoekstra et al. (2001) applied a mass model in which the dark matter surface density 
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is a scaled version of the observed HI surface density to a sample of 24 spiral galaxies, 
obtaining good fits for most galaxies, but not for those galaxies which show a rapid decline 
of the HI surface density in the outermost regions. 

Noordermeer (2006) provided a systematic study of HI rotation curves in spiral galaxies, 
finding that galaxy rotation curves have distinct features that may be traced back to the 
luminous components in the form of bumps and wiggles, and that the declines in the 
rotation curves at intermediate to large radii are rarely featureless. 

4.3.3 Core-modified dark matter 

It is important to notice that the E-map components of Equations (4.31) and (4.32) com- 
bine to produce maps with features that can be traced back to the luminous stellar disk 
component for r < r s and to the gaseous disk component for r « r 5 , whereas the dark 
matter halo dominates for r > r 5 , as shown by the surface mass density maps, plotted 
in Figure 4.2. The core-modified dark matter surface mass density distribution including 
the visible baryons is remarkably flattened in the galaxy core, as compared to the best-fit 
Newtonian core, in the absence of dark matter, as in §4.3.1. Moreover, for every galaxy in 
the sample, the central surface mass density, 

S = S(0), (4.54) 

is determined by the baryonic component alone, where the dark matter halo is rarified as a 
result of the core-modified model of Equation (4.2). This is precisely the reverse situation 
for the singular halo models of §4.1.1: Navarro- Frenk- White profile which dominate the 
Newtonian dynamic mass throughout the galaxy, including the core leading to artificially 
small stellar mass-to-light ratios, T C 1 as shown by the best-fit NFW and core-modified 
parameters, listed in Table 4.2. 

Therefore the reasoning that dominant dark matter erases the orphan features visible 
in galaxy rotation curves and the derived S-maps applies only to the NFW profile, and not 
to the core-modified profile. The excellent fits to the galaxy rotation curves, in Figure 4.1, 
confirm that the baryonic components are dynamically important, but the Newtonian force 
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law of Equation (4.7) without dark matter fits none of the galaxy rotation curves outside 
the Newtonian core. A comparison of Table 4.2, show a statistically significant reduction 
of the x 2 l v test in ~ 90% of the galaxies results from using the core-modified profile of 
Equation (4.2) instead of the NFW profile of Equation (4.1) because the orphan features 
in the dynamic data are correctly repatriated with the baryonic surface density maps. 

4.3.4 The Tully-Fisher relation 

The observational Tully and Fisher (1977) relation is an empirical relation between the 
measured total luminosity of a galaxy in a particular band (proportional to the stellar 
mass) and the amplitude of the rotation curve (the maximum velocity) of the form: 

L oc ^ ax where a « 3 - 4, (4.55) 

where the total luminosity 

L = 4^$d 2 , (4.56) 

may be inferred by measuring the isotropic flux, $, and knowing the distance, rf, between 
earth and the galaxy. Since all of the galaxies in the sample of §4.2 are at a common 
distance from the Milky Way, and because of the improvements identified in §4.2.1 from 
using the available near-infrared X-band, the large astronomical uncertainties are miti- 
gated, leaving an ideal laboratory to study the relationship between the luminosity of a 
galaxy, and theoretical predictions from §4.1, for each gravity theory, independent of the 
galactic mass distribution. 

Tully and Pierce (2000) showed that, although the exponent in Equation (4.55) depends 
on the wavelength of the measured luminosity, and increases systematically from B to K 
bands, there appears to be convergence in the near-infrared where 

a = 3.4 ±0.1. (4.57) 

McGaugh et al (2000); McGaugh (2005) studied a large sample of galaxies, with stellar 
masses ranging over five decades, and observed a change in slope in the Tully-Fisher relation 
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which disappears when using the total baryonic mass including both stellar and gaseous 
components, instead of just the luminous stellar mass, and concluded that the Tully-Fisher 
relation is fundamentally a relation between the total baryonic mass and the rotational 
velocity. Using a combination of X-band photometry and high resolution rotation curves, 
Noordermeer and Verheijen (2007) discovered a second change in slope at the high mass 
end of the Tully-Fisher relation which disappeared when using, in combination, the total 
baryonic mass and the asymptotic, outermost velocity instead of the velocity amplitude. 

Considering the asymptotic behaviour of the galaxy rotation curves of Figure 4.1, most 
high-resolution galaxy rotation curves are either slowly rising or slowly declining at large 
radii. Verheijen (2001) considered an alternate definition of the "flat rotation velocity", 
categorizing galaxies according to three kinds of behaviour depending on the shape of the 
rotation curve. 

Avila-Reese et al (2008) explored the variation in the Tully-Fisher relation using a 
large sample of 76 high and low surface brightness galaxies, and obtained a = 3.40 for the 
ordinary Tully-Fisher relation (where the stellar luminosity is taken proportional to the 
stellar mass), confirming Equation (4.57). However, the value of the exponent in Equation 
(4.55) may be as shallow as a = 3.00 for the baryonic Tully-Fisher relation, a = 2.77 for 
the actual £>-band, and a = 3.67 for the actual K-band, based on their sample. 
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In order to calculate the total X-band luminosity, apparent X-band magnitudes from 
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the 2MASS survey were used. Given an apparent K-band magnitude it is possible to 
calculate the i^-band luminosity as 

log 10 (L*:) = L364 - 2 -K T + log 10 (l + z) + 2 log 10 d, (4.58) 

where L K is the i^-band luminosity in units of 1O 1O L , K T is the i^-band apparent magni- 
tude and z is the redshift of the galaxy (determined from the NASA/IPAC Extragalactic 
Database), and d = 18.6 Mpc is the distance to the galaxy in the Ursa Major filament. 
The log 10 (l + z) term is a first order X-correction. 

The empirical K-band Tully-Fisher relation is plotted in Figure 4.7, with the ordinary 
relation in blue, including the best-fit power-law, 

L K oc *w 4 - 1±0 - 4 , (4.59) 

of the form of Equation (4.55). To consider the effect of identifying the velocity in the 
Tully-Fisher relation with the asymptotic velocity, v ont) instead of the maximum velocity, 
v max , the asymptotic K-band Tully-Fisher relation is plotted in red, including the best-fit 
power-law with results listed in Table 4.7. 

The empirical Tully-Fisher relation involves the total luminosity in a particular band, 
such as the K-band, which is proportional to the stellar disk mass through the stellar 
mass-to-light ratio, 

L = Y' (460) 

log(M) = a log(» + b - log (T) . (4.61) 

Thus, the effect of T 7^ 1 is to shift the log(M)-intercept; but does not affect the 
slope. Theoretical predictions may be may quantified by either computing the appropriate 
T values which depend on the particular band of the luminosity measurements, or by 
considering the respective curve fits to the actual Tully-Fisher relation: 



log(M) = a\og(v) + b. 



(4.62) 
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Table 4.7: Tully-Fisher relation 







V 


max 


^out 


Relation Type 


a 


b 


a 


b 


(i) 




(2) 


(3) 


(4) 


(5) 


Empirical 


K-band 


4.1 ±0.4 


-8.5 ±0.8 


4.6 ± 0.4 


-9.5 ±0.9 




STVG 


3.2 ±0.3 


-6.9 ±0.7 


3.5 ±0.4 


-7.4 ±0.9 


Stellar mass 


MSTG 


2.9 ±0.3 


-6.0 ±0.7 


3.1 ±0.4 


-6.4 ±0.9 


MOND 


4.1 ±0.4 


-8.7 ±0.8 


4.4 ±0.5 


-9.3 ± 1.0 




Dark Matter 


3.0 ±0.4 


-6.2 ±1.0 


3.1 ±0.5 


-6.5 ± 1.2 




STVG 


2.6 ±0.2 


-5.4 ±0.5 


2.8 ±0.3 


-5.7 ±0.7 


Baryonic mass 


MSTG 


2.5 ±0.2 


-5.0 ±0.5 


2.6 ±0.3 


-5.3 ±0.7 


MOND 


3.0 ±0.3 


-6.3 ±0.5 


3.3 ±0.3 


-6.7 ±0.7 




Dark Matter 


2.5 ±0.3 


-5.1 ±0.7 


2.7 ±0.4 


-5.4 ±0.9 


Total mass 


Dark Matter 


2.9 ±0.2 


-5.5 ±0.5 


3.1 ±0.3 


-5.9 ±0.7 



Notes. — The empirical and theoretical Tully-Fisher relation: Column (1) lists the relation type, 
where the empirical Tully-Fisher relation is the measured K-band luminosity, vs. velocity, 
plotted in Figure 4.7; and the theoretical Tully-Fisher relation identifies the total luminosity with 
either the stellar disk mass, in the ordinary case, or the baryonic mass including the stellar disk 
and HI (plus He) gas mass, each per gravity theory: Moffat's STVG and MSTG theories and 
Milgrom's MOND theory and the core-modified dark matter theory. The total mass relation 
identifies the combined masses of the stellar disk, HI (and He) gaseous disks and dark matter 
halo as the source of galactic dynamics. Columns (2) and (3) list the power-law index and 
proportionality constant of Equation (4.62), respectively, for the ordinary Tully-Fisher relation 
with v = % ax . Columns (4) and (5) list the power-law index and proportionality constant of 
Equation (4.62), respectively, for the asymptotic Tully-Fisher relation with v = v out . 

As regards dark matter, Figure 4.8 plots the total mass including the luminous baryonic 
components and the dark matter halo, M tota i = Mb ary0 n + ^fhaio? according to the best-fit 
core-modified dark matter theory, including the best-fit power-law, with results listed in 
the bottom row of Table 4.7, finding a minimum of scatter in the best-fit power-law, 

{ ? , 2.9±0.2 
^max -) , . \ 

4.63 
7 , 3.1±0.3 
^out 



The result is significant because it provides an empirical relation to determine the total 
mass of a galaxy, from a few simple dynamical velocity measurements. 

In the ordinary case, without dark matter, M is identified with Mdi S k- However, the 
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Figure 4.8: UMa — To- 
tal mass with dark matter vs. velocity relation 

The total mass, M tota i = M baryon + 
^haio, i n 1O 1O M , vs. velocity, v in km 
s- 1 for 19 HSB and 10 LSB galaxies, 
with M ha i of Equation (4.2). The 
velocity is identified according to 
Equation (4.65) as either v = ^ max , 
with blue circles, or as u = v out , 
with red crosses, with best-fits shown 
with solid blue and red dashed lines, 
respectively, with parameters listed 
in the bottom row of Table 4.7. 



7© 10 1 



10" 



50 



: : 



no . 1 1 
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theoretical underpinning of the Tully-Fisher relation suggests that the empirical relation 
is an approximation to the baryonic Tully-Fisher relation which identifies M with Mb ary0 n, 
instead, because in the absence of non-baryonic dark matter, galactic dynamics are sourced 
by the total baryonic mass, and not just the luminous stellar disk. Figure 4.9 plots four 
distinct Tully-Fisher relations, where the luminosity is taken proportional to the stellar 
mass in the left panels and the baryonic mass in the right panels, and where the velocity is 
taken as the velocity amplitude (maximum) in the top panels, and the asymptotic velocity 
at the position r out in the bottom panels: 



M 



Mdisk Ordinary Tully-Fisher relation, 
^baryon Baryonic Tully-Fisher relation, 

and for each, two depending on whether 



(4.64) 



v = 



^out 



Ordinary Tully-Fisher relation, 
Asymptotic Tully-Fisher relation. 



(4.65) 



The scatter in the theoretical Tully-Fisher relations is minimized in the case of the 
baryonic Tully-Fisher relation across all gravity theories, except dark matter, implying 
that the empirical Tully-Fisher relation - which involves only the luminous disk - is an 
approximate law. 
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Figure 4.9: UMa — Theoretical Tully-Fisher relations. 

The mass, M, in 1O 1O M , vs. velocity, v in km s" 1 for 19 HSB and 10 LSB galax- 
ies. The mass is identified according to Table 4.6 as either M = M^sk m 
the left panels, or as M = Mb ary0 n, in the right panels, for Moffat's STVG 
and MSTG theories and Milgrom's MOND theory and the core-modified dark 
matter theory. The velocity is identified according to Equation (4.65) as ei- 
ther v = u max , in the top panels, or as u = v out , in the bottom panels, with 
best-fits shown per gravity theory, with parameters tabulated in Table 4.7. 



Tully-Fisher relations 



Theoretical predictions 


with error bars 




Best-fit power-law 


with lines 


STVG 


black solid error bars 




STVG 


black solid line 


MSTG 


blue short dash-dotted error 


bars 


MSTG 


blue short dash-dotted line 


MOND 


cyan long dash-dotted error 


bars 


MOND 


cyan long dash-dotted lines 


Dark matter 


orange solid error bars 




Dark matter 


orange solid line 



Chapter 5 



Clusters of galaxies 

"A table, a chair, a bowl of fruit and a violin; what else does a man need to be 
happy?" 

Albert Einstein 



Smail et al. (1995) argued that as clusters of galaxies are the largest bound structures 
known in the universe, their mass-to-light ratios and baryonic fractions should approach 
that for the cosmos as a whole. Whereas X-ray luminosity measurements typically give 
temperature distributions an order of magnitude larger than observed from fits to observed 
isothermal gas spheres, there is a remarkable variation in the size and shapes of the X-ray 
distributions, to foil the search for a universal description of the phantom of dark matter, 
within the modified gravitation theory of Chapter 3. 

Brownstein and Moffat (2006b) applied the modified acceleration law of metric skew- 
tensor gravity (MSTG), as described in §3.3, obtained from the Yukawa skewon theory 
of §3.3.4 in which Einstein gravity is coupled to a Kalb-Ramond Proca field, as in §3.3.2. 
Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot 
intracluster medium (ICM) with King /3-models, the predicted X-ray surface brightnesses 
of the sample of 106 X-ray clusters were consistent without introducing a non-baryonic dark 
matter component. The sub-kiloparsec X-ray surface brightness distributions, as functions 
of radial distance, are well matched across the sample, including the correct shape, by 
the /?-model fits to the X-ray surface brightness distributions arising from the modified 
acceleration law. 

Prompted by the observed ring-like feature of the weak-lensing map of the galaxy 
cluster CL 0024+1654, Milgrom and Sanders (2008) argued that despite any underlying 
feature in the baryon distribution, the ring may be observed as the image of the MOND 
transition region. The possibility that these emergent features appear as phantom dark 



145 



146 



Chapter 5. CLUSTERS OF GALAXIES 



matter in the strong and weak lensing mass reconstruction of Jee et al (2007) indicates the 
degree to which MOG theories violate the strong equivalence principle in order to describe 
clusters of galaxies in the absence of dark matter. The same phenomenon applied to the 
Bullet Cluster 1E0657-558 produces the observed phantom dark matter in the strong and 
weak lensing mass reconstruction of Clowe et al. (2006a); Bradac et al. (2006); Clowe et al. 
(2007), in the form of spatially dislocated peaks. 

The physics of X-ray clusters is derived in §5.1. Dark matter distributions are computed 
and compared to actual gas mass measurements for each of the clusters of galaxies, with 
best-fit cluster models presented in §5.2. Direct evidence from the Bullet Cluster 1E0657- 
558 gravitational lensing experiment is presented in §5.3 which supports the necessity of 
dominant dark matter, or the modified gravity alternative. 

5.1 X-ray clusters 

The creation of X-ray mass profiles from astrophysical observations, as described in §5.1.1 
for clusters of galaxies, is subject to model dependent assumptions based on the isotropic 
isothermal model, as in §5.1.2. The road from measuring radial, X-ray temperature profiles, 
to surface mass density maps, as in §5.1.3, in some chosen gravity theory, takes its way 
through the dynamics of the isothermal sphere and dynamical mass computations, as 
described in §5.1.4, which result in best-fit cluster models, presented in §5.2. 

5.1.1 Astrophysical observations 

Clusters of galaxies have been known to require some form of energy density that makes 
its presence felt only by its gravitational effects since Zwicky (1933) analysed the velocity 
dispersion for the Coma cluster. The more than 1000 galaxies spherically distributed 
within the Coma cluster comprise a small fraction (10%) of the baryonic mass, the larger 
fraction consisting of a diffuse cloud of 100 million degree X-ray emitting plasma - the 
intracluster medium (ICM). The ICM itself comprises only a small fraction (10%) of the 
Newtonian dynamic mass as determined from X-ray luminosity measurements. 

Much closer to the Milky Way, the Virgo cluster forms the heart of the Local Super- 
cluster, and has a galaxy population as rich as Coma distributed in three groups. Messier 
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49 - an elliptical / lenticular galaxy - is the brightest member of the Virgo cluster and is 
the center of one of the subdominant groups. The ICM surrounding Messier 49 is a diffuse 
cloud of 10 million degree X-ray emitting plasma, but only accounts for a tiny fraction (1%) 
of the Newtonian dynamic mass as determined from X-ray luminosity measurements. The 
Fornax cluster is much smaller than Virgo, but at a similar distance from the Milky Way. 
The 15 million degree ICM surrounding the Fornax core - which is in the preliminary pre- 
heating stage of an imminent merger as determined from peculiar velocity measurements 
along a filament - comprises a similar fraction (3%) of the Newtonian dynamic mass. 

Abell 400 is an ongoing cluster-cluster merger, with multiple subclusters around a 
central main cluster containing the Dumbbell galaxy, which is the result of a galaxy- 
galaxy merger and is the topic of ongoing X-ray and radiowave analysis due to a pair of 
suspected supermassive black holes, bound and moving together. The 30 million degree 
X-ray emitting plasma accounts for 10% of the Newtonian dynamic mass. The Hydra- 
Centaurus supercluster contains two distinct X-ray clusters, each with 100 member galaxies 
near their respective centers, and have 50 and 40 million degree X-ray emitting ICM 
plasmas, respectively, which account for 10% of the Newtonian dynamic masses. At the 
center of the Great Attractor 65 Mpc distant, the Norma cluster is half the size of the 
Coma cluster, but larger than Centaurus, Hydra-A, Fornax and Messier 49 combined. It is 
in the process of swallowing a galaxy which shows a comet-like tail nearly twice as long as 
the galaxy itself, consisting of a 70 million degree X-ray plasma, and accounting for 10% 
of the Newtonian dynamic mass. 

Perseus is the brightest X-ray cluster in the sky and is nearly the size of the Coma 
cluster, but is not as rich in galaxies. The 80 million degree X-ray plasma accounts for 
20% of the Newtonian dynamic mass. Chandra has measured concentric ripples in the X- 
ray surface mass density surrounding a strong source of gravitation inside an X-ray cavity 
- a candidate for a 10 8 M black hole. Abell 2255 is only slightly larger than Perseus, and 
slightly less than the size of Coma, but the X-ray peak is offset from the brightest cluster 
galaxy, which has a large peculiar velocity, (1200 km s _1 ), indicating an ongoing cluster 
merger. The 80 million degree plasma accounts for 8% of the Newtonian dynamic mass. 
The giant Abell 2142 is one-and-a-half times larger than the Coma cluster, and is in the 
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later stages of a cluster-cluster merger showing bow-shock waves. The 110 million degree 
plasma accounts for 15% of the Newtonian dynamic mass. 

The Bullet Cluster 1E0657-558 is a merger between a giant main cluster with thousands 
of galaxies and a supersonic subcluster, aligned in the plane perpendicular to the line- 
of-sight. The 170 million degree main cluster ICM accounts for 10% of the Newtonian 
dynamic mass, and provides strong and weak gravitational lensing observations which 
show structure offset from the X-ray surface density map. 

The sample selection includes the Bullet Cluster 1E0657-558, Abell 2142, Coma, Abell 
2255, Perseus, Norma, Hydra-A, Centaurus, Abell 400, Fornax, and Messier 49, with 
cluster properties listed in Table 5.1 - ordered from the hottest X-ray emitting to the 
coolest of the clusters. The Newtonian dynamic masses and the ICM gas masses, for each 
cluster, are plotted in Figure 5.1, and the ratio of the Newtonian dynamic masses to the 
ICM gas masses are plotted as dynamic mass factors in Figure 5.2 - each as a function of 
radial position and compared to the theoretical predictions of core-modified dark matter 
halos, as in §5.2.1, Milgrom's MOND as in §5.2.2, and Moffat's MOG as in §5.2.3. 

The study of the Bullet Cluster 1E0657-558, in §5.3, is a detailed analysis of the X- 
ray gas surface density map in relation to the Newtonian dynamic mass inferred from the 
strong and weak gravitational lensing map. The missing mass in MOG is explained by 
galactic surface mass density maps, presented in Figure 5.9. The missing mass in terms of 
dark matter is presented in Figure 5.10. 

5.1.2 Isotropic isothermal model 

Recent observations from the XMM-Newton satellite suggest that the intracluster medium 
(ICM) is very nearly isothermal inside the region defined by the X-ray emission with tem- 
peratures ranging from ~ 1-15 keV (or 10 7 - 2 x 10 8 K) for different clusters (Arnaud 
et a/., 2001). The combination of the observed density profile, n e (r), and the temperature 
profile, T(r), obtained from X-ray observations of the galaxy cluster leads to a pressure 
profile, P(r), which directly leads to a mass profile, M(r), by assuming the gas is in nearly 
hydrostatic equilibrium with the gravitational potential of the galaxy cluster. Within a few 
core radii, the distribution of gas within a galaxy cluster may be fit by a King a /?-model". 
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The observed surface brightness of the X-ray cluster can be fit to a radial distribution 
profile (Chandrasekhar, 1960; King, 1966): 

-3)9+1/2 

(5.1) 

resulting in best-fit parameters, j3 and r c . A deprojection of the /3-model of Equation (5.1) 
assuming a nearly isothermal gas sphere then results in a physical gas density distribu- 
tion (Cavaliere and Fusco-Femiano, 1976): 

-3)9/2 

(5.2) 

where p(r) is the ICM mass density profile, and p denotes the central density. The mass 
profile associated with this density is given by 

Mir) = 4tt [ r p(r')r' 2 dr', (5.3) 
Jo 

where M(r) is the total mass contained within a sphere of radius r. Galaxy clusters 
are observed to have luminous distributions with finite spatial extent. This allows an 
approximate determination of the total mass of the galaxy cluster by first solving Equation 
(5.2) for the position, r out , at which the density, p(r out ), drops to « 10 -28 g/cm 3 , or 250 
times the mean cosmological density of baryons: 

1/2 

(5.4) 

Then, the total mass of the ICM gas may be taken as M gas w M(r out ): 

-3/3/2 

r f2 dr f . (5.5) 
Provided the number density, n, traces the actual mass, we may assume that n(r) oc 



I(r) = /o 



1 + 



p(r) = Po 



1 + 



' out 



Po 



2/3/3 



10~ 28 g/cm 



M gas = 4tt 



rout 



Po 
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p(r), which according to Reiprich (2001); Reiprich and Bohringer (2002) is explicitly 

p gas « 1.17n e m p , (5.6) 

and rewrite Equation (5.2) 



n e (r) = n 



1 + 



-3/3/2 



(5.7) 



For a spherical system in hydrostatic equilibrium, the structure equation can be derived 
from the collisionless Boltzmann equation 



where $(r) is the gravitational potential for a point source, a r and are mass- weighted 
velocity dispersions in the radial (r) and tangential (0, (j)) directions, respectively. For an 
isotropic system, 

a r = a e ^- (5.9) 
The pressure profile, P(r), can be related to these quantities by 

P(r) = a r 2 p(r). (5.10) 

Combining Equations (5.8), (5.9) and (5.10), the result for the isotropic sphere is 

dP(r) d<f>(r) 

-*r = ~ p{r) ^r- (5 - n) 

For a gas sphere with temperature profile, T(r), the velocity dispersion becomes 

o} = ^, (5.12) 
fi A m p 

where k is Boltzmann's constant, ha ~ 0.609 is the mean atomic weight and m p is the 
proton mass. We may now substitute Equations (5.10) and (5.12) into Equation (5.11) to 
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obtain 



d (kT(r) \ d*(r) 
Tr {— P(r))=-P(& 



(5.13) 



jLiAm p ' 'J x ' dr 

Performing the differentiation on the left hand side of Equation (5.11), we may solve for 
the gravitational acceleration: 



a(r) = 



d$(r) 
dr 
kT(r) 
ji A m p r 



dln(p(r)) d\n{T(r)) 
dhi{r) rfln(r) 



(5.14) 



For the isothermal isotropic gas sphere, the temperature derivative on the right-hand side of 
Equation (5.14) vanishes and the remaining derivative can be evaluated using the /3-model 
of Equation (5.2): 

<T) = -WV(*\ (5 . 16) 

fi A m p \r 2 + r 2 J 

5.1.3 Surface mass density map 

To make contact with the experimental data, we must calculate the surface mass density 
by integrating p(r) of Equation (5.2) along the line-of-sight: 



/Zout 
p(x,y,z)dz, 
~Znnt 



(5.16) 



where 



^out 



^out -x 2 -y 2 . 



Substituting Equation (5.2) into Equation (5.16), we obtain 



:(x,y) = PoJ_ 



Zout 
Zout 



1 + 



x 2 + y 2 + z 2 



-3/3/2 



dz. 



(5.17) 



(5.18) 



This integral becomes tractable by making a substitution of variables: 



2 x 2 + y 2 

U = 1 + = 



(5.19) 
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so that 



E(x,|/) 



po 



Po 

U 3P 



^out 
2>out 
Zo ut 



u 2 + 



21 "3/3/2 



^out 



1 + 



Z 



dz 

-3(3/2 



dz 



out 1 



I 3 

2' 2 



/3 



^out \ ' 

2/r r / 



(5.20) 



where we have made use of the hypergeometric function, F([a, 6], [c], z). Substituting 
Equation (5.19) into Equation (5.20) gives 



x 2 + y 2 " 



-3(3/2 



1 3 

2' 2 



^out 



x 2 + y 2 + r 



We next define 



(5.21) 



E = E(0, 0) = 2p z out F 



1 3 

2' 2 



/ ^OUtV 



(5.22) 



which we substitute into Equation (5.21), yielding 



E(x,y) = E 1 



x 2 + y 2 ^ 



-3(3/2 



' x 2 -\-y 2 -\-r 2 



Z 2 
^out 

^2 



(5.23) 



In the limit z out 3> r c , the Hypergeometric functions simplify to T functions, and Equations 
(5.22) and (5.23) result in the simple, approximate solutions: 



r (1/3) 



(5.24) 



and 



E(x,y) = E 1 + 



a; 2 + y2 X-(3/ 3 -l)/2 



(5.25) 



which we may, in principle, fit to the E-map data to determine the King /3-model param- 
eters, (3, r c and po- 
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5.1.4 Dynamical mass 

The Newtonian dynamical mass can be obtained as a function of radial position by equat- 
ing the gravitational acceleration of Equation (5.14) - derived in §5.1.2 for the isotropic 
isothermal model - with Newton's acceleration law: 



kT(r) 
fi A m p r 



dln(p(r)) dln(T(r)) 



dln(r) 



rfln(r) 



G N M N (r) 



with the solution, 



M N (r) 



r kT 



Gn ^Atn P [ rfln(r) 



dln(p(r)) dln(T(r)) 



d\n(r) 



and the isothermal (5- model result of Equation (5.14) can be rewritten as 



M N (r) 



3f3kT ( r 3 
li A m p G N \r 2 +r 



(5.26) 



(5.27) 



(5.28) 



5.2 Best-fit cluster models 

The study of X-ray clusters, according to §5.1, provides valuable information on their mass 
profiles and insight into their formation and evolution. It is no longer a matter of fitting 
the total masses of these systems, but a powerful means to model the spatial distribution of 
each component, as in the case of the Bullet Cluster 1E0657-558 presented in §5.3, which 
may further constrain cosmological models. Although X-ray luminosity measurements 
typically give temperature distributions an order of magnitude larger than observed from 
fits to observed isothermal gas spheres, this does not guarantee that the missing mass has 
the form predicted by A-CDM cosmological models, particularly because of the remarkable 
variation in the shapes and scales of the X-ray distributions. 

The core-modified dark matter halo, described in §5.2.1, provides an alternative which 
does not suffer from the cusp problem of the singular NFW fitting formula and does fit high 
and low surface brightness and dwarf galaxies to low x 2 . Alternatively, Milgrom's MOND 
as in §5.2.2 and Moffat's MOG as in §5.2.3 are to be explored, and the means by which 
each of the candidate theories provides best-fit cluster models are compared in §5.2.4. The 
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Table 5.2: Best-fit cluster model parameters 



Cluster 




Po 


r s 


X> 


Cluster 


a X 2 /" 




(10 6 M o /kpc 3 ) 


(kpc) 






(10- 8 cm s- 2 ) 


(1) 




(2) 


(3) 


(8) 


(l) 


(4) (8) 


Bullet 


2.73 ±0.47 


328.1 ±26.7 


4.5 


Bullet 


7 26 ± 42 653 9 


Abell 2142 


8.34 ±1.16 


141.8 ±7.7 


3.7 


Abell 2142 


3.14 ±0.22 1942.4 


Coma 


1.91 ±0.19 


262.1 ± 10.6 


3.9 


Coma 


6.05 ±0.25 162.6 


Abell 2255 


0.67 ±0.08 


403.0 ±21.0 


1.2 


Abell 2255 


5.04 ±0.22 311.8 


Perseus 


28.83 ±4.97 


58.7 ±3.9 


205.7 


Perseus 


1.63 ±0.15 2036.6 


Norma 


1.58 ±0.16 


222.7 ±8.3 


0.8 


Norma 


3.63 ±0.17 32.6 


Hydra-A 


22.79 ± 2.72 


59.3 ±3.0 


13.1 


Hydra-A 


2.79 ±0.16 4595.6 


Centaurus 


37.96 ± 2.48 


36.3 ±0.9 


104.3 


Centaurus 


3.99 ±0.15 626.9 


Abell 400 


2.09 ±0.18 


119.9 ±4.0 


3.7 


Abell 400 


2.51 ±0.1 575.0 


Fornax 


1.39 ±0.05 


117.4 ±1.8 


0.2 


Fornax 


17.46 ±0.18 1.7 


Messier 49 


36.75 ± 0.45 


21.3 ±0.1 


68.8 


Messier 49 


80.29 ± 0.27 229.0 


(a) 


Core- modified dark matter 




(b) Milgrom's MOND 


Cluster 


Mo 


ro 


X> 


Cluster 


£oo X> 






(1O 14 M ) 


(kpc) 






(G N ) 


(1) 




(5) 


(6) 


(8) 


(1) 


(7) (8) 


Bullet 


56.7 ±8.2 


116.8 ±8.5 


138.9 


Bullet 


8.7 ±0.5 56.7 


Abell 2142 


30.0 ±3.0 


56.8 ±2.9 


45.0 


Abell 2142 


8.7 ±0.5 375.0 


Coma 


30.7 ±2.6 


88.2 ±3.8 


14.3 


Coma 


12.3 ±0.5 56.6 


Abell 2255 


43.8 ±3.1 


157.4 ±5.6 


85.2 


Abell 2255 


12.4 ±0.5 42.1 


Perseus 


10.7 ± 1.0 


23.5 ± 1.1 


16.7 


Perseus 


7.5 ± 0.5 433.0 


Norma 


30.1 ± 1.9 


97.6 ±3.1 


2.9 


Norma 


11.8 ±0.5 12.5 


Hydra-A 


9.5 ±0.7 


23.9 ±0.9 


87.9 


Hydra-A 


10.8 ±0.5 740.5 


Centaurus 


10.0 ±0.3 


14.2 ±0.2 


3.3 


Centaurus 


18.4 ±0.5 107.5 


Abell 400 


6.0 ±0.3 


44.7 ± 1.3 


7.8 


Abell 400 


13.4 ±0.5 174.2 


Fornax 


13.7 ±0.2 


67.4 ±0.5 


1.3 


Fornax 


46.1 ±0.5 0.4 


Messier 49 


10.3 ±0.0 


10.8 ±0.0 


8.7 


Messier 49 


149.7 ±0.5 18.4 




(c) Moffat's MSTG 




(d) Moffat's STVG 



Notes. — Best-fitting parameters of the X-ray cluster sample for (a) Core-modified dark matter, 
(b) Milgrom's MOND, (c) Moffat's MSTG and (d) Moffat's STVG: Column (1) is the cluster 
name. Columns (2) and (3) list the best-fit parameters for the core-modified universal fitting 
formula of Equation (2.12). Column (4) list the best-fit MOND acceleration parameter. Columns 
(5) and (6) list the best-fit MSTG mass and range parameters, respectively. Column (7) list the 
best-fit STVG asymptotic gravitationally coupling. Column (8) list the reduced-x 2 statistic of 
Equation (4.20) per gravity theory. 



156 



Chapter 5. CLUSTERS OF GALAXIES 



common and unusual features of the fits, across the sample of clusters of galaxies, are 
presented in §5.2.4: Observations - with statistics provided in Table 5.2. 

5.2.1 Core-modified dark matter halos 

The simplicity of the dark matter paradigm allows predictions which justifiably ignore 
the X-ray surface mass. However, it has been known since Bahcall (1977); Hoffman and 
Shaham (1985) that the halo density profile of virialized clusters of galaxies cannot be 
fitted by a single power law, 

p(r) oc r~ 7 , (5.29) 

but 7 seems to increase with r. The NFW fitting formula of Equation (2.7) bridges this 
behaviour of Equation (5.29) with a singular, cuspy core, 7 — ► 1 at small r, and 7 —> 3 at 
large r. 

This point has led to controversy over the cuspy shape of the singular NFW profile, 
which appears as a robust prediction of N-body simulations without baryons, but is not 
actually observed (with low x 2 ) in X-ray cluster data. Conversely, Ettori et al (2002a) 
found 2 X-ray clusters out of a sample of 22 clusters that could not be x 2 fitted to the NFW 
profile at all. Further Chandra studies of the cores of clusters, such as Ettori et al (2002b), 
indicated that the NFW profile is only reliable outside the cluster core, as is the case in 
certain low surface brightness and dwarf galaxies. Sand et al (2004) studied the dark 
matter distribution in the central region of 6 clusters of galaxies by combining constraints 
from gravitational lensing and the stellar velocity dispersion profile of the brightest central 
galaxy, confirming that the core behaviour is statistically inconsistent with a singular NFW 
profile, and that the inclusion of baryonic matter affects the dark matter distribution not 
accounted for in conventional CDM simulations. 

In a gravitational lensing study of two X-ray clusters of galaxies, Smail et al (1995) 
found that the Newtonian dynamic mass is distributed similarly to the visible baryonic 
mass with the same core radius, but while it is more concentrated at the center, it is less 
cuspy than CDM predictions. Tyson et al (1998); Shapiro and Iliev (2000) argued that 
the singular density profiles based on NFW fitting formula are in apparent conflict with 
the observed mass distributions inside dark matter dominated halos on two extremes of 
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the halo mass function - dwarf galaxies and clusters of galaxies; each of which is better 
described by a less cuspy or constant density core. 

The core-modified fitting formula of Equation (2.12) bridges the behaviour of Equation 
(5.29) with a constant-density core, 7 — > at small r but 7 — > 3 at large r, which fits the 
high and low surface brightness galaxies in the Ursa Major sample of §4.2 including all of 
the dwarf galaxies. The dark matter power-law profile, plotted in Figure 4.5, confirms that 
the variation in the exponent of Equation (5.29) agrees with CDM predictions, provided 
the visible baryonic components are not neglected. It is therefore important to test the 
core-modified dark matter fitting formula of Equation (2.12) at the scale of X-ray clusters. 

Arieli and Rephaeli (2003) compared the best-fits to a sample of 24 X-ray clusters of 
galaxies, and concluded that a core-modified dark matter profile of the form of Equation 
(2.12) is statistically more consistent with ROSAT observational results than either the 
NFW profile of Equation (2.7) or a family of simple polytropic fitting formulae. 

Whereas, attempts to fit cluster mass distributions to NFW profiles lead to large uncer- 
tainties due to a parameter degeneracy between the central density parameter, p 0j an d the 
scale radius, r S) which prevented the computation of the best-fit p and r s from converging, 
regardless of the x 2 . Without numerical convergence, the NFW results either over-predict 
the density at the core or under-predict the total mass. However, the core-modified fitting 
formula of Equation (2.12) provides excellent fits with low x 2 to the mass profiles, plotted 
in Figure 5.1, and a reasonable explanation of the variation in the dynamic mass factors, 
plotted in Figure 5.2, providing one solution to the missing mass problem, presented in 
§5.2.4. 

5.2.2 Milgrom's MOND without dark matter 

Brownstein and Moffat (2006b) predicted convergent MOND X-ray surface brightness pro- 
files which did not match any observed distributions of a sample Reiprich (2001); Reiprich 
and Bohringer (2002) of 106 X-ray clusters. Without treating the MOND acceleration, a , 
as a free parameter as opposed to a universal constant, or considering improved but as yet 
undiscovered MOND interpolating functions, MOND cannot account for the observed X- 
ray luminosities without the addition of an unseen component to explain away the missing 
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mass. Sanders (2003, 2007) considered adding a neutrino halo, modelled as a nonluminous 
constant density rigid sphere, discussed in §5.3.5. 

Conversely, The and White (1988) were able to decrease the MOND discrepancy be- 
tween the X-ray observationally determined gas mass and the X-ray surface brightness of 
the Coma cluster by increasing the MOND acceleration by a factor of four greater than 
Equation (4.8). However, Aguirre et al (2001) presented evidence from the central 200 
kpc of three clusters which inflates the discrepancy in the MOND acceleration to a factor 
of ~ 10. More recently, Pointecouteau and Silk (2005) used X-ray data from the XMM- 
Newton satellite for eight clusters of varying temperature and masses to place constraints 
on the general use of MOND phenomenology. 

Furthermore, every galaxy rotation curve that produced a weak fitting MOND one- 
parameter best-fit by a variable stellar mass-to-light ratio, T, plotted in Figure 4.1 for the 
Ursa Major sample of §4.2, shows dramatic improvement and reduction in the reduced x 2 / l/ 
statistic using a two-parameter best-fit including a variable MOND acceleration parameter. 
The tabulation of a in Column (2) of Table 4.3 provides no statistical support that a is 
a universal constant due to gross uncertainties in the mean results of Equation (4.9). 

A varying choice of the MOND interpolating function, including those of Bekenstein 
(2004) and Famaey and Binney (2005), does not lead to significant improved behaviour 
since a(r) < a or x < 1 at all radii within clusters of galaxies. 

Therefore the alternatives for MOND are either add two additional parameters (or scal- 
ing relations) per system to include a dark matter component, or to determine if sensible 
fits are possible without dark matter using a one-parameter best-fit by a variable accelera- 
tion parameter. The absence of a universal acceleration parameter violates the notion that 
MOND is a fundamental theory, but the notion of a scale dependent acceleration parame- 
ter may be a dynamic, more natural effect of a covariant, but Lorentz- violating theory with 
a preferred frame, as in §3.2.3 and is not inconsistent with Bekenstein's TEVES action, as 
in §3.2.3. 
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5.2.3 Moffat's MOG with running couplings 

In the absence of non-baryonic dark matter, the modified gravity dynamical mass may 
be obtained as a function of radial position by substituting the MOG acceleration law of 
Equation (3.90) - with a varying gravitational coupling, G{r) - so that the result for the 
isothermal /?-model of Equation (5.28) becomes 



Brownstein and Moffat (2006b) predicted X-ray surface brightness profiles from X-ray 
luminosity observations consistent with the observed X-ray gas distributions of a sample 
of 106 X-ray clusters (Reiprich, 2001; Reiprich and Bohringer, 2002) using the modified 
acceleration law based upon metric skew-tensor gravity, as in §3.3. 

Metric skew-tensor gravity 

The MSTG dynamic mass is obtained by substituting G(r) of Equation (3.91) into Equa- 
tion (5.30) and may be written explicitly as a function of the Newtonian dynamic mass of 
Equation (5.28) and two parameters, M and r : 



which are MSTG mass and range parameters, respectively. However, it is not possible to 
fit any of the clusters of galaxies with the MSTG mass and range parameters of Equation 
(4.14), which were applied universally to high and low surface brightness galaxies including 
all of the dwarfs, in the Ursa Major sample of §4.2, with galaxy rotation curves plotted 
in Figure 4.1. Whereas every weak fitting MSTG one-parameter best-fit by a variable 
stellar mass-to-light ratio, T, shows dramatic improvement and reduction in the reduced 
X 2 /v statistic using a three-parameter best-fit including variable MSTG mass and range 
parameters, the tabulation of M and r in Column (2) and (3) of Table 4.4 provides no 
statistical support that the MSTG parameters are universal constants, but does provide 




(5.30) 



M MSTG (r) = M N (r) + M £(r) - ^M 2 f(r) 2 + 2M M N (r)^r) 



(5.31) 




(5.32) 
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very strong statistical support that the MSTG parameters are scale dependent. 

Brownstein and Moffat (2006b) used an empirically determined power-law scale relation 
to set the MSTG mass scale parameter, 

/ M \°- 39±0 - 10 
Mo = (60.4 ± 4.1) x 1O 14 M , (5-33) 

where M gas , given by Equation (5.5), is the mass of the ICM integrated to the distance at 
which the density drops to « 10 _28 g/cm 3 , or 250 times the mean cosmological density. 

In order to better determine the scale dependence of the parameters, it is reasonable 
to treat the MSTG mass and range parameters as variable and to perform two-parameter 
best-fits to the X-ray gas masses of the sample of 11 clusters of galaxies, using Equation 
(5.33) as initial value only. The mass profiles are plotted in Figure 5.1 according to the 
best-fit cluster model parameters tabulated in Panel (c) of Table 5.2, for MSTG. 

Scalar-tensor-vector gravity 

Moffat and Toth (2009) investigated a fundamental parameter-free solution to the running 
couplings using the modified acceleration law based upon scalar-tensor- vector gravity, as 
in §3.4. The STVG dynamic mass of Equation (3.164) may be written as a function of the 
Newtonian dynamic mass of Equation (5.28) and two functions a(r) and ji{r) which are 
derived from an action principle, with the equations of motion given by Equations (3.159) 
and (3.160), respectively, in terms of three constants of integration, L>, and G^. 

However, it is not possible to fit any of the clusters of galaxies with the values of 
Equation (4.19), which were applied universally to high and low surface brightness galaxies 
including all of the dwarfs, in the Ursa Major sample of §4.2, with galaxy rotation curves 
plotted in Figure 4.1. Whereas every weak fitting STVG one-parameter best-fit by a 
variable stellar mass-to-light ratio, T, show dramatic improvement and reduction in the 
reduced x 2 / l/ statistic using a four-parameter best-fit including variable parameters, the 
tabulation of D, E, and G^ in Columns (2), (3) and (4) of Table 4.5 provides no statistical 
support that the STVG integration constants are universal. 

For values of D sufficiently large and values of E sufficiently small, the STVG gravita- 
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tional coupling of Equation (3.143) simplifies to its asymptotic form, 



G(r) = G, 



(5.34) 



and is independent of r. Substituting this form of the gravitational coupling into Equation 
(5.30), we obtain the STVG dynamic mass for clusters of galaxies: 



Therefore, in order to determine the scale dependence of the STVG asymptotic coupling, it 
is reasonable to treat as variable and to perform one-parameter best-fits to the X-ray 
gas masses of the sample of 11 clusters of galaxies, plotted in Figure 5.1 according to the 
best-fit cluster models parameters tabulated in Panel (d) of Table 5.2, STVG. 

5.2.4 The missing mass problem 

The Newtonian dynamical mass of Equation (5.28) is a derived relation between the density 
profile for the X-ray gas component, according to the isotropic isothermal model of §5.1.2, 
and the measured isothermal temperature, T. The ICM gas mass is a spherical integral 
of the King /3-model of Equations (5.2) and (5.3). Table 5.1 includes the total ICM gas 
mass and total Newtonian dynamical mass within the position, r out , at which the density, 
p(^out), drops to « 10 -28 g/cm 3 , or 250 times the mean cosmological density of baryons. 
The total fraction of ICM gas mass is between 1% and 20% of the total Newtonian dynamic 
mass, and is typically 10%, as demonstrated in §5.1.1. Therefore, according to Newtonian 
dynamics, between 80% to 99% of the mass needed to explain the isothermal profiles is 
missing. 

However, whereas the solution that there is just enough dark matter to fill the total 
difference is consistent with the NFW fitting formula of Navarro et al (1996, 1997), the 
cusped profile does not correctly fit the shape of the dynamic mass profile. Arieli and 
Rephaeli (2003) suggested that there is a clear need to explore modifying the NFW pro- 
file, which has been adopted in hydrodynamic N-body simulations of the structure and 
evolution of A-CDM halos, or finding an alternative which provides a reasonable fit to the 




(5.35) 
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Figure 5.1: X-ray clusters of galaxies - Mass profiles. 



The radial mass profile, M(r) in M , vs. orbital distance, r in kpc for a sample of 
X-ray clusters. The dynamic data consist of the Newtonian dynamic mass of Equation 
(5.28), due to the measured isothermal temperature. The observed ICM gas masses 
are derived from Equations (5.2) and (5.3) using the best-fit King /3-model parameters 
listed in Table 5.1. The computed best-fitted results are plotted for Moffat's STVG and 
MSTG theories and Milgrom's MOND theory with variable parameters. Results are 
plotted for the best-fit core-modified dark matter theory including the X-ray gas mass 
component. The reduced-;^ 2 statistic is included in Table 5.2. The figure is continued. 

X-ray cluster masses. 

Similarly, it is not enough for any gravity theory to solve the missing mass problem, 
in the absence of dark matter, without providing a reasonable fit to the observed X-ray 
gas mass distribution for each cluster. The mass profiles of Figure 5.1 plot the Newtonian 
dynamical mass and observed ICM gas mass profiles, including the best-fits resulting from 
Moffat's STVG and MSTG theories, Milgrom's MOND theory with variable parameters. 
Results are plotted for the best-fit core-modified dark matter theory including the X-ray 
gas mass component. The reduced-;^ 2 statistic is included in Table 5.2, and reveals that the 
dark matter solutions of §5.2.1 and the modified gravity solutions of §5.2.3 are reasonable, 
although the MOND solution without dark matter of §5.2.2 is wrong - and a variable 
MOND acceleration parameter only allows a correct fit to the total cluster mass. This 
has prompted Sanders (2003, 2007) to consider the possibility of 2 eV neutrino halos as 
providing the missing 80% to 99% of cluster dark matter, but Angus et al (2008) showed 
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Figure 5.1 continued: X-ray clusters of galaxies - Mass profiles. 

The radial mass profile, M(r) in Mq, vs. orbital distance, r in kpc. 



M(r) 


Theory 


with lines 


Dark matter 


orange solid line 


Observation 


with error bars 


STVG 


black solid line 


Newtonian dynamic mass 


red crosses with error bars 


MSTG 


blue short dash-dotted line 


ICM gas mass 


green circles with error bars 


MOND 


cyan long dash-dotted line 



164 



Chapter 5. CLUSTERS OF GALAXIES 




Figure 5.1 continued: X-ray clusters of galaxies - Mass profiles. 

The radial mass profile, M(r) in Mq, vs. orbital distance, r in kpc for a sample of X-ray 
clusters. The dynamic data consist of the Newtonian dynamic mass of Equation (5.28), due 
to the measured isothermal temperature. The observed ICM gas masses are derived from 
Equations (5.2) and (5.3) using the best-fit King /3-model parameters listed in Table 5.1. 
The computed best-fitted results are plotted for Moffat's STVG and MSTG theories and 
Milgrom's MOND theory with variable parameters. Results are plotted for the best-fit 
core-modified dark matter theory including the X-ray gas mass component. The reduced- 
X 2 statistic is included in Table 5.2. 
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that MOND-neutrino-baryon models will not provide reasonable fits to the X-ray gas mass 
profile, particularly in the inner 100 to 150 kiloparsecs of the cluster. This neutrino halo 
hypothesis is explored in §5.3.5 as part of the analysis of the strong and weak lensing map 
of the Bullet Cluster 1E0657-558, presented in §5.3. 

Biviano and Salucci (2006) derived mass profiles of the different luminous and dark 
components of 59 X-ray clusters of galaxies and confirmed that the baryonic components 
are relevant to mass models of clusters of galaxies both near the center because of the 
substantial contribution from the central dominant galaxy and in the outer regions, because 
of the increasing mass fraction of the ICM gas - and the corresponding decreasing dynamic 
mass factor. Therefore the missing mass problem is most serious in the core of galaxy 
clusters, in complete opposition to the situation in the galaxy rotation curves of Chapter 4 
where the dynamical mass factors of Figure 4.4 show a maximum at the outermost observed 
radial position, for each gravity theory. 

The dynamical mass factors, plotted in Figure 5.2, show the ratio of the Newtonian 
dynamical mass to the observed ICM gas mass, including the best-fits resulting from 
Moffat's STVG and MSTG theories, Milgrom's MOND theory with variable a , and the 
best-fit core-modified dark matter theory including the X-ray gas mass component, where 
the results of Table 5.2 were used, respectively. 

For each cluster, the substitute of missing mass in MOND is the wrong shape, with 
only the correct total mass predicted due to a variable, best-fit MOND acceleration. For 
r < r out> the dynamic mass factor predicted by MOND is much smaller than observed 
leading to too great a predicted gas mass in these regions. For some of the clusters such 
as the Bullet Cluster 1E0657-558, Abell 2255 and Fornax, this trend is suddenly reversed 
for r < 100 kpc, where MOND predicts a dynamic mass factor which diverges strongly 
(as does the cuspy NFW profile not shown), but is not actually observed in the data even 
though the coolest of the clusters such as Messier 49 show a dynamic mass factor as large 
as T -> 400 as r -> 0. 

Unlike the NFW fitting formula of Navarro et al. (1996, 1997), the core-modified dark 
matter halos provide the means to fit X-ray masses with constant density cores. This 
solution provides missing mass in line with the observations plotted in Figure 5.2 at all 
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Figure 5.2: X-ray clusters of galaxies - Dynamic mass factors. 

The dynamic mass factors, T(r), vs. orbital distance, r in kpc for a sample of X-ray 
clusters. The dynamic data consist of the ratio of the Newtonian dynamic mass of Equa- 
tion (5.28), due to the measured isothermal temperature, to the integrated X-ray gas 
mass, derived from Equations (5.2) and (5.3) using the best-fit King /3-model parameters 
listed in Table 5.1. The computed best-fitted results are plotted for Moffat's STVG and 
MSTG theories and Milgrom's MOND theory with variable parameters. Results are 
plotted for the best-fit core-modified dark matter theory including the X-ray gas mass 
component. The reduced-^ 2 statistic is included in Table 5.2. The figure is continued. 



radial positions. 

Moffat's MOG theories provide a remarkable picture of the missing mass problem, 
even though the galactic mass components have been neglected, which are dynamically 
important in MOG due to the absence of dominant dark matter and the increased weight 
due to the larger than Newtonian gravitational coupling. These MOG effects due to visible 
baryons are explored in greater detail in §5.3.3 as part of the analysis of the strong and 
weak lensing map of the Bullet Cluster 1E0657-558, presented in §5.3. 

Observations 

In the case of the Ursa Major sample of high and low surface brightness galaxies, the 
dynamic mass factors of §4.2.6 



r(r) = 



Gir) 
Gn 



(5.36) 
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Figure 5.2 continued: X-ray clusters of galaxies - Mass profiles. 

The dynamic mass factors, T(r), vs. orbital distance, r in kpc. 
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Figure 5.2 continued: X-ray clusters of galaxies - Dynamic mass factors. 

The dynamic mass factors, T(r), vs. orbital distance, r in kpc for a sample of X-ray clus- 
ters. The dynamic data consist of the ratio of the Newtonian dynamic mass of Equation 
(5.28), due to the measured isothermal temperature, to the integrated X-ray gas mass, 
derived from Equations (5.2) and (5.3) using the best-fit King /3-model parameters listed 
in Table 5.1. The computed best-fitted results are plotted for Moffat's STVG and MSTG 
theories and Milgrom's MOND theory with variable parameters. Results are plotted for 
the best-fit core-modified dark matter theory including the X-ray gas mass component. 
The reduced-x 2 statistic is included in Table 5.2. 
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are monotonically increasing, nearly linear functions, plotted in Figure 4.4. This is a 
prediction of the modified dynamics at small accelerations, of §3.2.2, where the slope is 
determined by Equation (3.20) to be the inverse of the transition radius, 



where a is the transition acceleration. 

For the best-fit cluster models of §5.2, the dynamic mass factors plotted in Figure 5.2, 
show very different trends, never showing a monotonically linear rise as in Equation (5.37). 

For each of the clusters of galaxies in the sample, T{r) 3> r for all r, having the greatest 
magnitude in the cores of the smaller (cooler) clusters, in particular Messier 49. The slope 



is close to flat for the larger (hotter) clusters, in particular the Bullet Cluster 1E0657- 
558 and Abell 2255, but generally having the greatest magnitude outside the cores of the 
smaller (cooler) clusters. All of the cluster cores have a particularly slow varying slope. 
Therefore Equation (5.38) suggests that clusters of galaxies are observationally inconsistent 
with singular (cuspy) models. For the Bullet Cluster 1E0657-558, the relationship between 
the X-ray observed S-map and the gravitational lensing ft-map is discussed in §5.3.2. 




(5.37) 



5.3 Bullet Cluster 1E0657-558 




The surface density S-map peaks 
reconstructed from X-ray imaging 
observations are shown in red and 
the convergence ft-map peaks as re- 
constructed from strong and weak 
gravitational lensing observations are 
shown in blue. Image provided cour- 
tesy of Chandra X-ray Observatory. 



Figure 5.3: Bullet Cluster 1E0657-558 
False colour composite image. 
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The Chandra Peer Review has declared the Bullet Cluster 1E0657-558 to be the most 
interesting cluster in the sky. This system, located at a redshift z = 0.296 has the highest 
X-ray luminosity and temperature (T = 14.1 ±0.2 keV ~ 1.65 x 10 8 K), and demonstrates 
a spectacular merger in the plane of the sky exhibiting a supersonic shock front, with 
Mach number as high as 3.0 ±0.4 (Markevitch, 2006). The Bullet Cluster 1E0657-558 has 
provided a rich dataset in the X-ray spectrum which has been modelled to high precision. 
From the extra-long 5.2 x 10 5 s Chandra space satellite X-ray image, the surface mass 
density, E(x, y) ) was reconstructed providing a high resolution map of the ICM gas (Clowe 
et al, 2007). The E-map, shown in a false colour composite map (in red) in Figure 5.3 is 
the result of a normalized geometric mass model based upon a 16' x 16' field in the plane of 
the sky that covers the entire cluster and is composed of a square grid of 185 x 185 pixels 
(~ 8000 data-points). 

Based on observations made with the NASA/ESA Hubble Space Telescope, the Spitzer 
Space Telescope and with the 6.5 meter Magellan Telescopes, Clowe et al (2006a); Bradac 
et al (2006); Clowe et al (2007) reported on a combined strong and weak gravitational 
lensing survey used to reconstruct a high-resolution, absolutely calibrated convergence k- 
map of the region of sky surrounding Bullet Cluster 1E0657-558, without assumptions on 
the underlying gravitational potential. The ft-map is shown in the false colour composite 
map (in blue) in Figure 5.3. The gravitational lensing reconstruction of the convergence 
map is a remarkable result, considering it is based on a catalogue of strong and weak lensing 
events and relies upon a thorough understanding of the distances involved - ranging from 
the redshift of the Bullet Cluster 1E0657-558 [z = 0.296) which puts it at a distance of 
the order of one million parsecs away. Additionally, the typical angular diameter distances 
to the lensing event sources [z ~ 0.8 to z ~ 1.0) are several million parsecs distant. 

In most observable systems, gravity creates a central potential, where the baryon den- 
sity peaks. As exhibited in Figure 5.3, the latest results from the Bullet Cluster 1E0657-558 
show, beyond a shadow of doubt, that the E-map, which is a direct measure of the hot 
ICM gas, is offset from the ft-map, which is a direct measure of the curvature (convergence) 
of space-time. The fact that the ft-map is centered on the galaxies, and not on the ICM 
gas mass is certainly either evidence of "missing mass" , as in the case of the dark matter 
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paradigm, or evidence of a stronger gravitational coupling due to a modification to gravity, 
as supported by Brownstein and Moffat (2007). Clowe et al (2007) stated 

One would expect that this (the offset E- and ft-peaks) indicates that dark 
matter must be present regardless of the gravitational force law, but in some 
alternative gravity models, the multiple peaks can alter the lensing surface 
potential so that the strength of the peaks is no longer directly related to the 
matter density in them. As such, all of the alternative gravity models have to 
be tested individually against the observations. 

Clowe et al (2007) described this as a degeneracy between whether gravity comes from 

dark matter, or from the observed baryonic mass of the hot ICM and visible galaxies where 

the excess gravity is due to a fifth force modification to the potential. This degeneracy 

may be split by examining a system that is out of steady state, where there is spatial 

separation between the hot ICM and visible galaxies. This is precisely the case in galaxy 

cluster mergers such as the Bullet Cluster 1E0657-558, since the galaxies will experience a 

different gravitational potential created by the hot ICM than if they were concentrated at 

the center of the ICM. 

The data from the Bullet Cluster 1E0657-558 provides a laboratory of the greatest scale, 

where the degeneracy between "missing mass" and "extra gravity" may be distinguished. 

We are fortunate, indeed, that the Bullet Cluster 1E0657-558 is not only one of the hottest, 

most supersonic, most massive cluster mergers seen, but the plane of the merger is aligned 

with our sky! Brownstein and Moffat (2007) addressed the full-sky data product (Clowe 

et al, 2006b) for the Bullet Cluster 1E0657-558, and provide first published results for 

the King /5-model of the E-map, detailed in §5.3.1. Brownstein and Moffat (2007) utilized 

the metric skew-tensor gravity model of §3.3 to compute component mass profiles, and 

account for all of the baryons in each of the main and subclusters, including all of the 

galaxies in the regions near the main central dominant (cD) and the subcluster's brightest 

central galaxy (BCG), without non-baryonic dark matter, to fit the gravitational lensing 

convergence map, as in §5.3.2. The results of the analysis include a map of the visible 

baryon distribution, as in §5.3.3, and the dark matter distribution, as in §5.3.4. 
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Figure 5.4: Bullet Cluster 1E0657-558 - X-ray gas surface density map. 

Data reconstructed from X-ray imaging observations of the Bullet Cluster 1E0657-558, 
November 15, 2006 data release (Clowe et al, 2006b), showing E-map observed peaks 
(local maxima) and ft-map observed peaks. The central dominant (cD) galaxy of the 
main cluster, the brightest cluster galaxy (BCG) of the subcluster, and the MOG 
predicted gravitational center are shown. The colourscale is shown at the bottom, in 
units of 10 15 M /pixel 2 . The resolution of the E-map is 8.5 kpc/pixel, based upon the 
measured redshift distance ~ 260.0 kpc/arcminute (Bradac et al, 2006). The scale in kpc 
is shown at the top of the figure. J2000 and map (x,y) coordinates are listed in Table 5.3. 

5.3.1 X-ray gas map 

With an advance of the Clowe et al (2006b) November 15, 2006 data release, Brownstein 
and Moffat (2007) performed a precision analysis to model the gross features of the surface 
density S-map data in order to gain insight into the three-dimensional matter distribution, 
p(r), and to separate the components into a model representing the main cluster and the 
subcluster - the remainder after subtraction. 

The E-map is shown in false colour in Figure 5.4. There are two distinct peaks in the 
surface density S-map - the primary peak centered at the main cluster, and the secondary 
peak centered at the subcluster. The main cluster gas is the brightly glowing (yellow) 
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Table 5.3: J2000 sky coordinates of the Bullet Cluster 1E0657-558 



Observation 




J2000 Coordinates 


E-MAP 


K-MAP 






RA 


Dec 




0, y) 


(x, y) 


(1) 




(2) 






(3) 




(4) 


(5) 


Main cluster E-map peak 


06 


58 


31.1 


-55 


56 


53.6 


(89, 89) 


(340, 321) 


Subcluster E-map peak 


06 


58 


20.4 


-55 


56 


35.9 


(135, 98) 


(365, 326) 


Main cluster K-map peak 


06 


58 


35.6 


-55 


57 


10.8 


(70, 80) 


(329, 317) 


Subcluster K-map peak 


06 


58 


17.0 


-55 


56 


27.6 


(149, 102) 


(374, 327) 


Main cluster cD 


06 


58 


35.3 


-55 


56 


56.3 


(71,88) 


(330, 320) 


Subcluster BCG 


06 


58 


16.0 


-55 


56 


35.1 


(154,98) 


(375, 326) 


MOG Center 


06 


58 


27.6 


-55 


56 


49.4 


(105,92) 


(348, 322) 



Notes. — November 15, 2006 data release (Clowe et a/., 2006b): Column (1) provide the primary 
observational features. Columns (2) and (3) list the J2000 right ascension (RA) and declination 
(Dec) for each feature. Columns (4) and (5) provide the S-map and ft-map (x, y) coordinates 
using a resolution of 8.5 kpc/pixel, and 15.4 kpc/pixel, respectively, based upon the measured 
redshift distance ~ 260.0 kpc/arcminute of the Bullet Cluster 1E0657-558 (Bradac et al, 2006). 



region to the left of the subcluster gas, which is the nearly equally bright Shockwave region 
(arrowhead shape to the right). The ft-map observed peaks, the central dominant (cD) 
galaxy of the main cluster, the brightest cluster galaxy (BCG) of the subcluster, and the 
MOG predicted gravitational center are shown in Figure 5.4 for comparison. J2000 and 
map (x,y) coordinates are listed in Table 5.3. 

Since there is a multitude of source galaxies in a range of redshifts (z = 0.85 ± 0.15), 
the source distances, D S) may be averaged. For the Bullet Cluster 1E0657-558, Clowe et al. 
(2004) used 

J1, n ' 540 kpc, (5.39) 



l^ls 



and the Newtonian critical surface mass density (with vanishing shear) of Equation (4.28), 

c 2 A 



47TG N A As 



3.1 x 10 9 M /kpc' 



(5.40) 



is effectively constant. 
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Figure 5.5: Bullet Cluster 1E0657-558 
King /3-model fit to scaled E-map. 

A cross-section of the E-map of 
Figure 5.4 reconstructed from X-ray 
imaging observations (Clowe et al, 
2006b), shown in solid red, on a 
straight-line connecting the main X-ray 
cluster peak to the main central dom- 
inant (cD) galaxy. The King /3-model 
(neglecting the subcluster) of Equation 
(5.25) is shown in short-dashed blue, 
best-fit by Equations (5.41), (5.42) 
and (5.43). The unmodeled peak (at 
R ^ 300 kpc) is due to the subcluster. 
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King /?-model of the main cluster 

To calculate the best-fit parameters, /?, r c and p of the King /3-model of Equations (5.24) 
and (5.25), Brownstein and Moffat (2007) applied a nonlinear least-squares fitting routine 
(including estimated errors) to the E-map on a straight-line connecting the main cluster E- 
map peak (R = 0) to the main cD, and then extrapolated the fit to the entire map. This 
reduces the complexity of the calculation to a simple algorithm, but is not guaranteed 
to yield a global best-fit. However, the approximation provides a very low reduced x 2 
everywhere on the full sky map, except at the subcluster (which is masked for the best- 
fit). The X-ray gas surface density E-map data, and the King /3-model of Equation (5.25), 
best-fit to the scaled E-map, are shown in Figure 5.5, with the best-fit parameters, 

P = 0.803 ±0.013, (5.41) 
r c = 278.0 ± 6.8 kpc, (5.42) 



where the value of the E-map at the main cluster peak is constrained to the observed value, 

E = 1.6859 x 10 10 ^% ( a \Z X f V = 2.3181 x lO 8 M /kpc 2 , (5.43) 
pixel \ 8.528 kpc / 
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scaled by S c of Equation (5.40). Solving Equation (5.24) for the central density of the 
main cluster, 

po = \ 2 J = 0.334 x 10 6 M /kpc 3 , (5.44) 

which is between one and two orders of magnitude less than the dark matter central 
densities listed in Table 4.2 derived from the galaxy rotation curves of §4.2.3, proving that 
cluster scale dark matter does not affect the dynamics of galaxy rotation curves. The set 
of parameters, /?, r c and p ? completely determines the isotropic isothermal King /J-model 
for the density, p(r), of Equation (5.2) of the main cluster X-ray gas, and the isotropic 
isothermal model of §5.1.2 may be applied to measure the mass-luminosity relation in the 
main cluster and compute the ratio of the Newtonian dynamic mass to the X-ray gas 
(baryon) mass, per gravity theory. 

Brownstein and Moffat (2007) computed the Newtonian dynamic mass profile for the 
main cluster of the Bullet Cluster 1E0657-558, and determined the MSTG mass profile 
according to Equations (3.93) and (3.94), finding that the modified gravity mass profile 
is an excellent fit to the measured X-ray (baryon) mass profile, as shown in Panel (a) of 
Figure 5.1. Across the full range of the r-axis, and throughout the radial extent of the 
Bullet Cluster 1E0657-558, the la correlation between the gas mass, M(r) and the MOG 
dynamical mass, MmstgMj provides excellent agreement between theory and experiment. 

Substituting Equations (5.41), (5.42) and (5.44) into Equation (5.4), we obtain the 
main cluster outer radial extent, 

r out = 2620 kpc, (5.45) 

the distance at which the density, p(r out ), drops to w 10 _28 g/cm 3 , or 250 times the mean 
cosmological density of baryons. The total mass of the main cluster may be calculated by 
substituting Equations (5.41), (5.42) and (5.44) into Equation (5.5): 

M gas = 3.87 x 10 14 M , main cluster. (5.46) 

The MOG temperature prediction, from the MSTG best-fit, is increasingly consistent 
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with updated experimental values, shown in Table 5.4. 



Table 5.4: Isothermal temperature of the main cluster 







Year 


Source - Theory or Experiment 


T (keV) 


% ERROR 


2007 
2002 
1999 
1998 


Computed value 
accepted experimental value 
ASCA+ROSAT fit 
ASCA fit 


15.5 ±3.9 

14.5+f 
17.4 ± 2.5 


4.5 

6.5 
12.3 



Notes. — The computed isothermal temperature is consistent with the experimental values for 
the main cluster (Markevitch et al, 2002).: Column (1) and (2) list the year and source of the 
temperature result, respectively. Column (3) provides the temperature in keV, and Column (4) 
provides the percent error between the computed and experimental values. 

Model of the bullet subcluster 

Although the X-ray morphology of the main cluster is very regular, and well described by 
the King /3-model of the main cluster, Liang et al (2000) reported on a diffuse radio halo, 
which requires the acceleration of thermal electrons to ultra-relativistic energies, enhanced 
at the main X-ray gas peak and more focused at the densest part of the optical galaxy 
distribution. Since galaxies are collisionless, at the ~ 1 Mpc cluster scale, a merger with 
the subcluster - the bullet in the X-ray gas surface density map of Figure 5.4 - allows the 
galaxies to stream through the X-ray gas and generate the radio halo. 

Markevitch et al. (2002) reported on Chandra observations, providing evidence that 
the merger is ongoing and the subcluster is in a perturbed state far from hydrostatic 
equilibrium leading to an apparent increase in the X-ray temperature, 150 million years 
after its collision with the main cluster core. Barrena et al (2002) studied the dynamics of 
the collision, and determined that the subcluster is the remnant core of a moderate mass 
cluster of galaxies, whose properties have been strongly affected. Randall et al (2008) 
studied the prominent bow shock, estimating the supersonic merger velocity of 4700 km 
s _1 , and concluded that the subcluster X-ray gas mass distribution is significantly more 
peaked than a King profile. 

Brownstein and Moffat (2007) computed the surface mass density of the subcluster 
by subtracting the best-fit (x 2 < 0.2) King /3-model to the main cluster - which agreed 
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(a) The subcluster surface density S-map (b) The subcluster E-map su- 

perposed onto the best-fit King /3-model of the main cluster 

Figure 5.6: Bullet Cluster 1E0657-558 - X-ray gas surface density model. 

(a) Subtracted subcluster X-ray gas surface density map is shown in green. 

(b) X-ray gas surface density model. The blue surface represents the E- 
map due to the integrated (line-of-sight) King /?-model fit to main clus- 
ter. The green surface is the contribution to the E-map from the subcluster. 

with the main cluster surface mass E-map (data) within 1% everywhere - from the total 
X-ray surface mass density of Figure 5.4. The subcluster subtraction is accurate down to 
p = 10 -28 g/cm 3 ~ 563.2 M /pc 3 baryonic background density. After subtraction, the 
subcluster E-map peak takes a value of 1.30 x 10 8 M /kpc 2 , whereas the full E-map has 
a value of 2.32 x 10 8 M /kpc 2 at the subcluster E-map peak. Thus the subcluster (at 
its most dense position) provides only w 56% of the X-ray ICM, the rest is due to the 
extended distribution of the main cluster. 

Figure 5.6 is a stereogram of the subcluster subtracted surface density E-map and the 
subcluster superposed onto the surface density E-map of the best-fit King /?-model to the 
main cluster. 

Since the outer radial extent of the subcluster gas is less than 400 kpc, the E-map 
completely contains all of the subcluster gas mass. By summing the subcluster subtracted 
E-map pixel-by-pixel over the entire E-map peak, one is performing an integration of 
the surface density, yielding the total subcluster mass. Brownstein and Moffat (2007) 
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performed such a sum over the subcluster subtracted E-map data, obtaining 

M gas = 2.58 x 10 13 M©, subcluster. (5.47) 

for the mass of the subcluster gas, which is less than 6.7% of the mass of main cluster gas 
of Equation (5.46). This justifies the initial assumption that the subcluster may be treated 
as a perturbation in order to fit the main cluster to the King /?-model. The subsequent 
analysis of the thermal profile confirms that the main cluster X-ray temperature is nearly 
isothermal, lending further support to the validity of the King /?-model and the reliability 
of the isothermal temperatures of Table 5.4. 

5.3.2 Gravitational lensing convergence map 

The convergence ft-map of Figure 5.7 is a false colour image of the strong and weak gravi- 
tational lensing reconstruction (Clowe et al, 2006a; Bradac et a/., 2006; Clowe et a/., 2007) 
of all of the bending of light over the entire distance from the lensing event source toward 
the Hubble Space Telescope. The source of the ft-map is oc / G^p(r), along the line-of-site, 
as in the Newtonian case of Equation (4.27), but oc / G{r)p(r) as in Equation (4.29) of 
modified gravity with a spatially varying gravitational coupling. 

Modified gravity solution 

The lack of spherical symmetry in the ft-map, shown in Figure 5.7, is better visualized in 
Panel (a) of Figure 5.8, which demonstrates the importance of the subcluster's dynamic 
mass. 

Brownstein and Moffat (2007) utilized the metric skew-tensor gravity model of §3.3 to 
compute the weighted surface mass density, E of Equation (5.50), of the X-ray gas mass of 
the main cluster using the King /3-model of Equation (5.2) with the best-fit parameters of 
Equations (5.41), (5.42) and (5.44). This is shown as E(r)/E c by the black surface of the 
ft-model of Panel (b) of Figure 5.8. Including the galaxies is accomplished by Equation 
(5.48) which is shown by the red surface of the ft-model. 

Brownstein and Moffat (2007) proceeded to account for the spherical symmetry break- 
ing effect of the subcluster on the dynamic mass of Equation (3.94): Remarkably, as the 
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Figure 5.7: Bullet Cluster 1E0657-558 - Gravitational lensing convergence map. 

Data reconstructed from strong and weak gravitational lensing of the Bullet Clus- 
ter 1E0657-558, November 15, 2006 data release (Clowe et al, 2006b), showing 
convergence ft-map observed peaks (local maxima) and E-map observed peaks. The 
central dominant (cD) galaxy of the main cluster, the brightest cluster galaxy (BCG) of 
the subcluster, and the MOG predicted gravitational center are shown. The scale in kpc 
is shown at the top of the figure. J2000 and map (x,y) coordinates are listed in Table 5.3. 



MOG center was separated from the main cluster E-map peak, due to the gravitational 
effect of the subcluster, the centroid naturally shifted toward the ft-map peak, and the 
predicted height of the ft-map decreased, flattening the peak and dimpling the core and 
skewing the distribution in the direction opposite to the shift in the MOG center. Although 
Moffat (2006a) demonstrated that the integration of the ft-map, assuming a constant sur- 
face mass density for the galaxies, produced a peak offset from the X-ray peak, the effect 
alone was insufficient to fit the Bullet Cluster 1E0657-558 ft-map data. However, the 
difference can be entirely accounted for by including the surface mass density of the galax- 
ies, S ga i ax (r)/S c , as indicated by the red surface of the best-fit k model of Figure 5.8. 
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(a) ft-map, November 15, 2006 data release (b) Best-fit ft-model 

Figure 5.8: Bullet Cluster 1E0657-558 - Best-fit model to the gravitational lensing convergence map. 

(a) Gravitational lensing convergence ft-map November 15, 2006 data release (Clowe 
et al : 2006b) is shown in gold. The twin peaks are due to the main and subcluster, 
respectively, (b) Gravitational lensing convergence model. The black surface is the 
best-fit ft-map to the main cluster X-ray gas component. The red surface represents 
the excess ft-map of the galactic component beyond the best-fit X-ray gas component. 

Combining the black surface and the red surface, we obtain the best-fit model, 



K,(r) 



E(r) + S ga i ax (r) 



(5.48) 



which is equivalent to the ft-map data illustrated by the gold surface on the left hand side 
of Figure 5.8. 

As introduced in §4.2.4, predictions for the ft-map of high resolution sub-kiloparsec 
galaxy-galaxy lensing, plotted in Figure 4.2, are computed by 



K(r) = J 



(5.49) 



where 



G(r) 
Gn 



p{r)dz, 



(5.50) 



is the weighted surface mass density. For the multiple source Bullet Cluster 1E0657-558 
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reconstruction, Clowe et al. (2004) used E c « 3.1 x lO 9 M /kpc 2 , without estimate of the 
uncertainty. The precision of the ft-map depends on the validity of the assumption of small 
variation AE C across the lens, which depends on the variation in the ratio, D\ s /D s . 
Substituting Equation (5.36) into Equation (5.50), 



In the Newtonian limit, G(r) —> G/v, and therefore the factor T(r) — > 1 is removed from 
the integral, and E —> E, recovering the E-map of Equation (5.16). 

Whereas T(r) is a maximum value for the outer radial positions of galaxies, contributing 
more weight to the integral of Equation (5.51) in the galactic halo, T(r) contributes more 
weight to the cores of spherically symmetric clusters of galaxies. 

5.3.3 Visible baryon distribution 

The galaxies contribute a "measurable" surface mass density, 



which we may interpret as the difference between the ft-map and the scaled contribution 
from the weighted surface density of the ICM gas. The result of the galaxy subtraction 
of Equation (5.52) is shown as the galactic surface mass density map, in Figure 5.9. The 
surface mass density of the visible baryons is taken to be the sum of the ICM gas component 
and the galaxies, as shown in the left panel of Figure 5.10. 

The total mass of the galaxies is determined by integrating over the E-map, 



Brownstein and Moffat (2007) performed the integration within a 100 kpc radius aper- 
ture about the main cluster cD and subcluster BCG, separately, the results of which are 
listed in Table 5.5, where they are compared with the upper limits on galaxy masses set 
by HST observations. If the hypothesis that the predicted M ga i ax is below the bound set 




(5.51) 




(5.52) 




(5.53) 
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Figure 5.9: Bullet Cluster 1E0657-558 - Galactic surface density map. 

The predicted S-map due to the galaxies as computed by the difference between 
the ft-map and the MOG ft-model, scaled as surface mass density according to Equation 
(5.52). E-map and ft-map observed peaks are shown for comparison. The central 
dominant (cD) galaxy of the main cluster, the brightest cluster galaxy (BCG) of the 
subcluster, and the MOG predicted gravitational center are shown. J2000 and map (x,y) 
coordinates are listed in Table 5.3. Component masses (integrated within a 100 kpc 
radius aperture) for the main and subcluster, the MOG center and the total predicted 
baryonic mass, Mb aryon , for the Bullet Cluster 1E0657-558 are shown in Table 5.5. 

by HST observations is true, then it follows that 



M bary = M t 



gas 



M, 



galax ) 



(5.54) 



requires no addition of non-baryonic dark matter. The results of our best-fit for M gas , M ga i ax 
and Mb ar y of Equation (5.54) are listed in Table 5.5. The result of M ga i ax /M gas w 0.4% in 
the central ICM is due to the excellent fit in MOG across the hundreds of kpc separating 
the main and subcluster. 
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predicted galaxies, S ga i ax /E c , and the MOG predicted density of dark matter, Hdm/S c , compared to ICM gas. 

visible baryons, Ebary/^c, compared to the ICM gas. 

Figure 5.10: Bullet Cluster 1E0657-558 - Distribution of visible and dark matter 

Plot of the scaled surface density, E(r)/E c , of the ICM gas is shown in short- 
dashed green on each plot along the line connecting the main cluster E-map peak with 
the main cD. In the left panel, the MOG prediction for the galaxies of Equation (5.52) 
is shown in long-dashed magenta, and the prediction of the visible baryonic mass due 
to the combined ICM gas mass and galaxies is shown in solid brown. The calculation 
of Equation (5.55) for dark matter is plotted in the right panel in dash-dot black. 

5.3.4 Dark matter distribution 

In the absence of modified gravity, dark matter is hypothesized to account for all of the 
"missing mass" which results in applying Newton/Einstein gravity. This means, for the 
November 15, 2006 data release (Clowe et al, 2006a; Bradac et al, 2006; Clowe et al, 
2007, 2006b), that the "detected" dark matter must contribute a surface mass density, 



E D m(^, y) ~ k(x, y)E c - E(x, y), 



(5.55) 



and is plotted in the right panel of Figure 5.10. 

The total mass of the dark matter distribution with an associated total mass, 



M DM = / E BM (r)dxdy. 



(5.56) 
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Upon substitution of Equation (5.55), the integral of Equation (5.56) becomes: 

where we have neglected M ga i ax in Equation (5.57), because the contribution from the 
galaxies is < 1 - 4% of M tota i due to dark matter dominance. 

Brownstein and Moffat (2007) computed M DM in Equation (5.57) by a pixel- by-pixel 
sum over the convergence ft-map data and surface density E-map data, within a 100 kpc 
radius aperture around the main and subcluster ft-map peaks, respectively. The result of 
the computation, including the mass ratios, M ga i ax /M gas , for the main and subcluster and 
central ICM are provided in Table 5.5. 



Table 5.5: Component mass predictions. 



Component 


Main cluster 


Subcluster 


Central ICM 


Total 


(i) 


(2) 


(3) 


(4) 


(5) 


M gas 


7.0 x 10 12 M Q 


5.8 x 10 12 M Q 


6.3 x 10 12 M Q 


2.2 x 10 14 M Q 




1.8 x 10 12 M 


3.1 x 10 12 M© 


2.4 x 10 10 M 


3.8 x 10 13 M 


Mjary 


8.8 x 10 12 M 


9.0 x 10 12 M 


4.9 x 10 12 M Q 


2.6 x 10 14 M 




2.1 x 10 13 M© 


1.7 x 10 13 M 


1.4 x 10 13 M 


6.8 x 10 14 M 


M galax /M gas 


26% 


53% 


0.4% 


17% 


M gas /M DM 


33% 


34% 


45% 


32% 



Notes. — Column (1) specifies the component masses and mass fractions. Columns (2) and (3) list 
the component masses integrated within a 100 kpc radius aperture for the main and subcluster, 
respectively. Columns (4) lists the component mass integrated within a 100 kpc radius aperture 
for the central ICM located at the MOG center. Column (5) lists the total of each component 
masses integrated over the full H-map. 



The dark matter result of M gas /MvM ~ 45% in the central ICM implies that the 
evolutionary scenario does not lead to a spatial dissociation between the dark matter 
and the ICM gas, which confirms that the merger is ongoing. In contrast, the MOG 
result shows a true dissociation between the galaxies and the ICM gas as required by 
the evolutionary scenario. The baryon to dark matter fraction over the full S-map is 
32%, which is significantly higher than the A-CDM cosmological baryon mass- fraction of 
17±i;l% (Spergel et a/., 2007). 
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5.3.5 Neutrino halos 

Sanders (2003) postulated a two component model for the Coma cluster, adding a nonlu- 
minous rigid sphere to include the contribution of finite mass neutrinos - a candidate for 
hot non-baryonic dark matter - with a constant density core, 

^.SxlO-gcm-^Q 1 . (5.58) 

Sanders (2003) assumed that the constant density cores have finite radii that scale as 

r v = 2r c , (5.59) 

where r c is the gas core radius of the isothermal King /3-model of Equation (5.2). For a 2 
eV neutrino and the accepted experimental value of Table 5.4 suggest 

p v < 3.9t°ol x lO 6 M kp C - 3 . (5.60) 

Substituting Equation (5.42) into Equation (5.59) gives the constant density neutrino core 
radius, 

r v = 556 ± 14kpc, (5.61) 
which has an integrated mass within an aperture of 100 kpc of 

M 100 = l-6±g;l x 1O 13 M , (5.62) 

which is within 30% of the required value according to Table 5.5. However, the total 
integrated core mass of 

M y = 2.8i°J x 1O 15 M , (5.63) 

exceeds the mass of the Bullet Cluster 1E0657-558 by a factor of three, implying that 
while Equation (5.60) is reasonable, Equations (5.59) and (5.61) may be overestimating 
the extent of the neutrino halos by a factor of two. 

Sanders (2007) elaborated on the MOND neutrino-baryon model of clusters, confirming 
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the result of §5.2.4 that the need for dark matter appears to decrease with increasing 
temperature, suggesting Equation (5.58) is opposite to observation, but argued that the 
observed trend is caused by the cooling and inflow of baryons. Sanders (2007) provided 
an improved neutrino halo model for clusters, with the constant density core of Equation 
(5.58), and included a theoretically derived scaling relation instead of Equation (5.59), 



O.TMpc^)" 7 forT<3keV( 



r„= V - 2 —i " V ' (5-64) 



2eVy 



l.lMpc(f-) 2 (^) 1 otherwise, 



which assumes three flavours of neutrinos, each of which have comparable velocity disper- 
sions to the baryons and maintain their cosmological density ratio, 

^Maryon = 2 -8 (^J . (5.65) 

Now, substituting the accepted experimental value of Table 5.4 of T = 14.5^7 keV into 
Equation (5.64) gives the constant density 2 eV neutrino core radius, 

TV = 575t?? kpc, (5.66) 

which is consistent with Equation (5.61) and therefore too large by a factor of 2 to explain 
the total mass of the Bullet Cluster 1E0657-558. 

Angus et al. (2007) confirmed that a simple model of 4 dominant constant density cores 
of 2 eV neutrinos can supply the missing mass in the peaks of the gravitational lensing 
convergence ft-map of the Bullet Cluster 1E0657-558, provided the neutrino cores have 
radii r v < 50 kpc. Brownstein and Moffat (2007) provided a comparison of the surface 
density of dark matter to the surface density of the X-ray emitting ICM gas, shown in 
Panel (b) of Figure 5.10, consistent with 2 extended, overlapping halos centered at the 
galactic regions, which may have constant density cores in the inner 50 kpc, but then 
declining more rapidly. 

In the absence of scaling relations, such as those investigated by Sanders (2003), each 
neutrino halo requires additional free parameters to specify the shape of the density pro- 
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file, which may be better described by other possibilities including the King /3-model, or 
the core-modified dark matter profile of §5.2.1 which fits clusters of galaxies without the 
necessity of MOND. The core-modified dark matter fit to the main cluster of the Bullet 
Cluster 1E0657-558 from §5.2, with best-fit parameters listed in the top row of Panel (a) 
of Table 5.2, provides 

po = (2.73 ± 0.47) x lO 6 M kp C - 3 , (5.67) 

r s = 328.1 ±26.7kpc, (5.68) 

Identifying Equation (5.67) with Equation (5.58), we may solve for the upper limit on the 
neutrino mass, 

m v < 1.8±0.1eV, (5.69) 

which is below the Mainz/Troitsk experimental limit on the electron neutrino, m^ e < 2.2 
eV, but is falsifiable in the near future. 

Angus et al (2008) decomposed the mass profiles of 26 X-ray systems according to 
MOND, with temperatures ranging from 0.5 keV to 9 keV, and concluded that whatever 
the equilibrium distribution, 2 eV neutrino halos cannot explain the inner 100 to 150 
kiloparsecs of clusters within MOND. This issue is seen in the dynamic mass factors plotted 
in Figure 5.2, since each plot is maximized in the inner region of every cluster in the sample, 
where the missing mass problem is most pronounced. 
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Solar system 

"I was like a boy playing on the sea-shore, and diverting myself now and then 
finding a smoother pebble or a prettier shell than ordinary, whilst the great 
ocean of truth lay all undiscovered before me. " 

Sir Isaac Newton 



The motion of the planets and planetoids, their satellites, and the chunks of matter 
that comprise the asteroids and the comets are along paths derived from matter's response 
to gravity. The opportunity to discover new celestial physics in the solar system provides a 
challenge to form deeper understandings of Kepler's eponymous laws, from which Newton's 
theory of universal gravitation is founded. Precise observation of orbits of the many bodies 
in the solar system suggest Kepler's three laws require subtle corrections: 
Kepler's first law: The path of planets and bodies about the sun are near elliptic 
in shape, with a focus near the center of the sun, but changing in time under the 
influence of Jupiter and the other solar bodies. 
Kepler's SECOND law: An imaginary line drawn from the center of a body to the center 
of a body in orbit will sweep out nearly equal changing areas in equal intervals of 
time, where the change in area slightly increases if orbital angular momentum is 
transferred to the orbiting body from the spin of the central body, and decreases if 
angular momentum is transferred in the opposite direction. 
Kepler's third law: The ratio of the squares of the periods of any two planets is nearly 
equal to the ratio of the cubes of their average distances from the sun, where the 
difference in this near equality is most significant at the orbit of Jupiter. 
Because the sun is not the only source of gravity in the solar system, and since so few 
solutions to Einstein's gravity theory are known, modelling gravity in the solar system is 
a managed process, such as the astronomer's ephemerides, which are datacentric solutions 
without the elegance and utility of a theoretical prediction, which does not need daily 
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updates to correct for unmodelled physics, deemed unnecessary. Jupiter adds a significant 
source of gravity to the solar system, with mass = O.OO95M , which is more than 
twice the total mass of all the other smaller bodies, combined. 

New physics beyond the orbit of Jupiter must be nearly consistent with Kepler's three 
laws, and should make quantitative predictions of the necessary amendments. However, 
all terrestrial and solar system attempts to falsify Moffat's nonsymmetric gravity theory 
(NGT) have led only to upper bounds on the possible strength of the modified gravity 
fifth force, including predictions for the Gravity Probe B experiment. Moffat and Brown- 
stein (1990) considered spinning test particles and the motion of a gyroscope, finding that 
the difference between the NGT correction to the gyroscope precession, and the Einstein 
correction, would be smaller than the Gravity Probe B experiment could detect in orbit 
about Earth. 

Brownstein and Moffat (2006c) considered the motion of the Pioneer 10 and 11 space- 
craft in the metric skew-tensor gravity theory, as in §3.3, proving that the unexpected 
sunward acceleration can be explained by modified gravity without leading to disagree- 
ment between the predicted and actual orbits of the outermost planets. According to the 
Pioneer Explorer Collaboration, the most likely explanation is that there is a systematic 
origin to the effect, such as a thermal recoil force investigated by Toth and Turyshev (2009), 
using a simulated Pioneer 10 data set, but neither has NASA ruled out the modified gravity 
solution, presented in §6.1.2. 

6.1 Pioneer 10/11 Anomaly 

The radio tracking data from the Pioneer 10/11 spacecraft during their travel to the outer 
parts of the solar system have revealed a possible anomalous acceleration. The Doppler 
data obtained at distances r from the Sun between 20 and 70 astronomical units (AU) 
showed the anomaly as a deviation from Newton's and Einstein's gravitational theories. 
At this time, NASA continues to support the search for a gravitational solution, as in 
§6.1.2, but the Pioneer Explorer Collaboration may eventually be able to rule out modified 
gravity as the origin of the effect, once the recovered data sets have been formatted and 
a comprehensive model can be applied, as progressing according to Toth and Turyshev 
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(2009). 

Brownstein and Moffat (2006c) applied the metric skew-tensor gravity theory of §3.3, in 
which Einstein gravity is coupled to a Kalb-Ramond-Proca field, as in §3.3.2, and provided 
a fit to the available anomalous acceleration data or the Pioneer 10/11 spacecraft consistent 
with all current satellite, laser ranging and observations for the inner planets. 

The Pioneer anomalous acceleration observations are described in §6.1.1, and the fit is 
presented in §6.1.2. The effect of modified gravity in the solar system on Kepler's law of 
motion and the planetary ephemerides are explored in §6.1.3 and §6.1.4, respectively, and 
the constraints set by observations of the anomalous perihelion advance are identified in 
§6.1.5. 

6.1.1 Pioneer anomalous acceleration 

Anderson et al (1998, 2002); Turyshev et al (2006a) observed the Doppler residuals data as 
the differences of the observed Doppler velocity from the modelled Doppler velocity, and 
computed the anomalous acceleration directed towards the Sun, with an approximately 
constant amplitude over the range of distance, 20 AU < r < 70 AU: 

a P = (8.74 ± 1.33) x lO^cm^ 2 . (6.1) 

After a determined attempt to account for all known sources of systematic errors, 
Anderson et al (1998, 2002); Turyshev et al (2006a) reached the conclusion that the 
Pioneer anomalous acceleration towards the Sun could be a real physical effect that requires 
a physical explanation. Turyshev et al (2006b) reviewed NASA's efforts to recover the 
extended Pioneer doppler data set, emphasizing that the apparent difficulty to explain the 
anomaly within standard physics is a motivation to look for new physics, including the 
model of Brownstein and Moffat (2006c). 

In NASA's official statement, Turyshev (2007) reported, 

"As of March 2007, the existence of the anomaly is confirmed by seven 
independent investigations using different navigational codes - the signal is 
present in the Doppler data received from both Pioneers 10 and 11. The most 
important question now is to identify the cause of this anomalous frequency 
drift discovered in the Pioneer data. 



192 



Chapter 6. SOLAR SYSTEM 



. . Our thermal modelling of the Pioneer vehicles is progressing very well. 
We finished the development of the geometric mathematical models of the 
spacecraft that include geometry and properties of most of the important space- 
craft components and surfaces. We are now working on the thermo-dynamical 
model of the vehicles. At this stage, we have a very good understanding of 
heat re-distribution within the craft and soon will be ready to compute the 
heat flow to the outside of the craft. Soon, we will be able to tell whether or 
not heat contributes to the formation of the anomaly." 

MOND is not considered a viable mechanism because the value of the MOND universal 

acceleration of Equation (4.8) that provides good fits to galaxy rotation curves, as in 

Chapter 4, is orders of magnitude smaller than the acceleration of the Pioneer satellites, 

ap, until the satellite reaches the MOND transition radius of Equation (3.21), 



V &o 

and is not likely to be observed on the scales of the solar system. 

Galaxy scale dark matter cannot affect the solar system, since the density of the Milky 
Way dark matter halo in the vicinity of the solar system is 



and therefore a galaxy scale dark matter globe, equivalent in mass to Earth, would have a 
radius of greater than 10,000 AU. 

6.1.2 Gravitational solution 

The acceleration law of Equations (3.90) and (3.91), derived from the metric skew-tensor 
gravity theory of §3.3.5, can be written 




(6.2) 



5.0 x 1O" 19 M /AU 3 , 



(6.3) 



a(r) 



G(r)M 



(6.4) 



where 




(6.5) 
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Brownstein and Moffat (2006c) postulated a gravitational solution that the Pioneer 10/11 
anomaly is caused by the difference between the running G(r) of Equation (6.5) and the 
bare value, Gn- So the Pioneer anomalous acceleration directed towards the center of the 
Sun is given by 



5G(r)M e 







where 



SG(r) = G(r) - G N = G N a(r) 
The dynamic mass factor is defined as 



1 - exp(-r/A(r)) 1 + 



A(r), 



r(r) = G(r)/G N = 1 + 



6G(r) 
G 



N 



G(r)/G N = 1 + a(r) 



1 - exp(-r/A(r)) 1 + 



A(r) 



(6.6) 



(6.7) 



(6.8) 
(6.9) 



and measures the degree to which the observed acceleration of the Pioneer satellite differs 
from the Newtonian acceleration, 



a P (r) = (r(r) - l)a N (r), 



(6.10) 



where the Newtonian acceleration is 



a N (r) 



G N M 



(6.11) 



Therefore, a measurement of T{r) ~ 1 dismisses the Pioneer anomaly, whereas modified 
gravity predicts a monotonically increasing T{r) due to Equation (6.9). 

Brownstein and Moffat (2006c) proposed the following parametric representations of 
the running of a(r) and A(r): 



<*(r) ^(1 - exp(-r/r)) 6/2 , 



A(r) 



(1 — exp(— r/r)) b 



(6.12) 
(6.13) 
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Figure 6.1: Pioneer 10/11 anomalous acceleration 

Pioneer 10/11 anomalous acceleration, ap(r) in 10 -8 cm s _1 , vs. orbital distance, r 
in AU, of Equation (6.6), is plotted on a linear scale out to r = 50 AU, in the 
left panel, and on a logarithmic scale out to r = 5,000AU, in the right panel. 



Theoretical best-fit 


with lines 


Data 


with error bars 


a P (r) 


red solid line 


Pioneer 10 
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mean dp 


horizontal black dotted lines 


Pioneer 11 


blue closed circles 



Here, f is a non-running distance scale parameter and b is a constant. 

In Figure 6.1, we display a best-fit to the Pioneer 10/11 anomalous acceleration data 
from Nieto and Anderson (2005, Figure 4) obtained using a nonlinear least-squares fitting 
routine including estimated errors from the Doppler shift observations (Anderson et a/., 
2002). 

The best-fit parameters are: 

= (1.00 ± 0.02) x 10" 3 , 
= 47±1AU, 
f = 4.6±0.2AU, 

b = 4.0. (6.14) 



The small uncertainties in the best-fit parameters are due to the remarkably low variance 
of residuals corresponding to a reduced x 2 per degree of freedom of 0.42 signalling a good 
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Figure 6.2: Dynamic mass factor in the Solar system 

Pioneer 10/11 anomaly dynamic mass factor, T(r), vs. orbital distance, r in AU, 
of Equation (6.8), is plotted on a linear scale out to r = 50 AU, in the left 



panel, and on a loj; 


^arithmic scale out to r 


= 5,000AU, 


in the right panel. 


Theoretical best-fit 


with lines 


Data 
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Pioneer 10 
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Pioneer 11 


blue closed circles 



fit. 

In Figure 6.2, we display the Pioneer 10/11 data in the form of dynamic mass factors: 

a N (r) + a P (r) 

T r = — , 6.15) 

a N (r) 

which is a rearrangement of Equation (6.10), and compare to the MOG prediction of 
Equation (6.9) for the parametric values of a(r) and A(r) of Equation (6.12) and Equation 
(6.13), respectively, using the best-fit values for the parameters given in Equation (6.14). 
The behaviour of G{t)/Gn is closely constrained to unity over the inner planets until 
beyond the orbit of Saturn (r > 10 AU) where the deviation in Newton's constant increases 
to an asymptotic value of G^/Gn —> 1.001 over a distance of hundreds of AU. 

Although MOND is not expected to provide a viable solution to the Pioneer 10/11 
anomaly because of Equation (6.2), the variation in the dynamic mass factor, T(r), shown 
in Panel (a) of Figure 6.2, is consistent with the deep MOND linear relation of Equation 
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(3.20) with a best-fit MOND acceleration of 



a = (3.0 ±0.3) x 10" n cm s 



-2 



(6.16) 



provided the MOND interpolating function is so gentle that the onset of the deep MOND 
regime occurs at 



instead of r t , which is not consistent with MOND, and improbable to explain using any 
generalized theory involving a preferred frame, as in §3.2.3, including Bekenstein's TEVES 
theory of §3.2.3. 

Since the density of the Milky Way dark matter in the vicinity of the solar system, 
according to Equation (6.3), is at least 10 orders of magnitude too small to affect the 
acceleration of spacecraft, then the Solar System must have its own halo for dark matter 
to provide a viable solution to the Pioneer 10/11 anomaly. Frere et al (2008) calculated 
the bound on the dark matter density of a spherical halo centered about the sun from 
high precision Solar System measurements, finding that a dark matter halo around the 
Solar System may be as much as 5 to 6 orders of magnitude more dense than the Milky 
Way's dark matter halo, but this is still at least 4 orders of magnitude too low to affect 
the acceleration of spacecraft in the Solar System. 

6.1.3 Kepler's laws of motion 

A consequence of a variation of G and GM for the solar system is a modification of 
Kepler's third law: 



where T PL is the planetary sidereal orbital period and clpl is the physically measured semi- 
major axis of the planetary orbit. For given values of clpl and T PL) Equation (6.18) can 
be used to determine G(r)M e . The standard method is to use astrometric data to define 
GM Q for a constant value, 



r = 10.5 ± 0.5 AU 



(6.17) 




(6.18) 



G(r)M = G(a e )M = k\ 



(6.19) 
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Figure 6.3: Solar system — Kepler's third law. 



The cube of the orbital distance, a 3 PL in AU 3 , vs. the square of the orbital period, 
Tp L in day 2 for the Solar system. The orbital data consist of the measured semi-major 
axis of the planateary orbit, clpl and the measured planetary sidereal orbital period, 
Tp L , listed in Table 6.2, respectively. Corrections due to MOG in Kepler's third law to 
Equation (6.8) is plotted using the result for T(r) of Figure 6.2. The figure is continued. 

where a e is the semi-major axis for Earth's orbit about the Sun, and k is the Gaussian 
gravitational constant given by 1 



For several planets such as Mercury, Venus, Mars and Jupiter there are planetary ranging 
data, spacecraft tracking data and radiotechnical flyby observations available, and it is 
possible to measure apL directly. For a distance varying GM e we derive (Fischbach and 

1 http : / / ssd . jpl . nasa . gov/? const ant s 



k = 0.01720209895 AU 3/2 /day. 



(6.20) 



We obtain the standard semi-major axis value at 1 AU: 




(6.21) 
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Figure 6.3 continued: Solar system — Kepler's third law. 

The cube of the orbital distance, o? PL in AU 3 , vs. the square of the orbital period, Tp L in day 2 
for the planets of the Solar system. The orbital data consist of the measured semi-major 
axis of the planateary orbit, clpl and the measured planetary sidereal orbital period, T PL) 
listed in Table 6.2, respectively. Corrections due to MOG in Kepler's third law to Equation 
(6.8) is plotted using the result for T(r) of Figure 6.2. 
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Talmadge, 1999; Talmadge et al, 1988): 



dpi 



1 + Vpl 



G(a PL )M^ 1/3 



K 2 



(6.22) 



Here, it is assumed that GM Q varies with distance such that tjpl can be treated as a 
constant for the orbit of a planet. We may substitute the Gaussian gravitational constant 
of Equation (6.19) into Equation (6.22) and obtain 



Vpl = 



GM _L (6 - 23) 



6.1.4 Planetary Ephemerides 

For the nine planets, we obtain the values of t]pl shown in Table 6.1. We see that we are 
able to obtain agreement well within the bounds of possible variation of GM consistent 
with the data (Fischbach and Talmadge, 1999; Talmadge et al, 1988) for Mercury, Venus, 
Mars and Jupiter. No observational limit on r]pp for Saturn or the outer planets has yet 
been established; but this is precisely where the deviation 5G(t)/Gn leads to a sizeable 
contribution in the theoretical prediction for r\pp. 

The reason for the uncertainty beyond the orbit of Saturn and the lack of observational 
limits on tjpl is that the ephemerides for the outer planets is based on optical measure- 
ments. Even in the context of Newton's theory, the extrapolation of Kepler's third law of 
Equation (6.18) using the Gaussian gravitational constant of Equation (6.20) which fits 
the inner planets missestimates the semi-major axis, ap^, or the orbital period, T PL , of the 
outer planets resulting from Newtonian perturbations due to Jupiter and the gas giants and 
their satellites, the Kuiper belt and hundreds of asteroids. The latest version of the plan- 
etary part of the numerical ephemerides is a numerical integration of the post-Newtonian 
metric. It attempts to account for these perturbations from Kepler's law beyond Saturn 
by a least squares adjustment to all the available observations including the CCD opti- 
cal astrometric observations of the outer planets. These values (without uncertainty) are 
available from the Solar System Dynamics Group (SSD) of the Jet Propulsion Laboratory 
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Table 6.1: Planetary predictions and observational limits 



Planet 


r 


Prediction 


Observational Limit 


ap 




(AU) 


Vpl (10- 10 ) 


Vpl (10- 10 ) 


(10" 8 cm/s 2 ) 


(i) 


(2) 


(3) 


(4) 


(5) 


Mercury 


0.38 


-6.55 x 10" 5 


+40 ± 50 


1.41 x 10" 10 


Venus 


0.72 


-6.44 x 10" 5 


-55 ± 35 


5.82 x 10" 8 


Earth 


1.00 


0.00 x 10° 





1.16 x 10" 6 


Mars 


1.52 


4.93 x 10" 3 


-0.9 ±2.1 


4.42 x 10" 5 


Jupiter 


5.20 


4.19 x 10 2 


+200 ± 400 


2.76 x 10" 1 


Saturn 


9.54 


1.67 x 10 4 




3.27 x 10° 


Uranus 


19.22 


1.84 x 10 5 




8.86 x 10° 


Neptune 


30.06 


4.39 x 10 5 




8.65 x 10° 


Pluto 


39.52 


6.77 x 10 5 




7.72 x 10° 



Notes. — Theoretical predictions of the values of t]pl and the best-fit theoretical predictions for 
the Pioneer Anomaly, ap, for the planetary bodies and observational limits. Planetary bodies are 
listed in Column (1), with their mean distance, r, from the Sun shown in Column (2). Column (3) 
is the theoretical prediction oiriPL of Equation (6.23), and may be compared to Column (4) for the 
observational limits taken from Talmadge et al (1988). No observational limits were computed 
beyond Saturn in Talmadge et al. (1988) due to uncertainty in opical data. Beyond the outer 
planets, the theoretical predictions for 77 (r) approaches the asymptotic value 7700 = 3.34 x 10 -4 . 
Column (5) lists the anomalous accelerations, at the planetary positions, predicted by the best-fit 
to the Pioneer 10/11 anomaly. 

(JPL) through the Horizon's ephermeris DE410 online 2 . The Russian Academy of Sciences 
has also placed their latest values known as EPM2004 online 3 . Because the perturbations 
change daily due to the motion within the solar system, the planetary ephemerides quoted 
values for clpl and T PL change daily. In order to compute deviations from Kepler's third 
law for the outer planets, we have listed today's best known values in Table 6.2. 

The uncertainty in the EPM2004 deduced values for the semi-major axes of the planets, 
AapL, have been studied in Pitjeva (2005) and the quoted values are listed in Table 6.2. 
Pitjeva (2005) warned that the real errors may be larger by an order of magnitude. The 
uncertainty in the periods for the outer planets are not quoted in either EPM2004 or 
DE410, and so we have assumed small uncertainties based on the precision provided by 
the JPL Horizon's online ephemeris. 



2 http: / / ssd.jpl.nasa.gov/horizons.html 

3 ft p : / / quasar . ipa. nw . ru / incoming / EPM2004 
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We may calculate the uncertainty, Ar]p L) by propagating the errors AapL and AT PL 
according to Equations Equation (6.21) and Equation (6.22), neglecting any uncertainty 
in the Gaussian gravitational constant of Equation (6.20): 



Although according to Table 6.1 we are consistent with the observational limits of t]pl 
for the inner planets to Jupiter, the computation of Talmadge et al (1988) attempted to 
set model-independent constraints on the possible modifications of Newtonian gravity. The 
procedure was to run the planetary ephemerides numerical integration with the addition 
of t]pl as free parameters. Because there was one additional parameter for each planet, 
they were only able to find observational limits for the inner planets including Jupiter. 
In order to compute the observational limit for tjpl for the outer planets, it would be 
necessary to compute the planetary ephemerides using the modified acceleration law of 
Equations Equation (6.4) and Equation (6.5). Although this is beyond the scope of the 
current investigation, we may approximate here the observational limit of t]pl for the outer 
planets as the uncertainty Atjpl from Equation (6.24), for the perturbations of Figure 6.2, 
SG(r) /Gat, are small compared to the Newtonian perturbations acting on the outer planets. 
The results for Atjpl due to the uncertainty in the planetary ephemerides are presented 
in Table 6.2 for the nine planets and exceed the predictions, t]pl, of Table 6.1. 

Modified gravity can explain the Pioneer anomalous acceleration data and still be 
consistent with the accurate equivalence principle, lunar laser ranging and satellite data 
for the inner solar system as well as the outer solar system planets including Pluto at a 
distance of r = 39.52 AU = 5.91 x 10 12 meters. The ephemerides for the outer planets are 
not as well know as the inner planets due to their large distances from the Sun. 

The orbital data for Pluto only correspond to the planetoid having gone round 1/3 of 
its orbit. It is important that the distance range parameter lies in the region 47 AU < 
A(r) < oc for the best-fit to the Pioneer acceleration data, for the range in the modified 
Yukawa correction to Newtonian gravity lies in a distance range beyond Pluto. Further 
investigation of fifth force bounds obtained by an analysis of the planetary data for the 




(6.24) 
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outer planets, based on the modified gravity theory is required. We are predicting that 
measurements of a fifth force in the solar system will become measurable at distances 
r > 10 AU from the Sun where as shown in Figure 6.2, 5G(r)/Go (and tj P l) become 
potentially measurable. The likely possibility that the Pioneer 10/11 anomaly is caused 
by thermal effects would cause these predictions to be treated as bounds on the effects of 
MOG in the solar system. Moffat and Toth (2009) give good agreement with solar system 
bounds using the scalar-tensor-vector modified acceleration law of §3.4.4. 

6.1.5 Anomalous perihelion advance 

The relativistic equation of motion for a test particle in our gravitational theory may be 
solved perturbatively in a weak field approximation for the anomalous perihelion advance 
of a planetary orbit: 

6ttG n M q 

ACJ PL = — 2~-(l ~ a PL), ( 6 - 25 ) 



where we have assumed as with Kepler's third law that GM e and a vary with distance 
such that they can be treated as constants for the orbit of a planet, where we have made 
use of the approximation G(r) ~ Gn Moffat (2006b), which is the case from the fit to 
the Pioneer 10/11 anomalous acceleration data. We may rewrite Equation (6.25) as the 
perihelion advance in arcseconds per century: 

A^pl 3G N M e 
2ttT pl c 2 a PL (l - e z PL )T PL 

where T PL is the planetary orbital period, and e P L is the planetary orbital eccentricity. 
We may separate Equation (6.26) into the usual Einstein anomalous perihelion advance, 
and a prediction of the correction to the anomalous perihelion advance: 

vpl =vo+vu (6.27) 

where 

3G N M e 

= ~2 T\ 2 \rr > 6 - 28 

c 2 a PL (l - e 2 PL )T PL 

Cji = -a PL u) Q) (6.29) 
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are the Einstein anomalous perihelion advance, and the predicted retrograde, respectively. 
Note the minus-sign in the predicted retrograde of Equation (6.29) as compared to the 
Einstein anomalous perihelion advance of Equation (6.28). The measured perihelion pre- 
cession is best known for the inner planets (for Mercury the precession obtained from 
ranging data is known to 0.5% (Will, 2006)). For each of the planets in the solar system, 
we find that apL « 1, so that our fit to the Pioneer anomalous acceleration is in agree- 
ment with the relativistic precession data. The results for the Einstein perihelion advance, 
and our predicted retrograde for each planet, and the observational limits set by the recent 
ultra-high precision ephemeris are listed in Table 6.3. 

The validity of the bounds on a possible fifth force obtained from the ephemerides of the 
outer planets Uranus, Neptune and Pluto are critical in the exclusion of a parameter space 
for our fits to the Pioneer anomaly acceleration. Beyond the outer planets, the theoretical 
prediction for r)(r) approaches an asymptotic value: 

r]^ = lim r](r) = 3.34 x 10~ 4 . (6.30) 

We see that the variations (running) of a(r) and A(r) with distance play an important 
role in interpreting the data for the fifth force bounds. This is in contrast to the standard 
non-modified Yukawa correction to the Newtonian force law with fixed universal values of 
a and A and for the range of values < A < oc, for which the equivalence principle and 
lunar laser ranging and radar ranging data to planetary probes exclude the possibility of 
a gravitational and fifth force explanation for the Pioneer anomaly. 

Perhaps, a future deep space probe can produce data that can check the predictions ob- 
tained for the Pioneer anomaly from modified gravity theory. Or perhaps utilizing Mars or 
Jupiter may clarify whether the Pioneer anomaly is caused by the gravitational field. (Page 
et a/., 2006). An analysis of anomalous acceleration data obtained from earlier Doppler 
shift data retrieval will clarify in better detail the apparent onset of the anomalous accel- 
eration, or support the thermal recoil explanation of Toth and Turyshev (2009), perhaps 
to as low a x 2 as the modified gravity solution. 



Part III 



Conclusions 



Chapter 7 



Summary 

"Any intelligent fool can make things bigger and more complex and more vio- 
lent. It takes a touch of genius - and a lot of courage to move in the opposite 
direction. " 

Albert Einstein 



The mysteries of the gravitational field continue to challenge mankind as our physi- 
cal models of the universe evolve. Isaac Newton's great contribution was to deduce the 
analytical form of the force of gravity exerted by an isolated object. However, the prin- 
ciple of superposition is not exact, as Newton assumed, because the gravitational field is 
non-linear, and models of the gravity internal to astrophysical matter distributions are 
Newtonian approximations. Albert Einstein's great contribution was to deduce the geo- 
metric form of the relativity principle, reinterpreting the force of gravity as a geometric 
distortion of space and time, but this sets the gravitational field apart from the three other 
known forces. Gravity is the only force that couples universally to matter and energy. 
However, the strong equivalence principle, which Einstein first assumed, does not hold for 
stable gravity theories that include scalar, vector and tensor modifications to the met- 
ric with associated couplings, even when additional charged quantum numbers associated 
with new symmetries are suppressed. Such modified gravity theories suggest the presence 
of a fifth force, which is assumed to couple universally to matter and energy, and gains 
in strength at astrophysical scales to become the dominant force. This dominant force, if 
neglected by means of a Newton-Einstein approximation, emerges as the phantom of dark 
matter in galaxies and clusters of galaxies. 

To address the hypothesis, stated in §1.1.3, to the missing mass problem in galaxy 
rotation curves and clusters of galaxies, the following theories were studied: 

1. Cold non-baryonic dark matter (CDM), 

2. Milgrom's modified Newtonian dynamics (MOND), 
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3. Moffat's metric skew-tensor gravity theory (MSTG), 

4. Moffat's scalar tensor vector gravity theory (STVG). 

Conclusions drawn upon identification of the missing mass as CDM is presented in §7.1.1. 
Some common lessons learned from the modified gravity theories are summarized in §7.1.2. 
Corresponding to each of the astrophysical scales in Table 1.2, conclusions based upon 
galactic-scale and cluster-scale astrophysics are summarized in §7.2 and §7.3, respectively. 
Some possible directions for future astrophysical tests are suggested in Chapter 8. 

7.1 Lessons learned 

Whether identified as dark matter halos with profiles defined in Chapter 2, or the massive 
fifth-force fields of a modified gravity theory with acceleration laws derived in Chapter 3, 
conclusions are drawn in §7.1.1 and §7.1.2, respectively. 

7.1.1 CDM halos 

According to Table 4.2 from the Ursa Major sample of galaxies of §4.2, every one of the 
galaxy rotation curves presented in Figure 4.1 have excellent fits, within Einstein-Newton 
gravity including a non-baryonic dark matter halo density described by the core-modified 
dark matter profile of Equation (4.2), with reasonable mean stellar mass to light ratios T, 
as compared to baryon suppressed fits with T = 0. 

Provided baryons are included in the core-modified dark matter halo, the total mass 
with dark matter vs. velocity relation, plotted in Figure 4.8, showed the least scatter of 
any of the Tully-Fisher relations, even when compared to the modified gravity alternatives 

Shown in the same table, in the leftmost columns, the Navarro-Frenk- White (NFW) 
profile of Equation (4.1) showed low x 2 when including the visible baryons, but because 
the cusp problem occurs precisely where baryons are important, there were exceptions in 
the sample of Chapter 4, involving high and low surface dwarf galaxies, that could not be 
X 2 best-fit for any non-zero value of T. For those galaxies that the NFW formula produced 
excellent best-fits, the halo mass function overcompensated for the baryons by suppressing 
the best-fit stellar mass to light ratio, T < 1. 

The dynamic mass factors show no indication of a cusp at small r, but unlike for 
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clusters of galaxies, take on a global minimum, with lim r ^ r(r) = 1, whereas the core- 
modified dark matter formula provides excellent fits to the dynamic mass factors, T(r) for 
all positions, r, plotted in Figure 4.4, for every galaxy in the sample including the dwarfs. 

In the CDM hierarchical structure formation scenario, galactic halos are considered 
subhalos to the larger structure which is their cluster (or filament). These subhalos are 
self-similar to the cluster halo if they are describable by similar fitting formulae. Shown 
in Panel (a) of Table 5.2, the core-modified dark matter profile of Equation (4.2) showed 
the lowest x 2 best-fit cluster model parameters for the sample of Chapter 5, including 
X-ray clusters of varying mass, scale radius, and central temperature, as compared to the 
modified gravity alternatives. 

The NFW profile could not be x 2 best-fit to the sample of X-ray clusters because the 
data does not exhibit the cusp in the core, as shown by the variation in the dynamic mass 
factor, plotted in Figure 5.2. 

Therefore, the core-modified dark matter formula describes a self-similar halo profile for 
both high and low surface brightness galaxies, including the dwarfs, and also for the clusters 
of galaxies, including the dwarf clusters and the Bullet Cluster 1E0657-558, provided the 
baryons are not suppressed in the fits. 

7.1.2 Modified gravity theories 

Unlike the core-modified dark matter or NFW profiles, each of the modified gravity theories 
produced low x 2 fits to the galaxy rotation curves of Chapter 4 with universal parame- 
ters, averaged over the sample, shown in Table 4.3 for MOND, Table 4.4 for MSTG, and 
Table 4.5 for STVG. The monotonic (near-linear) rise of the dynamic mass factor, plotted 
in Figure 4.4 across all gravity theories, suggests that the missing mass problem is most 
pronounced at the edge of the luminous disk. 

Only the STVG acceleration law of §3.4.4 is derived without any phenomenological 
input from the Tully-Fisher relation. The Tully-Fisher relations are compared in Table 4.7, 
for each of the gravity theories, and plotted in Figure 4.9. 

For the sample of X-ray clusters of galaxies of Chapter 5, it was confirmed that MOND 
provides poor best-fits to X-ray cluster mass profiles, as plotted in Figure 5.1, even with a 
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varying MOND acceleration, <2q, as shown in Panel (b) of Table 5.2. The possibility, within 
MOND, of including a non-luminous component, such as neutrino halos, was considered in 
the case of the Bullet Cluster 1E0657-558 in §5.3.5, although the models of Sanders (2003, 
2007) may be overestimating the extent of the neutrino halos by a factor of two to provide 
a universally consistent explanation. 

Unlike MOND, both MSTG and STVG theories provided excellent best fits to the 
X-ray cluster mass profiles of Figure 5.1, as shown in Panels (c) and (d) of Table 5.2, 
respectively. Whereas the STVG theory produced a constant dynamic mass factor, which 
is an approximate feature present in the data, the MSTG produced dynamic mass factors 
which mimicked the core-modified dark matter result and produced excellent fits, as shown 
in Figure 5.2, from the smallest of the clusters to the largest Bullet Cluster 1E0657-558. 

Effectively, the common feature of the modified gravity paradigm is given by Equation 
(3.93), in which 

M(r) = GM^ {71) 

where G(r) is the spatially varying gravitational coupling. The modified gravity hypothesis 
of phantom dark matter suggests that, for sufficiently large r, G(t)/Gn 3> 1 so that the 
luminous, baryonic mass, M(r), is less than the observed Newtonian dynamic mass, M N {r), 
by the same factor. 

The metric skew-tensor gravity (MSTG) theory, presented in §3.3, identifies the phan- 
tom of dark matter with the massive skewon of the Kalb-Ramond-Proca field, as described 
in §3.3.2, with a Yukawa interaction that leads to motion under the fifth force, as in §3.3.3. 
For MSTG astrophysical predictions, the gravitational coupling, G(r), of Equation (3.91), 
was substituted into Equation (7.1) leading to a nonlinear equation because G(r) itself is 
a function of M(r), which was solved exactly in Equations (3.95) and (3.96). This solu- 
tion permits the analytic computation of baryon masses according to MSTG from dynamic 
measurements such as velocity rotation curves or X-ray temperature distributions, without 
approximation. 

Using a simpler field structure, the MSTG skewon can be replaced by a massive vec- 
tor phion in the scalar-tensor-vector gravity (STVG) theory, presented in §3.4, which then 
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identifies the Maxwell-Proca field , as in §3.4.2, as the phantom of dark matter. The STVG 
action includes three dynamical scalar fields which lead to gravitationally strong interac- 
tions, as shown in the Yukawa phion theory of §3.4.4, which leads to the modified central 
force law of §3.4.5. For STVG astrophysical predictions, the gravitational coupling, G(r), 
of Equation (3.143), was substituted into Equation (3.93) leading to Equation (3.164), 
where the Yukawa strength and phion mass, a and /x, are given by Equations (3.159) and 
(3.160), respectively. 

To compare the predictions of the gravity theories relevant to astrophysical scales, 
Equation (7.1) can be written in terms of the dynamic mass factor, of Equation (2.22), 

mr) = (7.2) 

It was shown that MOND may be the weak-field limit of certain Lorentz violating theories, 
including the family of theories discussed in §3.2.3, which have non-metric field structure 
that are subject to violations of the strong equivalence principle, as described in §3.1.3, 
and may thereby provide the phantom of dark matter detected in galaxy rotation curves 
and X-ray cluster masses. In §3.2.2 it was shown that the dynamic mass factor due to mod- 
ified dynamics at small accelerations, is precisely the inverse of the MOND interpolating 
function, according to Equation (3.15), 

r(r) = l/fi(x(r)). (7.3) 

At sufficiently small accelerations, the dynamic mass factor is a linear function in r ac- 
cording to Equation (3.20), for all MOND interpolating functions, and is proportional to 
y/ao. Whereas this fits the observations in galaxy rotation curves, it does not correspond 
to the observations in X-ray cluster masses, even with larger values of the MOND universal 
acceleration. 

The dynamic mass factor, 

is the ratio of the Newtonian dynamic mass to the observed, baryonic mass. For CDM, 
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Equation (7.4) is constrained by the particular choice of the dark matter fitting formula. 
7.2 Galactic astrophysics 

In Chapter 4, a core-modified fitting formula was derived in §4.1: Core-modified profile, 
and found to fit the sample of high and low surface brightness galaxies, including all of 
the dwarfs. The NFW fitting formula led either to fits with very poor x 2 ) if baryons 
were omitted, or to a suppression of the stellar mass-to-light ratio if photometric data was 
included, particularly for the dwarf galaxies. The worst of these dwarf galaxies could not 
be fitted using the NFW profile with any non-zero value of the stellar mass-to-light ratio, 
as discussed in §4.1.1: Observations. 

The core-modified fitting formula of Equation (2.12) produced the lowest reduced-;^ 2 
best-fits to the galaxy rotation curves for the sample of §4.2, plotted in Figure 4.1, with two 
parameters, p and r s , which varied across the sample, tabulated in Table 4.2. The surface 
mass densities, plotted in Figure 4.2, show the baryon dominated cores transitioning to 
dark matter dominated halos. The best-fitting mass profiles are plotted in Figure 4.3, 
showing the halo component negligible in the galaxy cores, but adding up to the dominant 
mass at the outermost radial points. At large distances from the center of each galaxy in 
the sample, the density profile of the dark matter halo is well described by a steep power- 
law, with power-law index 7^3, whereas at distances toward the center of the galaxy an 
increasingly shallow power-law is observed, as plotted in Figure 4.5. For distances less than 
the dark matter halo core radius, r < r s ,the total density profile including baryons shows 
a universal 7 — ► 1 power-law index, and the density profile of the dark matter component 
alone approaches a rarified, constant density core. 

§4.1.2: Observations discusses the MOND best-fits to the galaxy rotation curves for 
the sample of §4.2, plotted in Figure 4.1, with universal acceleration, a , tabulated in 
Table 4.3. 

§4.1.3: Observations discusses the MSTG and STVG best-fits to the galaxy rotation 
curves for the sample of §4.2, plotted in Figure 4.1, with universal parameters depending 
on the theory, tabulated in Table 4.4 and Table 4.5, for MSTG and STVG respectively. 

All of best-fits showed a central disk dominated by the Newtonian potential, where 
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r(r) ~ 1, outside of which, the dynamical mass factor increased approximately linearly 
with distance, plotted in Figure 4.4. The results, plotted (with lines) are nearly equivalent 
for all of the gravity theories studied. However each theory's data (with error bars) has a 
dependence on the result of the best-fit stellar mass-to-light ratio, T, provided in Table 4.2 
for the NFW profile and the core-modified dark matter profile, Table 4.3 for MOND, 
Table 4.4 for MSTG, and Table 4.5 for STVG. 

The analysis of the Ursa Major sample of §4.2 involved a series of calculations using a 
variety of computational resources. The error analysis was a fruitful exercise in measuring 
the properties of the variation of the stellar mass-to-light ratio within high or low surface 
brightness galaxies, and the quality of the best-fit. The stellar mass-to-light ratio, T(r), 
varied strongly within every galaxy, in comparison to the variation across different galaxies 
of the best-fit stellar mass-to-light ratio. 

The Newtonian core was calculated from the plot of the stellar mass-to-light ratio, T(r) 
of Figure 4.6 (vertical lines), at the position, r c , where the Newtonian dynamics induce 
a rapid increase in the slope of T(r), for every galaxy including the dwarfs. The best-fit 
stellar mass-to-light ratio for the Newtonian core is plotted (horizontal line) to r = r c . 
Within the luminous disk, the stellar mass-to-light ratio never exceeds a value of ten in 
any galaxy. This naturally constrains the total amount of dark matter required and allows 
the best-fit dark matter theory, with two parameters p and r s , and the best-fit stellar 
mass-to-light ratio, T, to be simultaneously varied toward minimum x 2 . 

The stellar mass-to-light ratios, T(r), plotted in Figure 4.6 (with lines) are nearly 
equivalent for MOND, MSTG and STVG, with similar mean values, near T = 1, with 
larger values for the best-fit Newtonian core model, in all galaxies. The best-fit values, 
determined by a non-linear least squares fitting algorithm, are plotted (with horizontal 
lines). Since the galaxies of Ursa Major are at a common distance from the Milky Way, 
the variation in the actual stellar mass-to-light ratio from galaxy to galaxy is not expected 
to be large. The results for the best-fit stellar mass-to-light ratio, with a core-modified 
dark matter halo, are plotted (with horizontal lines), with mean values close to predicted 
values according to MOND and the MOG theories. This is an example of the importance 
of the luminous baryons in the computation. It is the variable stellar mass-to-light ratio 
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which includes the data, by force, and allows the ultra-low reduced x 2 test for the core- 
modified dark matter model. The computational results are provided in Table 4.2 for the 
NFW profile and the core-modified dark matter profile, Table 4.3 for MOND, Table 4.4 
for MSTG, and Table 4.5 for STVG. 

Qualitative assessment of each theory's predictions for the Ursa Major galaxies is 
provided in §4.3, and includes the predictions of the best-fit Newtonian core model of 
§4.3.1, which is a base-line for any improvement. Whereas repatriating orphan features, 
as described in §4.3.2, provides a reasonable test for theories which fit galaxy rotation 
curves, particularly the core-modified dark matter halos which properly include the lumi- 
nous baryons, as in §4.3.3. The power-law profile, for either the NFW or core-modified 
models is derived as the logarithm slope of Equation (4.52), which relates 

which depends on the baryon distribution through the mass-to-light ratio, according to 
Equations (4.38) and (4.40). The halo component is computed by substituting M = Mh a io 
into Equation (7.5) and the dark matter logarithm slope is computed by substituting 
M = |Mhi + TMdi S k + Mh a io, thereby including the gaseous and luminous stellar disks 
into the Newtonian dynamic mass as discussed in §4.2.7: Solution to the dark matter cusp 
problem. 

The implication of Occam's razor, that the total mass of a galaxy should be less in a 
theory without non-baryonic dark matter, depends on how large the halo is taken to be, 
beyond the luminous disk, where there is data. The results for the variation in the total 
integrated mass, M(r), to the outermost radial position, r out , are plotted in Figure 4.3 
(with lines), per theory. Components are plotted for the actual HI gas, and the stellar disk 
is normalized with T = 1 for relevance across theories, each with a best-fit T. The dark 
matter halo component is plotted, which is a small part of the total mass in the dark matter 
model in the core of every galaxy studied. However the dark matter halo component is the 
fastest rising mass in the galaxy because of the spherical distribution, compared to both 
the exponentially-thin gaseous disk, and the best-fit luminous disk, becoming dominant 
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outside the Newtonian core, r c , per galaxy. 

The Tully-Fisher relation, as in §4.3.4, confronts the dynamical importance of the 
luminous baryons compared to the ordinary Tully-Fisher relation, plotted in Figure 4.7 
(blue) with power-law index of a — 4.1 ± 0.4 according to the best-fit of Equation (4.59). 
This result, although familiar, is today considered too large with larger samples providing 
a K-band Tully Fisher relation power-law index of a = 3.4 ±0.1 (Tully and Pierce, 2000). 
Table 4.7 provides the best-fit logarithm slopes and intercepts for six relations of the 
form M oc v a ) where either M = TMdi S k for the stellar Tully-Fisher relation, or M = 
|M H i + TMdi S k for the baryonic Tully-Fisher relation, or M = |M H i + TMdi S k + M^aio for 
the total mass-velocity relation, including the best-fit core-modified dark matter halo, as 
in Figure 4.8. The best-fits (with lines) to the theoretical Tully-Fisher relations are plotted 
(with error bars) for the two stellar and two baryonic relations in Figure 4.9, per gravity 
theory. 

In MOND, the relation with the least scatter is the baryonic mass to outermost velocity, 
v out , with a logarithm slope ofa = 3.3±0.3 whereas the ordinary, ^ max , stellar relation has 
a logarithm slope of a = 4.1 ± 0.4. However, the STVG, MSTG and core-modified dark 
matter theories show less scatter for the ordinary baryonic relation with smaller logarithm 
slopes of a = 2.6 ± 0.2, a = 2.5 ± 0.2, and a = 2.5 ± 0.3, respectively. Whereas overall, the 
relation which shows the least scatter is that of the total mass, including luminous baryons 
and core-modified dark matter, vs. velocity as plotted in Figure 4.8, with a = 2.9 ± 0.2 
for the ordinary case. This restores the Tully-Fisher relation to the dark matter solution, 
provided the baryons are included, and dismisses the notion that the Tully-Fisher relation 
is unnatural due to dark matter dominance. 

Some possible directions for future galactic astrophysical tests are presented in §8.1. 

7.3 Cluster-scale astrophysics 

In Chapter 5, the missing problem is studied with X-ray clusters, as described in §5.1, 
using the astrophysical sample of §5.1.1. The King /3-model of the X-ray gas distribution 
of Equation (5.2) is presented in §5.1.2, and the collisionless Boltzmann equations are 
derived, and the solution is shown in Equation (5.15). In §5.1.3, a computation expresses 
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the surface mass density map, E(x,y), as the simple analytical result of Equations (5.24) 
and (5.25), in terms of the King /?-model best-fit values (p , A r c ). This was used in the 
analysis of the Bullet Cluster 1E0657-558, in §5.3, as the initial study of the X-ray gas 
map, described in §5.3.1. The Newtonian dynamical mass of Equation (5.28) is derived in 
§5.1.4. 

Each of the theories that were tested at galactic-scale using galaxy rotation curves are 
studied at cluster-scale using the best-fit cluster models of §5.2. The NFW fitting formula 
of Equation (2.7) generated large uncertainties due to a parameter degeneracy between 
the central density parameter, p 0j an d the scale radius, r 5 , and could not be x 2_ fitted. 
Without numerical convergence, the NFW results either over-predicted the density at the 
core or under-predicted the total mass. Alternatively, the core-modified dark matter model 
of §5.2.1 provided excellent fits using the fitting formula of Equation (2.12), with results 
provided in Panel (a) of Table 5.2. 

Although MOND does not fit X-ray cluster masses, there have been studies that claim 
improvements using a larger value of the MOND universal acceleration, ag, or to include 
a non-luminous component. Both avenues were considered in this thesis. In §5.2.2, the 
MOND universal acceleration was treated as a variable parameter, and found to lead to 
very poor fits with larger than galaxy-scale values of ag, as shown in Panel (b) of Table 5.2. 
Furthermore, the problem of too large a dynamic mass factor, for small r, and too small a 
dynamic factor for intermediate r, was not corrected although the increase in the value of 
a did lead to the correct dynamic mass factor for r ~ r out . This means that the best-fit 
MOND solution in clusters of galaxies without dark matter does not fit the shape of the 
X-ray mass profile, except at r out so that the total mass is corrected. 

Both MSTG and STVG theories, as discussed in §5.2.3, provided excellent fits to the 
X-ray gas masses, with results provided in Panels (c) and (d), respectively, of Table 5.2. 

The missing mass problem at cluster-scale, presented in §5.2.4, is best demonstrated 
by the dynamic mass factor of Equation (7.4), which is plotted (with lines) in Figure 5.2, 
per gravity theory. The observations are plotted (with error bars) as the ratio of the 
Newtonian dynamical mass, M N {r) of Equation (5.28), to the best-fit King /3-model to 
the gas mass, M(r), of Equations (5.2) and (5.3), and is therefore a ratio of two very 
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large masses, increasing with separation, r, as plotted (with red crosses, and green circles, 
including error bars). 

The study of the Bullet Cluster 1E0657-558, in §5.3, includes a detailed analysis of the 
X-ray gas map, of §5.3.1, with the subcluster masked-out, and a best-fit to the King (3- 
model for the main cluster is derived in Equations (5.41), (5.42) and (5.44). The subcluster 
subtracted from the X-ray data was added to the best-fit King (3 model of the main cluster, 
as shown in Figure 5.6. The study of the gravitational lensing convergence map, according 
to §5.3.2, applies the derivation of Equations (4.29) and (4.30) with running gravitational 
coupling, G(r), leading to the best-fit ft-model of Figure 5.8. The best-fit MOG model 
is used to compute the visible baryon surface mass profile in §5.3.3, and the result of 
the galaxy subtraction of Equation (5.52) is shown in Figure 5.9. Predictions for the 
component masses of the main cluster, subcluster and central ICM are listed in Table 5.5. 
The dark matter distribution is computed as the difference between the ft-map and the 
scaled S-map by a pixel by pixel subtraction in §5.3.4. The distribution of visible and dark 
matter, plotted in Figure 5.10, provides a comparison of the distribution of galaxies, gas 
and total baryons according to MOG in Panel (a), as compared to the distribution of gas 
and dark matter in Panel (b). Although the NFW profile does not fit the main cluster, 
the excellent fit to the core-modified profile is shown in Panel (a) of Figure 5.1. 

Some possible directions for future cluster-scale astrophysical tests are presented in 
§8.2. 



Chapter 8 



Future astrophysical tests 

"If I have seen further than others, it is by standing upon the shoulders of 
giants. " 

Sir Isaac Newton 

The work of Part II explored the importance of directly measuring the dynamics of 
the gravitational field and comparing to the observed galactic and gaseous components, 
within the gravity theories of Part I. The results show a strong interplay between current 
astrophysical measurements and the resulting surface density maps predicted by gravity 
theory, with predictions which are testable and falsifiable. Some possible directions in 
future space observations are presented in §8.1 for galactic-scale astrophysics, and in §8.2 
for cluster-scale astrophysics. 

8.1 Galactic astrophysics 

Testing whether galaxy dynamics are dominated by a distribution of cold non-baryonic dark 
matter, like the models of Chapter 2, or whether galaxy dynamics are dominated by a 
modified gravity theory which violates the strong equivalence principle, such as any of the 
candidates of Chapter 3, is possible through a combination of sub-kiloparsec resolution 
luminous disk observations, galaxy rotation curve measurements, and future directions in 
galaxy-galaxy lensing, depending on the next generation of space observatories. Future 
directions for rotation curve methods, and galaxy-galaxy lensing are explored in §8.1.1, 
and §8.1.2, respectively. 

8.1.1 Galaxy rotation curves 

The rotation curves for the Ursa Major filament of galaxies, studied in §4.2, are already 
sufficiently detailed to provide features challenging to any candidate theory's best-fit. As 
the state of the art of computational models of gravity theories continues to improve, 
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astrophysical observations of nearby filaments and clusters of the Local Group offer the 
best possible future laboratory. 

Within the Local Group, the Milky Way is the most well studied galaxy, but the 
available rotation curve remains poorly known. Brownstein and Moffat (2006a) performed 
a MOND and MSTG best-fit to the Milky Way galaxy rotation curve supplied by Sofue 
(1996) using a parametric model for the surface mass density, independent of photometric 
observations, whereas the method using photometry, as described in §4.2.1, for the surface 
mass computation of §4.2.4, and best-fitting by means of the stellar mass to light ratio, 
as described in §4.2.8, is preferred. According to our best-fit parametric model, the total 
baryonic mass of the Milky Way was determined to be 



A new study of the rotation curve for the Milky Way should include better high- 
resolution rotation curve measurements, particularly for orbital distances greater than the 
Solar system's, and an updated model of the luminous stellar and gaseous components using 
available photometry, providing a better answer of the mass of our galaxy as determined 
by each gravity theory. The 21-cm line emission of neutral hydrogen, as traced by the 
Leiden/ Argentine/Bonn (LAB) galactic HI survey (Kalberla et al, 2005), is a full sky map 
which may be used to provide the gas component in Equation (4.34), assuming a big-bang 
nucleosynthesis relation such as Equation (4.23). Spitzer Space Telescope's encompassing 
infrared view of the plane of the Milky Way consists of nearly one million images which 
have been composed into a 120° mosaic, and may be used to compute the vertical scale 
height, z , and provide the stellar disk component in Equation (4.35), assuming a Van der 
Kruit and Searle law such as Equation (4.24). 

The Milky Way has three large dwarf satellite galaxies. Sagittarius is a 6 kpc dwarf 
spheroidal, located just below the galactic plane, twice as close as the Large and Small 
Magellanic Clouds, which are between 4 and 8 kpc in diameter. In addition, the Milky Way 
has sixteen dwarf spheroidal galaxies each less than 2 kpc in diameter. These satellite dwarf 
galaxies are an important part of the future laboratory for rotation curve astrophysics. 
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8.1.2 Galaxy- Galaxy lensing 

Although galaxy rotation curves, as in Figure 4.1, are the primary research tool for mea- 
suring dark matter or predicting modified gravity line-of-sight surface mass densities, as 
in Figure 4.2, and their galactic mass profiles, as in Figure 4.3, it is the dynamic mass fac- 
tor of Figure 4.4 which best shows where the missing mass according to Einstein-Newton 
gravity is distributed. For every high and low surface brightness galaxy, including the 
dwarf galaxies, the dynamic mass factor increases monotonically with distance from the 
galactic center. It is this factor which amplifies the different predictions between cold non- 
baryonic dark matter and each of the modified gravity solutions. E-maps, as in Figure 4.2 
for the Ursa Major filament of galaxies, are predictions for future galaxy-galaxy lensing 
measurements. 

Galaxy-galaxy lensing is a unique method to probe the dynamic mass of the foreground 
(lens) galaxy out to large radii. Strong galaxy-galaxy lensing, in which multiple background 
(source) galaxies appear in the image, allow direct measurements of the dynamic mass on 
scales < 10 kpc. The Center for Astrophysics at Harvard - Arizona Space Telescope Lens 
Survey (CASTLES) has identified (~ 100) strong gravitational lens galaxies from Hubble 
Space Telescope images. However, in order to directly measure the line-of-sight surface 
mass distribution across the foreground galaxy out to ~ 100 kpc, precision weak lensing 
is used, whereby the slight, coherent, gravitational lensing induced shape correlations can 
be averaged over multiple (~ 1000) source galaxies. The Cosmological Evolution Survey 
(COSMOS) has observed a 2 square degree equatorial field with the Advanced Camera 
for Surveys aboard the Hubble Space Telescope, and Leauthaud et al (2007) report on 
the source catalogue constructed from COSMOS, containing almost 400,000 galaxies with 
shape measurements and uncertainties. Refregier and Douspis (2008) report on plans for 
the Dark Universe Explorer (DUNE) wide-field mission concept, consisting of a 1.2 m tele- 
scope designed to carry out an all-sky survey in one visible and three near infrared bands, 
optimised for weak gravitational lensing. With the future direction of galaxy-galaxy lens- 
ing measurements, sub-kiloparsec distributions of the dynamic and visible components are 
attainable, providing high resolution images of the distribution of matter within galaxies. 

Whether identified as CDM or the massive fifth-force fields of a modified gravity theory, 
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the solutions to the missing mass problem are testable and falsifiable by current galactic 
observations, as in Chapter 4, and through future directions in measuring galaxy rotation 
curves, with infrared photometry on the stellar and gaseous disks, as in §8.1.1, and through 
techniques such as gravitational lensing for directly measuring the dynamic mass, as in 
§8.1.2. 

8.2 Cluster-scale astrophysics 

Some tests involving X-ray space observatories for measuring the intercluster medium are 
introduced in §8.2.1, which should be expected to generate higher resolution and greater 
bandwidth E-maps, which are correlated to the ionized electron density, and thereby a 
measurement of the visible baryon distribution. Some directions involving the Hubble and 
Spitzer space telescopes for the future of gravitational lensing measurements are introduced 
in §8.2.2, which provide T-maps correlated to the missing mass, whether identified as dark 
matter, or the massive fifth-force fields of a modified gravity theory, according to the 
lessons learned in §7.1. 

8.2.1 Intercluster medium 

The intercluster medium (ICM) is X-ray measurable with temperatures as high as 15.5 ± 
3.9 keV, from our 2007 value shown in Table 5.4 for the main cluster of the Bullet Cluster 
1E0657-558 (Brownstein and Moffat, 2006b), based on the Chandra observatory's very high 
resolution surface mass density S-map, shown in Figure 5.4. However, Chandra's band- 
width limitation prevents the Bullet Cluster 1E0657-558 being identified as the hottest 
cluster in the sky, until the next generation of X-ray space observatories, and a combina- 
tion of the high bandwidth, but low resolution, data available from methods involving the 
Sunyaev-Zeldovich effect. New S-maps of the ICM will therefore be of even greater preci- 
sion and resolution, and more complete at the hottest concentrations of ionized electrons, 
where the baryon densities are the greatest. 

Unlike the high precision, high resolution data analysis for the Bullet Cluster 1E0657- 
558 of §5.3, the X-ray mass profiles of Figure 5.1 are based on the high precision, but low 
resolution, data analysis of Reiprich (2001); Reiprich and Bohringer (2002), and therefore 
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updating the results with Chandra data is an important future effort, particularly as the 
next generation of high resolution Sunyaev-Zeldovich imaging also becomes available. 

8.2.2 Gravitational lensing 

Weak and strong gravitation lensing surveys of large clusters, such as presented in §5.3 
for the Bullet Cluster 1E0657-558, may be used to map the scaled surface mass density 
by means of the convergence ft-map, discussed in §5.3.2. The November 15, 2006 data 
release (Clowe et al : 2006b), shown in Figure 5.7, is a computational reconstruction based 
an a survey of both strong and weak gravitational lensing, using images provided by the 
Advanced Camera for Surveys aboard the Hubble Space Telescope, providing a much more 
complete convergence ft-map than was previously known. Gravitational lensing surveys, 
particularly the CASTLES and COSMOS surveys discussed in §8.2.2, have provided cat- 
alogues of lensing events, which may lead to as yet undiscovered distributions of mass, 
whether it be non-baryonic dark matter or the massive fifth-force fields of a modified 
gravity theory. 

The solutions to the missing mass problem in clusters of galaxies are testable and falsi- 
fiable by current observations, as in Chapter 5, and through increasingly complete, more 
precise, and higher resolution measurements of the intercluster medium will lead to a bet- 
ter identification of the dynamic mass measured by gravitational lensing reconstructions. 
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